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ABSTRACT 


The  object  of  this  thesis  is  to  survey  and  collate  the  extensive 
literature  relating  to  axially  symmetric  vacuum  solutions  of  Einstein's 
gravitational  field  equations,  and  to  develop  a  number  of  new  results. 

"Space-time"  considerations  have  imposed  fairly  severe  restrictions 
on  the  range  of  topics  covered.  For  instance,  the  large  amount  of  work 
on  approximate  solutions, on  axially  symmetric  gravitational  waves  and  on 
solutions  of  combined  Einstein-Maxwell  field  equations  has  been  left 
entirely  out  of  account.  Within  these  limitations  however,  the  author 
has  aimed  to  be  reasonably  comprehensive. 

We  itemize  briefly  the  main  original  contributions  in  the  thesis. 

§  2.5.  The  theorem  of  page  J2  is  stated  here  for  the 
first  time.  While  a  result  of  this  type  is  implicit 
in  the  work  of  several  earlier  authors,  it  has  not 
previously  been  given  a  precise  and  general  formulation. 

§  page  50  and  Appendix  IV,  page  118,  the  fields 

of  oblate  and  prolate  spheroids  given  here  are  based 
on  the  original  calculations  of  the  author. 

§§  3*4,  3*5>  3*6.  Field  of  a  set  of  collinear  particles. 

§  4.3.  Isolation  of  a  new  class  of  axially  symmetric 


stationary  fields. 
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Chapter  V.  Derivation  of  an  explicit  solution  of 
Bondi's  problem. 

A  number  of  other  new  results  are  scattered  through  the  thesis. 
Where  these  involve  novelty  of  substance  rather  than  merely  of  presentation 


they  have  been  explicitly  indicated  as  such. 
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CHAPTER  I 


HISTORICAL  REVIEW 


This  chapter  is  intended  to  review  the  work  done  by  authors  on 
vacuum  gravitational  fields  which  possess  axial  symmetry  in  general 
relativity. 

Such  a  review  of  work  is  appropriate  at  the  present  time  because 
the  renaissance  of  interest  in  general  relativity  is  yielding  a  flood  of 
new  solutions,  by  using  recently  developed  techniques  (e.g.  algebraically 
special  fields).  At  the  same  time,  a  considerable  body  of  solutions  were 
discovered  and  rediscovered  by  many  authors  in  the  1920's  and  1930* s,  and 
appeared  in  numerous,  but  not  easily  accessible  journals.  One  of  the  purposes 
which  this  thesis  will  serve,  will  be  to  collect  some  of  these  results  under 
one  roof. 

In  the  extraordinary  burst  of  activity  in  Germany  and  other  countries 
which  followed  publication  of  Einstein's  gravitational  theory,  the  theory  of 
static  axially  symmetric  fields  was  one  of  the  earliest  and  most  important 
achievements.  Einstein's  definitive  formulation  of  his  theory  of  gravitation 
appeared  in  the  "Sitzungsberichte  der  Preussischen  Akademie  der  Wissenschaften" 
in  late  1915  [1].  Just  two  or  three  months  later  the  same  journal  carried 
the  first  exact  non-trivial  solution  of  the  fields  equations  -  Karl 
Schwarzschild ' s  spherical  symmetric  solution  representing  the  static  exterior 
field  of  a  massive  sphere  or  a  particle  [2], 

In  a  classic  paper  published  in  1917>  Weyl  [3]  showed  that  by 
introducing  suitable  coordinates  -  "canonical  cylindrical  coordinates" 
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r,  z,  Cp  - 
functions 


?  could  be  reduced  to  the  determina 
A 

X(r,  z),  v(r,  z)  ,  the  first  satisfying  Laplace's  equation 


the  vacuum  field  equations  could  be  reduced  to  the  determination  of  two 


a2x  i  ax  a2x 

- -  +  -  v~  +  - 7;  =  0 

3r2  r  57  3z2 


and  the  second  a  Poisson  equation. 

Levi-Civita  [4]  working  at  almost  the  same  time,  further  simplified 
the  theory  of  static,  axially  symmetric  gravitational  fields.  He  showed  that 
v  is  obtainable  from  X  by  a  straight  quadrature. 

In  this  way,  a  large  class  of  exact  solutions  of  the  field  equations 
was  made  available.  Some  idea  of  Weyl's  mathematical  power  and  of  the 
difficulties  of  the  problem  can  be  gained  from  the  fact  that  no  significant 
widening  of  this  class  has  been  achieved  in  the  almost  half-century  since 
the  appearance  of  Weyl's  paper.  What  little  is  known  about  exact  non¬ 
static  solutions  has  often  been  won  by  using  some  simple  variant  of  the 
Weyl-Levi-Civita  formalism.  For  instance,  Einstein  and  Rosen  [5]  and  later 
workers  [6]  applied  the  complex  transformation  z  -»  it  ,  t  — >  iz  to  Weyl's 
line-element  to  discuss  cylindrical  gravitational  waves.  In  this  category 
belongs  too,  the  essay  of  Bondi  [7]  on  the  motion  of  two  masses  +  m  ,  -  m 
discussed  in  detail  in  Chapter  V. 

In  1922,  Weyl's  pupil  R.  Bach  [8]  used  his  techniques  to  derive  a 
number  of  explicit  axially  symmetric  solutions,  in  particular  the  field  of 
a  circular  ring  and  the  static  field  of  two  particles.  Already  in  1917> 

Weyl  had  indicated  that  the  Schwarzschild  solution  is  obtained  if  one 
formally  puts  X  proportional  to  the  Newtonian  potential  of  a  uniform  rod. 
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To  obtain  his  2-particle  solution,  Bach  chose  A  equal  to  the  potential  of 
two  such  rods  placed  along  the  z-axis. 

In  a  supplement  to  Bach's  paper,  Weyl  [9]  pointed  out  a  pit-fall 
into  which  duly  fell  a  succession  of  later  authors  [ 10 — 1 3 ]  who,  unaware  of 
this  work,  independently  rediscovered  Bach's  formal  result  or  variants  of 
it.  It  is  possible  that  the  singularity  of  the  cylindrical  co-ordinates 
on  the  axis  of  symmetry  r  =  0  may  mask  a  true  singularity  of  the  space- 
time  manifold  on  this  axis.  By  paying  due  attention  to  this  point,  Bach  and 
Weyl  were  able  to  show  that  there  is  in  fact  a  singularity  of  the  manifold 
along  the  line  joining  the  two  particles.  This  singularity  is  nowadays 
interpreted  as  a  strut  which  holds  the  two  particles  apart,  and  keeps  them 
at  rest. 

In  his  paper  [9],  Weyl  established  the  general  result  that  Einstein's 
field  equations  admit  no  axially  symmetric  solution  representing  two  static 
bodies  completely  separated  by  a  surface  z  =  constant,  which  is  non-singular 
everywhere  outside  the  two  bodies.  This  was  one  of  the  first  anticipatory 
clues  to  the  most  remarkable  and  significant  feature  of  Einstein's  theory, 
a  feature  which  still  distinguishes  it  from  all  other  present-day  field 
theories.  The  non-linearity  of  the  field  equations  means  that  i  f  <px  .  <p2 
are  solutions  each  representing  a  static  singularity,  it  is  impossible  to 
superimpose  these  to  get  a  static  2-body  solution  cp^  +  cp^  •  A  solution 
having  two  singularities  exists  only  under  special  circumstances  -  provided 
the  two  singularities  move  in  the  right  way!  The  equations  of  motion  are 
in  fact  built  into  the  field  equations.  There  is,  for  instance,  no  need  to 
assume  as  a  separate  postulate  that  test-particles  follow  geodesics  in 
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space-time,  as  Einstein  did  in  1915  -  this  is  a  consequence  of  the  field 
equations  [ 1 6 ] . 

However,  it  was  left  to  Einstein  in  1927  [17]  to  extract  the  full 
physical  implications  of  Weyl ' s  result  and  with  his  collaborators,  to  develop 
the  detailed  theory  of  gravitational  motion  [18].  In  1922,  Weyl  merely  stated 

[9,  P.145] 


"Die  physikalische  Bedeutung  dieses  Resultats  wird  man  nicht 
ubertreiben  durfenj  fur  die  Losung  des  wirklichen  Zweikorper 
-  Problems,  die  Bestimmung  der  Bewegung  zweier  sich  anzie- 
hender  schwerer  Massen,  ist  damit  nichts  gewonnen.  Aber  es 
ist  doch  ein  physikalisch  sinnvoller  exakter  Ausdruck  gewonnen 
fur  die  Kraft  mit  der  sich  zwei  Massenpunkte  nach  der 
Einsteinschen  Theorie  anziehen" . 


Perhaps  because  of  post-war  communication  difficulties,  the  work 


of  Bach  and  Weyl,  and  even  the  earlier  work  of  Weyl  and  Levi-Civita,  remained 
for  some  time  largely  unknown  outside  the  countries  of  their  authors.  In 
particular,  the  work  of  most  English  writers  seems  to  have  been  completely 
independent.  W.  Wilson  [19]  obtained  in  1920  the  external  field  of  an 
infinite  rod  or  cylinder.  In  1924,  H.  E.  J.  Curzon  [12]  formally  derived  a 
static  solution  representing  two  "particles"  in  which  X  was  given  by 


This  agrees  with  Bach's  solution  only  asymptotically  at  large  distances,  and 
unlike  it,  does  not  reduce  to  Schwarzschild ' s  solution  when  one  of  the  masses 
vanishes.  Curzon  also  failed  to  recognize  that  his  solution  becomes  singular 
on  the  axis  of  symmetry  between  the  particles.  In  1936,  in  a  paper  entitled 
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"Two-centres  solution  of  the  gravitational  field  equations  and  the  need 
for  a  reformed  theory  of  matter"  L.  Silberstein  [ 1 5 ]  rediscovered  Curzon's 
result,  again  overlooking  the  axial  singularity  (c.f.  the  reply  by  Einstein 
and  Rosen  [14]).  In  the  meantime,  Bach's  formal  result  had  been  rediscovered 
by  Chazy  [10]  in  France,  and  Palatini  [11]  in  Italy.  A  rather  obscure  paper 
on  the  same  subject  by  E.  Trefftz  [ 1 5 ]  appeared  almost  simultaneously  with 
Bach's  in  1922. 

Recently,  Marder  [20]  has  given  an  interesting  new  interpretation 
to  Curzon's  solution:  it  may  be  regarded  as  the  exterior  field  of  a  body 
encased  in  a  second  hollow  body.  Bach's  solution  can  of  course  be  similarly 
interpreted  [see  Chapter  III]. 

The  extension  of  Curzon's  solution  to  the  static  field  of  n 
"collinear  Curzon  particles"  was  first  explicitly  carried  out  quite  recently 
by  Hoffmann  [21],  It  was  shown  by  Bergmann  [22]  that  such  a  solution 
necessarily  involves  extraneous  singularities  if  the  masses  are  all  assumed 
positive.  This  corresponds  to  the  fact  that  there  is  no  Newtonian  equilibrium 
configuration  for  n  collinear  positive  masses.  However,  Hoffmann  pointed 
out  that  static  solutions  without  struts  exist  if  the  particles  are  allowed 
to  have  masses  of  either  sign.  Similar  remarks  have  been  made  by  Synge  [23]. 
Even  in  this  case,  it  seems  that  no  strut-free  static  solution  exists  if  n 
is  even.  The  corresponding  generalization  of  Bach's  solution  to  n  collinear 
"Schwarzschild  particles"  is  given  for  the  first  time  in  this  thesis  (Chapter 
III). 


Closely  related  to  this  is  a  problem  first  considered  by  Bondi  [7] 
who  has  attempted  to  examine  the  dynamics  of  a  pair  of  equal  and  opposite 
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masses.  Simple  arguments  based  on  the  principle  of  equivalence  lead  to 
the  expectation  that  the  two  masses  if  started  with  the  same  velocity 
will  move  as  a  rigid,  uniformly  accelerated  system*-  the  acceleration  being 
directed  from  the  negative  to  the  positive  mass.  He  obtained  a  vacuum 
line-element  whose  domain  of  validity  does  not  cover  the  whole  space-time, 
but  only  covers  a  sector,  one  quarter  of  the  whole  space-time  bounded  by 
light-like  asymptotes.  He  was  able  to  prove  the  global  existence  of 
another  type  of  solution  involving  four  particles.  His  solution  possesses 
reflexional  symmetry  and  represents  two  identical  Bondi  dipoles  moving  in 
radially  opposite  directions.  We  discuss  this  problem  in  detail  in  Chapter 
V. 


Very  recently  Bonnor  and  Swaminarayan  [24]  have  presented  exact 
global  solutions  of  Einstein's  equations  representing  four  Curzon-type 
particles  in  motion.  One  of  their  metrics  is  an  explicit  global  solution 
of  Bondi's  4-body  problem  for  the  case  of  four  Curzon  particles  (c.f.  our 
simpler  solution  for  four  Schwarzschild  particles).  They  have  also  presented 
solutions  corresponding  to  the  case  where  all  four  masses  are  positive.  In 
this  case,  struts  along  the  z-axis  are  necessary  to  maintain  the  motion. 

On  the  basis  of  Weyl-Levi-Civita  formalism,  Chou  [25]  investigated 
the  static  axially  symmetric  fields  of  prolate  and  oblate  spheroidal  homeoids. 

Erez  and  Rosen  [26]  have  given  the  static  field  of  a  particle 
possessing  a  mass  quadruple  moment. 

In  i960,  Misra  [27]  using  Weyl's  technique,  determined  the  external 
gravitational  field  of  an  oblate  spheroid. 


*  -  called  a  Bondi  dipole  in  this  thesis. 
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Unaware  of  Misra's  work,  the  present  author  carried  out  calcula¬ 
tions  for  oblate  and  prolate  spheroids.  These  calculations  are  presented  in 
Appendix  IV.  The  final  result  agrees  with  Misra's  work  in  the  case  of 
oblate  spheroids. 

We  pass  from  static  solutions  to  stationary  solutions  with  axial 
symmetry,  corresponding  physically  to  the  field  of  an  axisymmetric  body  in 
steady  rotation.  Approximation  methods  for  dealing  with  this  case  have 
received  considerable  attention.  The  first  approximate  solution  was  given 
by  Lense  and  Thirring  [28].  Bach  [29],  Andress  [30],  Akeley  [31],  Das, 
Florides  and  Synge  [32]  and  Florides  and  Synge  [33]  have  also  used  the 
approximation  method.  However,  this  lies  outside  the  scope  of  this  thesis. 

The  field  equations  for  stationary  axially  symmetric  vacuum  fields 
emerge  as  a  complicated  set  of  non-linear  partial  differential  equations. 

As  regards  exact  solutions,  it  has  not  yet  proved  feasible  to  obtain  more 
than  a  few  particular  results.  Lewis  [34]  found  a  solution  with  cylindrical 
symmetry  which  he  interpreted  as  the  external  field  of  an  infinite  rotating 
cylinder . 


A  corresponding  interior  solution  was  derived  by  Van  Stockum  [35]* 
Some  other  particular  solutions  have  been  given  by  Papapetrou  [3^1,  Tiwari 
and  Misra  [37],  but  they  do  not  seem  amenable  to  a  realistic  physical 
interpretation.  In  general,  it  is  fair  to  say  that  the  approach  of  Andress, 
Lewis  and  their  successors  has  not  to  date  yielded  anything  of  real  physical 
value . 


Very  recently,  Roy  Kerr  [38]  investigating  the  class  of  empty  space 
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times  which  are  characterized  by  the  existence  of  a  geodesic,  shear-free 
null  congruence,  obtained  a  remarkable  particular  solution  which  is 
stationary,  possesses  axial  symmetry  and  contains  the  Schwarzschild  solution 
as  a  special  case.  He  has  suggested  that  his  solution  represents  the  field 
of  a  spinning  particle.  Still  more  recently,  Newman  and  Janis  [39]  have 
put  forward  arguments  supporting  their  interpretation  of  Kerr's  solution 
as  the  field  of  a  uniform  circular  ring  spinning  about  its  axis. 
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CHAPTER  II 


AXIALLY  SYMMETRIC  VACUUM  FIELDS 
INDEPENDENT  OF  TIME:  GENERAL  THEORY 


§2.1  Einstein's  theory  of  gravitation. 

In  its  mathematical  aspects,  general  relativity  is  the  study  of 

a  four-dimensional  Riemannian  space  ("space-time").  In  terms  of  an 

admissible  co-ordinate  system,  the  squared  metric  interval  between  two 

i  j.  i 

neighbouring  points  or  "events"  x  and  x  +  dx  is  given  by  the  quadratic 
differential  form 


$  =  g.  .(x)  dx1  dxJ 
ij 

with  repeated  latin  indices  summed  from  1  to  4.  The  signature  of  the  metric 
is  assumed  to  be  +2  ,  i.e.  g  is  locally  reducible  at  any  event  to  the 
diagonal  (l,  1,  1,  -1).  The  separation  dx1  of  two  events  is  called 
space-lilce,  time-like  or  null  according  as  $  is  positive,  negative  or 
zero . 


The  contravariant  metric  tensor  gX~*  is  defined  as 


(2.1) 


ij 


where  G(i,j)  is  the  cofactor  of  g^  in  the  determinant  g  =  |g„  |  ^  0 
The  two  kinds  of  Christoffel  symbols  formed  from  the  fundamental 


tensor  g. .  are  defined  as 
ij 
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(2>£)  rij,k  =  i(ajglk  +  Si8jk  -  dk8ij)  >  uhore  aj=^J 


and 


r*1  =  ghk  r 

ij  8  ij.k 


From  these  symbols  we  define  a  fourth  order  tensor 


(2.3)  Rh  =  d.rh,  -  d  r™ -  rj1  r™. 

.ijk  j  lk  k  ij  jm  ik  km  ij 


called  the  "Riemann-Christof fel  tensor". 


By  contraction  of  this  tensor,  we  have 


(2.4)  R. .  =  Rh,h.  =  \rh .  -  c).rh  +  r*1.  -  rh.  r™ 

ij  *ihj  h  xj  j  lh  hm  xj  mj  lh 


This  symmetric  tensor  R^  ,  which  plays  a  fundamental  role  in 
relativistic  gravitational  theory,  is  known  as  the  Ricci  tensor. 


Contracting  the  Ricci  tensor  we  obtain  the  invariant 


(2.5) 


R  =  R*  =  g1JR.  . 

1  1  j 


called  the  "curvature  invariant". 


The  Einstein  tensor  is  defined  by 


(2.6) 


G. .  =  R. .  -  ig. .R  . 
ij  ij  2  ij 


Tile  conditions 
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(2.7) 


0 


are  necessary  and  sufficient  for  flatness  of  the  space-time  manifold,  i.e. 
reducibility  of  the  g^  everywhere  to  diag  (l,  1,  1,  -l). 

For  the  physical  interpretation  of  space-time,  we  say  that  a 
physical  "observer"  or  a  "frame  of  reference"  is  identified  with  an  orthogonal 
triad  of  space-like  vectors  propagated  by  Fermi  transport  [23]  along  a  time¬ 
like  curve  -  the  "world  line"  or  the  "historyMof  the  observer.  The  observer 
correlates  arc-length  along  his  world-line  with  time  registered  by  his  clocks 
and  the  orthogonal  triad  with  a  set  of  fixed  axes  in  his  laboratory  to  which 
he  refers  incoming  and  outgoing  light-signals,  test-particles,  etc.  If 
the  separation  of  two  neghbouring  events  is  time-like,  then  both  "happen" 
to  some  observer  (i.e.  his  world-line  encounters  both  the  events)  and  the 
metric  interval  between  them  is  his  measure  of  the  time-lapse  between  their 
occurrence.  If  their  separation  is  space-like,  then  observers  can  be 
found  for  whom  they  are  simultaneous,  and  the  metric  interval  is  the  spatial 
distance  (as  measured  by  such  observers)  between  their  places  of  occurrence. 


The  dynamics  of  particles  and  photons  in  a  gravitational  field  is 
governed  by  the  geodesic  hypothesis,  actually  a  consequence  of  the  field 
equations  (2.10):  the  world-lines  of  spinless  particles  subject  only  to 
gravitational  influence  are  time-like  geodesics.  The  history  of  a  light 
flash  is  a  null  geodesic. 

To  complete  the  theory,  it  is  necessary  to  specify  how  the 
geometrical  structure  of  space-time  is  determined  by  the  distribution  of 
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matter . 

The  fie ld-equat ions  of  general  relativity  (analogous  to  Poisson's 
2 

equation  ^  V  =  4jtp  in  Newton's  theory)  are 

(2.8)  G. .  =  k  T.  .  , 

ij  ij 


where  k  is  a  contant,  and  T.  .  is  the  energy  tensor  of  the  material 

ij 

generating  the  gravitational  field.  For  a  peifjct.  fluid  with  density  p  , 


pressur  P  ,  and  the  4-velocity  u^  ,  has  the  form 


(2.9) 


T.  . 
ij 


(p  +  P)  u  .u  . 

i  J 


+ 


We  shall  for  the  most  part  confine  our  attention  to  solutions 
valid  in  regions  free  from  matter.  This  means  we  have  to  solve  the  vacuum 
field  equations  under  certain  suitable  initial/boundary/junction  conditions. 
One  of  the  conditions  is  normally  taken  to  be  asymptotic  flatness  at  spatial 
infinity . 


The  vacuum  field  equations  are 

(2.10)  G. .  =  0 

or,  equivalently, 

(2.11)  R. .  =  0  . 

This  is  a  set  of  ten  partial  differential  equations  for  the  ten  g^  , 
quasi-linear  in  the  second  derivatives  and  quadratic  in  the  first  derivatives. 
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The  existence  of  four  differential  identities 


(2.12) 


SjG 


ij  _m  j  i  im  1 

J  +  G  JT  ,  +  G  r, 
mj  mj 


0 


between  the  equations  (2,10)  means  that  four  arbitrary  functions  will 


appear  in  the  solution  for  g_  .  This  is  precisely  the  amount  of  arbitrari¬ 
ness  to  be  expected  on  account  of  the  arbitrariness  of  the  coordinate  system. 


To  sum  up 8  In  general  relativity,  a  gravitational  field  in  empty 
space  is  characterized  by  a  normal  hyperbolic  Riemannian  space  of  four 
dimensions  with  vanishing  Ricci  tensor. 


§2.2  Symmetries  and  groups  of  motions. 

(For  a  detailed  discussion  of  the  material  of  this  section,  the 
reader  is  advised  to  see  [40],  [41],  [42]  in  the  Bibliography.) 

Suppose  that  for  each  point  P(x^)  of  space-time,  there  is  a 
corresponding  point  P(x^)  which  is  physically  equivalent  to  it  in  the 
sense  that  all  intrinsic  properties  of  the  world  as  viewed  by  an  observer 
at  P  are  indistinguishable  from  those  as  viewed  by  an  observer  at  P  . 
Then  the  geometrical  structure  of  space-time  is  unchanged  by  going  from 
P  to  P  ,  i.e.  the  point  transformation 

(2.13)  =  ^(x) 


is  a  rigid  motion  of  space-times 


i.e. 


®(P,  Q)  =  ®(P,  Q) 
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for  two  neighbouring  points  P(x*),  Q(x^  +  dx^),  or 
(2.14)  g^(x)  dx1  d*~'  *  Sj  j (x)  dx1  dx^ 


(2-13), 

(2.15) 


Let  x  =  <pM(x)  denote  the  inverse  of  the  point  transformation 
and  introduce  the  new  co-ordinates 


M 


=  /(x) 


Then  P  has  the  same  numerical  co-ordinate  values  in  the  primed  system 

i 

(x^  =  x^)  as  P  has  in  the  original  system  of  co-ordinates  (i,e, 

co-ordinates  "dragged  along"  with  the  motion). 


(2.16) 


Since  0  is  invariant  under  co-ordinate  transformations, 

o  *  c  8  g 

g  (x)  dx1  dx~*  =  g!.(x')  dx1  dx~* 

lj  lj' 

=  g! .(x)  dx1  dxj  . 
ij 


Here,  g|  (x)  means  %  compute  g^.(x') 
of  transformation 


from  (2.15)  and  the  tensor  law 


(2.17) 


gab(x) 


and  then  evaluate  the  left  hand  side  at  x1^  =  x1'*  .  Comparing  (2.14)  and 
(2.16),  we  have  as  the  condition  for  a  rigid  motion 


15 


(2.19)  [Bu(x,)1x  *  sljW 


Hence  space-time  admits  a  rigid  motion  if  and  only  if  the  following  condition 
is  satisfied: 


There  exists  a  co-ordinate  transformation  =  cp^(x)  such 

I 

that  g| ^(x1)  is  the  same  function  of  x*3  as  g^ ^(x)  is  of  x^ 


If  we  now  assume  that  P 

consisting  of  points  equivalent  to 

g.,(x)  dx*  dx^  =  g,.(x)  dx^-  dx^ 

ij  ij 

tion 


can  be  reached  from  P  by  a  path 
P  ,  then  the  condition 
will  be  satisfied  for  the  point-transforma- 


(2.20)  =  /(x,  cr) 

containing  a  continuous  parameter  cr  ,  for  some  value  of  which  (say  cr  =  0), 

(2.20)  becomes  the  identity  transformation  x^  =  x^  . 

The  set  of  transformations  (2.20)  form  a  group,  for  it  is  clear 
from  (2.14)  that  if  P^  is  equivalent  to  P  and  P  is  equivalent  to  P  , 
then  P^  is  equivalent  to  P  .  For  the  infinitesimal  transformation  of  the 
group 

(2.21)  x11  =  +  ^(x)  5cr  , 


where 


- 

' 


5-  .. 
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(2.14)  yields 
(2.22) 


i.e. 

(2.23) 


i  (dx^+df-  ^60")  (dx^+di- ^5cr) 

=  g^x)  dx1  dxJ 


a.g.. 

k  ij 


+  +  Sjk3/ 


=  0 


which  can  be  written  in  the  form 


(2.24) 


?i/j  +  V 


0 


These  are  Killing's  equations  for  the  infinitesimal  generator  £  of  a 
continuous  group  of  rigid  motions.  As  is  clear  from  (2.23),  they  also 
impose  restrictions  on  the  geometric  structure  of  space -time . 

Let  the  four  functions  cp~*(x)  be  independent  solutions  of 

Sl(x)  ^4  =  0  (j  =  1,  2,  3) 

ax1 

=  i  (j  =  k) 


and  make  the  co-ordinate  transformation  x^  =  cp~*(x)  .  In  the  new  co¬ 

ordinate  system  we  have 


Hence  (see  [4.1])  : 
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Co-ordinates  can  always  be  found  which  make  the  components  of 
the  infinitesimal  generator  of  a  one-parameter  group  reduce  to 

S1  =  I2  =  I3  =  o  ,  =  1  . 

In  these  co-ordinates,  the  finite  equations  of  the  group  are  [c.f.  (2.21)] 

(2.25)  x  ^  =  x4  (|a  =  1,  2,  3),  x1'  =  +  a  , 

and  the  Killing  equations  (2,23)  reduce  to 

(2.26)  =  0  , 

a  result  immediately  obvious  from  (2.14)  and  (2.25). 


§2.3  Stationary  and  static  space-time  with  axial  symmtery. 

A  space-time  manifold  will  be  called  stationary  if  it  admits  a 

one-parameter  group  of  time -like  rigid  motions,  that  is,  there  exists  a 

vector  satisfying  (2.23)  and  <  0  .  According  to  the  last- 

result  of  §2.2  ,  this  means  that  there  exists  a  co-ordinate  system 
12  3  4 

x  ,  x  ,  x  ,  x  =  t  ,  with  t  a  time-like  co-ordinate,  such  that  the  metric 

12  3 

tensor  g_(x  ,  x  ,  x^)  is  independent  of  t  .  This  co-ordinate  system 
is  fixed  to  within  an  arbitrary  "spatial  transformations" 


n  Ji,  1  2  3n 

<r  =  cp  (x  ,  x  ,  x  ) 


where  (a  =  1,  2,  3  and  to  arbitrary  time  translations,  t 1  =  t  +  constant  , 
since  these  preserve  the  relations 
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S1  =  (o,  o,  o,  l)  , 

x1'*  =  x^fii  =  1,  2,  3),  it1'  =  xlh  +  o r 

Vij  ■  0  • 


A  stationary  space-time  is  called  axially  symmetric  if 


who 


(i)  there  exists  one-parameter  group  of  space-like  rigid  rotations 
se  generating  vector  is  orthogonal  to  ^  ; 


(2.27) 


>  0  ,  S^V  =  0 


i  J 


\/k  +  v 


+  T)  =  0 


The  orthogonality  of  I;1  ,  rj1  implies  that  r]1  defines  a  rigid  motion  in 
each  hypersurface  t  =  constant.  Hence  by  a  spatial  transformation  we  can 

1  g  -z  12 

introduce  co-ordinates  x  ,  x  ,  xJ  =  cp  such  that  g_(x  ,  x  )  is 
independent  of  cp  . 


reversal  o 


(2.28) 


(ii)  Space-time  geometry  is  invariant  under  the  simultaneous 

f  cp  and  t  ; 

cp  — >  cp  =  -  cp  ,  t-*t  =  -  t  . 


According  to  (2.14)  this  gives 


v.;  \r A 


■ 

p  To  V'  3'.  cik 

f" 
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(2.29) 


^ABdxAdxB  "  2SA3dxAdcP  “  28A4dxA(lt  +  8^dC(/  +  2g  ^dcpdt  +  §44dt? 
=  8ABdxAdxC  +  2gA^dxAd9  +  2g^dxAdt  +  g^dcp2  +  2g^dcpdt  +  g^dt2 


inhere  ^  =  g^  (-cp,  -t)  and  (A,B  =  1,2)  .  But  since  are  independent 

of  (p  and  t  ,  we  have  ,  and  we  infer  from  (2.29)  that 


(2.30) 


8^4  -  Sol.  -  8i  -z  ~  Sox  _  0 


>2k 


13  &23 


!7e  briefly  indicate  the  intuitive  motivation  for  the  above 
definition.  A  stationary  axially  symmetric  gravitational  field  is  the  kind 
oi  field  produced  by  an  axially  symmetric  distribution  of  matter  rotating 
steadily  about  its  axis  of  symmetry.  An  electromagnetic  analogue  would  be 
the  field  of  an  axially  symmetric  charge  distribution  in  steady  rotation. 
Imagine  a  compass  inserted  to  show  the  direction  of  the  magnetic  field,  and  a 
film  taken  of  the  phenomenon.  The  film,  if  run  backwards  (t  -  t)  will 
net  give  an  accurate  picture  of  the  field  produced  by  the  same  distribution 

rotating  in  the  opposite  sense,  since  the  magnetic  field  ought  to  change 

direction  but  does  not.  However,  if  we  somehow  reverse  the  sense  of 
rotation  (cp  — >  —  cp)  while  running  the  film  backwards,  we  return  to  an 
accurate  picture  of  the  original  field.  Condition  (ii),  thus  gives 
expression  to  the  hypothesis  that  a  stationary  axially  symmetric  gravitational 
field  is  indifferent  to  time  inversion  if  we  at  the  same  time  reverse  the 

sense  of  rotation  of  the  masses  producing  the  field.  It  is  well  to  add 

that  while  this  (or  some  equivalent)  hypothesis  has  usually  been  made  in 


/ 

■ 

t 


. 

■ 


:  1  ,  :  :  .  ■■■■  1 
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discussing  such  fields,  the  case  for  it  is  not  overwhelmingly  strong. 


Summing  up,  in  a  stationary  axially  symmetric  space-time,  it  is 

12  3  4 

possible  to  introduce  co-ordinates  x  ,  x  ,  x  =  cp  ,  x  =  t  in  which  the 

metric  reduces  to 

(2.31)  ®  =  gAfidxAdxB  +  2g  ^dq>dt  +  g^dcp2  +  gUidt2 


’44 


where  (A,  B  =  1,  2)  , 


1  2, 


A,  B 


8ij=  gij(*  >  x  )’  SU4  <  0  >  s33  £  0  >  gABdX“dX”  >  0 


A  further  simplification  is  possible.  We  recall  that  any  two- 
dimensional  line -element 


/  1  2x  A  B 
§AB(X  *  x  ’  dx  dx 


is  reducible  to  "isothermal  form"  [A(x^  ,  x2  )]2[(dx^  )2+(dx2  )]  by  a 
co-ordinate  transformation.  This  is  immediately  seen.  Consider  the  two 
families  of  (complex)  null-curves 


f  1  2x  _  ,  /  1  2N  _  . 

u(x  ,x  )  =  constant,  v(x  ,x  )  =  constant  , 


whose  infinitesimal  tangent  vectors  satisfy 


,  A  .  B  . 
gABdx  dX  =  0  » 


v  is  the  complex  conjugate  of  u  .  In  terms  of  the  null  co-ordinates 
u,  v  the  metric  is 


gAgdxA  dxB  =  [A(u, v) ]2  du  dv 


r  ■  Qj  .  :  v  1.2  t:.r  , 

■  -  '  ■  ■  -  :.'.s  .  :  i-  ■:  .  tfj 

(T\ 

' 


‘ 


,-V:-  ■; 


- 


. 

'■  :  :  o:  ,a  p  ■ 


.  'fCi-  "  r  .  .  .. 

-  -  '  ■  '  -  r'sysyji-.^j  .i.v.v:  ■ :  .1  ;  f  -  ..  ,!?* 
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Defining  the  real  co-ordinates 


1  i/\  2  1 /  \ 

c  =  £(u+v),  x  =  2i(«-v) 


the  metric  becomes 


,2r/,  1\2 


2'  n2. 


A  [(dx"  )  +  (d*  )  ] 


Thus  we  have  succeeded  in  reducing  the  metric  tensor  for  a  stationary 
axially  symmetric  field  to  a  form  in  which  the  only  non-vanishing  off- 
diagonal  element  is  g  ^  . 

We  shall  now  make  a  further  specialization; 

Definition;  -  A  space-time  is  called  STATIC  if  it  is  STATIONARY 
and  invariant  under  time-inversion  t  =  -  t  . 


Intuitively,  a  static  field  is  produced  by  matter  in  a  state  of 


rest . 


By  an  argument  similar  to  that  leading  to  (2.30),  we  find  that 
the  static  condition  is  equivalent  to 

(2.32)  Sl)+  =  =  8^4  =  0 

This  means,  geometrically,  that  the  t-lines  are  orthogonal  to  the  isometric 
space-like  hypersurfaces  t  =  constant. 

We  have  now  shown  that  in  a  static,  axially  symmetric  space-time 

12  3  4 

co-ordinates  x  ,  x  t  x  =  cp  9  x  =  t  can  be  found  which  reduce  the  metric 

to  the  form 


^  ,  f,-,.  ’  : 


:  ..  "  ■-  : 
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(2.33)  <&  =  A^fdx1)2  +  (dx2)2]  +  B2(dx5)2  -  C2(dx1')2 

1  2 

where  A,  B,  C  are  functions  of  x  ,  x 


§2.4  Weyl  fields. 

Following  Weyl  [3]  and  Levi-Civita  [4],  we  now  consider  the  solution 
of  Einstein's  field  equations  for  static,  axially  symmetric  vacuum  fields  - 
often  called  Weyl  fields. 

From  (2.33) 

.  .2r  ,  1 \2  2x2,  _2 /  ,  3n2  n2fA  4x2 

<t>  =  A  [dx  )  +  (dx  )  ]  +  B  (dx  )  -  C  (dx  ) 

1  2 

where  A,  B,  C  are  functions  of  x  ,  x 


From  the  components  of  Ricci  tensor  corresponding  to  (2,33)>  which 


we  have  given  in  the  Appendix  I  for  convenience,  it  follows  that 


(2.34) 


Rk  -  ^ 

h  -  b2 


1  f  S2(bc)  S2(bc)  1 

„2„„  1  ,,  1,2  2,2  f 

A  BC  (ox  )  (ox  )  J 


In  a  matter-free  domain,  the  field  equations  are 


(2.35)  R.j  =  0  . 


Hence 


(2.56) 


3  4 

R  i  +  R, 
3  4 


2 

A  BC 


a2(Bc) 

(ax1)2 


+ 


a2(Bc)  1 

(ax2)2  J 


0 


or 


•  -  ■  :;r  ^  ■ 


■ 


*•  faff' :  ::  S  fjT 
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(2.37) 


a.2(Bc)  _a2.(B_c)  _ 

/*  1x2  +  2x2  " 

(dx  )  (dx  ) 


1  2 

i.e.  BC  is  a  harmonic  function  of  (x  ,  x  )  .  Let 


(2.38) 


BC  =  r (x1,  x2) 


and  suppose  z  is  the  conjugate  function  of  r  ,  then 


(2.39) 


•  e  f  1  .  2  \ 

r  +  xz  =  f(x  +  ix  ) 


12  12 
where  f(x  +  ix  )  is  an  analytic  function  in  x  +  ix 


1  2- 


In  making  the  conformal  transformation  from  (x  ,  x  )  to  (r,  z) 


we  have 


(2.40) 


.2r  /  ,  1\2  2\2-i  2r  ,  2  ,  2 

A  [ (dx  )  +  (dx  )  ]  =  a  [dr  +  dz  ] 


where  a  is  a  function  of  (r,  z) 


From  (2.38) 


(2.U1) 


B  = 


Writing 

(2.42)  x^  =  r  ,  x2  =  z  ,  x^  =  cp  and  x^  =  t 
the  metric  form  (2.33)  can  be  written  as 

/o  ,  *  2/ ,  2  ,2,  2*2. _2  n2.  2 

(2.43)  a  (dr  +  dz  ]  +  r  C  dcp.  »  C  dt 

where  A  and  C  are  functions  of  (r,  z)  . 


i'-  "  ;  !  Si- 


n-o.t  u  •  is^vxls^o  $./L  $ijj .!««  jjJ 


-  .}  ■  '  t-  ■:  S';  o/J-r 


' 


“TC  i, 


j, 


*#)■  -S'c:  a  .  V  "  -a  :  ■.  .,  w 
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Hence  we  may  state  that  in  any  static  axially  symmetric  domain, 

■5  4 

to  reduce  the  metric  to  the  form  (2.43),  it  is  sufficient  that  =  0 

In  order  to  avoid  cumbersome  mathematical  manipulations  to 
determine  the  vacuum  equations  both  Weyl  [3]  and  Levi-Civita  [4]  assume  that 

v~A  -A  A 

A  =  e  ,  B  =  re  and  C  =  e 


A  and  v  are  functions  of  (r,  z)  . 


Hence  (2.43)  becomes 


(2.44) 


2(v"A) ,  2  1  2s  2  -2A  2  2A  2 

i  v  (dr  -t-  dz  )  +  r  e  dtp  -  e  dt 


The  expressions  for  the  non^vanishing  components  of  the  Ricci  tensor 
corresponding  to  (2.44)  are 

R11  -  V11  +  V22  ’  Au  -  *22  +  +  -V1 


2  Vl”^] 

r  tz  v  +v  -A  -A  +  2A  -  — — — — • 
22  11  22  11  22  2  r 


2A, 


(2.45) 


R10  =  2  A  A  -  — 

12  1  2  r 


R 


33 


”2v  2  r*  -1 

e  r  [\l  +  7*1  +  A221 


R 


44 


-2v  +  4A,-»  l-x  >  i 
t\l  +  7*1  +  *22] 


where  A,  = 


3A  ,  ^A 

ch:  3  2  Bz 


and  so  on  . 


•; 
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where 


Av  = 


a!v 

3r2 


we  have 

A 

‘  +  — )  + 

■\  2 

A,  + 

II 

OJ  ( 
,< 

1 

2 

t*  ;  m 

A 

A 

Of 

bz* 

or2 

Sz2 

4(Rn  -  R22) 


X2  -  X2 


=  0 


R12=  2V2 


=  0 


R33  ”6  r  ^11  +  r  \  +  A22 


=  0 


3  4 

R  i  +  R,  =  0 

3  4 


i.e . 

(2.46) 

<i 

ro 

—  -  -f-  ^  —  0 

H  r  1  22 

(2.47) 

V1 

■  ^  -  44 

(2.48) 

V2 

=  2r\A2 

(2.49) 

Av 

+  A2  +  A2  *  0  . 

Thus  A  satisfies  the  Laplace's  equation  (2.46)  in  cylindrical 
co-ordinates  (r,  cp,  z)  in  Euclidean  J-space  and  v  satisfies  the  Poisson 
equation . 


. 


. 

• 

■ 
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Hence  by  the  introduction  of  "canonical  cylindrical  co-ordinates 
(r,  z,  cp)",  Weyl  [3]  and  Levi-Civita  [4]  reduced  the  problem  of  obtaining 
any  static, axially,  symmetric  vacuum  field  equations  to  the  determination 
of  two  functions 


A(r,  z)  and  v(r,  z) 

v(r,  z)  is  obtained  from  A  by  a  straight  quadrature 


(2.50) 


v 


r 


dr  +  2' 


<bA 

dr 


dA 

dz 


§2.5  Elementary  flatness. 

We  come  now  to  a  slightly  subtle  point  vHih ,  as  we  have  explained 
in  Chapter  I,  was  overlooked  by  many  early  authors.  Weyl's  canonical 
cylindrical  co-ordinates  become  singular  on  the  axis  of  symmetry  r  =  0  . 

To  discover  whether  this  co-ordinate  singularity  conceals  a  true  singularity 
of  the  manifold  at  a  point  z  =  z^  of  the  axis,  one  can  check  whether  the 
2-space  S  formed  by  drawing  all  geodesics  normal  to  the  z-axis  at 
has  finite  Gaussian  curvature  at  z^  . 

We  recall  that  the  Gaussian  curvature  K  of  a  2-space  at  any 
point  can  be  computed  by  surrounding  the  point  by  a  small  geodesic  circle 
of  circumference  C  and  a  geodesic  radius  p  ;  then  [43] 

(2.51)  K  =  limit 

p  — ►  o  —  rtP 


■5  "  1  '  -  ]  V/,  /  -  V  !  •  ;  '!(  ,  . 

I 


•c.  aif ;  ;-iOi  c  :  ~ '  to 


r 


'■  1  .  ■■■■fir  " 


■  wr.ppt:  •  «ot  iaol*  y~ 

:  7  ■  ta  -  -  :o*os  ai  :i  .  7  ;>? 


it&r  -i  -  ■  i.  J'ftr.v  UV.'S-x  ,7' 


:>  S'  .r/  rtmr/pr.  .  ,  • 
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Thus,  for  the  absence  of  singularities,  it  is  necessary  that 


2jt  -  —  =  0(p  )  p  — ►  0 

P 


for  every  geodesic  circle. 

Intuitively,  (2.5I)  expresses  the  condition  that  the  2-space 
should  possess  a  tangent  plane  at  the  point  p  =  0  . 

In  particular,  for  a  geodesic  circle  in  S  about  r  =  0  ,  z  =  z^  , 
we  have  from  (2.44) 

p  ^  re^V  ^  ;  C  ~  2«re^  (r  -»0)  . 

Hence  the  point  r  =  0,  z  =  is  a  true  singularity  of  the  manifold  (2,44) 
unless 


(2.52) 


v(0,  zx)  = 


0 


Now  a  singularity  of  K  on  the  z-axis  implies  a  singularity  of 

the  Riemann  tensor.  This  implies  the  violation  of  the  vacuum  equations 

R. ,  =  0  on  the  z-axis ,  and  therefore  the  presence  of  a  material  strut  on 
ij 

the  z-axis. 


It  was  Weyl  [9]  who  first  pointed  out  the  remarkable  repercussions 
of  elementary  flatness  condition  on  the  existence  of  axially  symmetric  static 
2-body  solutions  in  general  relativity.  We  shall  discuss  this  in  detail 
with  respect  to  particular  solutions  in  Chapter  III.  The  general  argument  which 
follows  is  a  simplification  and  generalization  of  those  due  to  Weyl  [9] 
and  Bondi  [7 ] . 


■ 


■  '  -1  ■  ii:  t£  9  ■'!  :  ■  {  l  ;  ,  O  ■  I  i 


■ 
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Let  the  metric 


2(v-A)  /  ,  2  ,  2X  2 -2 A  ,  2 

<t>  =  e  '  '(dr  +dz  )  +  r  e  dcp 


2A  .2 

-  e  dt 


represent  the  static,  axially  symmetric  field  of  a  set  of  discrete  bodies 


B,  B^  >  B0,  .... 

.  This 

means  the 

vacuum  field  equations 

(2.53) 

v2a 

=  0 

(2.54) 

dv 

dr 

- 

(2.55) 

dv 

-  2r—  . 

’  dr 

dA 

dz 

are  satisfied  everywhere  except  within  the  shaded  regions  B,  B^  .... 

in  the  Euclidean  map  of  the  cylindrical  co-ordinates  (r,  z,  cp)  (see  Fig.  l) 


Let  us  suppose  that  all  matter  is  contained  within  a  finite 
distance  from  the  origin.  Then  (like  the  Newtonian  potential  which  it 
asymptotically  approaches)  A  behaves  at  spatial  infinity  according  to 


(2.60) 


A (r ,  z)  =  0(i)  (p  ->oo) 


(2.61) 


,b A  6A, 

dr  ’  ciz' 


0(™) 

P 


where 


2  2  2 
p  =  r  +  z 


From  (2.55)  we  have  v  =  constant  on  segments  of  z-axis  which  are  free  of 
matter.  Let  us  suppose  the  arbitrary  constant  in  v  is  adjusted  to  make 


■  a  !■ 


. 

I  -i  A  : 


. 


i 


.  '  ..  .  .  .  ,  v  '  •  ■  .  ;  V  jnu: 


Figure  1. 


Figure  2. 
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v  =  0  on  the  segment  KL  of  z-axis,  extending  down  to  z  =  -  oo  (Fig.  2). 

To  compute  v(r,  z)  at  any  other  point  P(r,  z,  cp),  we  connect 
P  to  a  point  a  on  KL  by  a  path  C-j  and  evaluate  the  line  integral 


(2.62) 


/dA,2  ,dA.,2 

(dr}  " 


dr  + 


dA 

dz 


dz 


Since  (2.55)  is  the  integrabi lity  condition  o  f  (2.54)  and  (2.55), 
it  therefore  expresses  the  condition  that  the  line-integral  is  invariant 
if  the  path  is  deformed  continuously  into  another  path  (such  as 

or  C^)  joining  P  to  KL  without  passing  through  matter  in  the  process 
of  deformation.  This  proves  immediately  that  v  =  0  on  the  segment  RS  . 
For  a  path  joining  KL  to  RS  can  be  deformed  into  an  infinite  semi-circle 
and  (2.61)  ensures  that  the  line-integral  over  this  semi-circle  vanishes. 


To  compute  v  on  the  segment  ST  of  the  axis,  we  connect  KL 
and  ST  by  a  path  C  in  empty  space  lying  in  a  meridian  plane  cp  =  constant 
Then 


(2.65) 


ST 


r(ds  x  A) 


where  ds 


=  (dr,  dz,  0)  is  the  element  of  arc-length  along  C 


(2.64) 


A  .  2^ 

r  dr 


dA 

dz  ’ 


r dA.2  /dAx 2 

^dr  'clz 


Acp  0 


and  the  subscripts  r,  z,  cp  denote  physical  components  in  the  Euclidean 


map . 


h,  ij  o  0  -  -r 

■‘Is 
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Let  S  be  the  surface  of  revolution  formed  by  revolving  C 
about  the  axis.  If  n  is  the  unit  outward  normal,  we  have 

'•o 


(2.65) 


rdcp(ds  x  A) 


A  •  n  dS 


and  hence 


(2 .66) 


v 

ST 


dcp 


(2.67) 


A  •  n  dS 


This  integral  is  independent  of  the  shape  of  the  surface  S  enclosing  the 
body  B  :  this  follows  because 


(2.68) 


8a  .  8a 

, .  .  r  1  .  z 

div  A  =  ■*r— -  +  —  A  +  ’r— — 

~  or  r  r  dz 


v2a 


) 


which  vanishes  in  matter-free  space  by  virtue  of  (2.53)» 


It  is  easy  to  show  that  a  function  A  ,  satisfying  (2,53)>  (2.60) 
and  (2.6l)  in  the  domain  exterior  to  B,  ...  can  be  uniquely 

decomposed  into  a  sum  of  functions 


A  +  A  -f  A  +0.0. 
o  1  2 

2 

such  that  A^  satisfies  7  A,  =  0  everywhere  outside  . 

We  now  construct  a  Newtonian  analogue  system  for  the  bodies  B^ 
by  imagining  each  domain  B^  to  be  filled  with  a  distribution  of  mass 
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whose  classical  Newtonian  potential  is  precisely  V  (r,  z)  exterior  to 


the  region  .  This  mass  distribution  is  not  unique  -  e.g.  any  spherically 

symmetric  distribution  with  a  given  total  mass  has  the  same  external 
potential . 

Converting  the  surface  integral  (2.67)  into  a  volume  integral  taken 
over  the  Newtonian  analogue  of  B  ,  we  have 


v 


ST 


div  A  dv 


(2.69) 


3—  V  A  dv 


dz 


=  -4 


F  p*  dv 
z 


J- 

where  p  =  ^  V  A  is  the  density  of  the  Newtonian  analogue  of  B  ,  and 
Fz  =  -  H  is  the  z-component  of  the  Newtonian  gravitational  force  field. 


We  have  thus  reached  the  following  remarkable  result: 


The  difference  between  the  values  of  v  on  two  free  segments  of 


the  axis  of  symmetry  separated  by  a  material  body  B  is  equal  (apart  from 
a  constant  factor)  to  the  Newtonian  force  exerted  on  any  Newtonian  analogue 
of  body  B  . 


It  follows  that  we  can  satisfy  the  condition  of  elementary  flatness 


on  every  free  portion  of  the  axis  if  and  only  if  we  can  find  a  Newtonian 

analogue  system  in  which  each  body  is  in  dynamical  equilibrium. 
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From  this  result,  and  familiar  results  of  Newtonian  theory  we  can 
infer  immediately  that  there  is  no  static,  axially  symmetric  solution  of 
Einstein's  field  equations  which  corresponds  to  two  bodies  on  opposite  sides 
of  a  plane  z  =  constant,  each  having  non-zero  total  mass.  [The  relativistic 
mass  m^  of  body  is  defined  by 

2A  2m 

(2.70)  -  g44  =  e  =  1  -  ~~  +  0(— g)  (p  -*»).] 

P  p‘‘ 

On  the  other  hand,  static  2-body  solutions  are  possible,  if,  following 
Marder,  we  allow  one  body  to  be  enclosed  by  another. 


Body  1 


Figure  3. 
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CHAPTER  III 

STATIC  AXIALLY  SYMMETRIC  VACUUM  FIELDS: 
PARTICULAR  SOLUTIONS 


§  3*1  Introduction. 

In  the  previous  chapter,  we  found  that  any  static,  axially 
symmetric  gravitational  field  in  vacuo  ("Weyl  field" )  is  characterized  by 
the  canonical  line-element. 

(3.1)  $  =  e  '  '(dr  +dz  )  +  r  e  dcp  -  e  dt  , 


where  the  functions  A(r,z),  v(r,z)  satisfy 


L 


(3.2a) 


hv 


=  r 


(3.2b) 


cv  _  „  5^ 

bz  ~  br  ‘  bz  ’ 


(3.3) 


vh 


b2~K  1  a* 

Sr2  +  7  Sr 


0  . 


(3.2b). 


Equations  (3.3)  is  the  integrability  condition  for  (3.2a)  and 


On  segments  of  the  z-axis  free  of  matter,  it  is  necessary  that  v 


satisfy  the  elementary  flatness  condition 


'  :  -rr  .  •- 

•  . .  ...  ......  .... 
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(3.4)  v(0,z)  =  0  . 

For  the  exterior  fields  of  bounded  sources,  it  is  customary  to 
impose  the  boundary  condition  of  asymptotic  flatness  at  spatial  infinity: 

(3.5)  A  =  0(i)  ,  v  =  0(~)  ,  (p2  s  r2  +  z2  -.»)  . 

P  C 


The  paths  of  infinitesimal  test-particles  in  the  field  are  governed 
by  the  geodesic  equations 


(3.6) 


,2  i  .  ,  p  q 

d  x  _  pi  dx  dx 

,  2  ~  pq  ds  °  ds 


where  the  4-velocity 


is  normalized  by  the  identity 


(3-7) 


g 


dx* 
pq  ds 


dx^ 

ds 


-1  . 


For  particles  moving  slowly  (as  compared  with  the  velocity  of  light)  with 

respect  to  an  observer  at  rest  in  the  static  field,  we  can  neglect  squares 

dx1  dx2  dx5 

and  products  of  the  spatial  components  .  From  (3.6)  and 

(3.7)  we  obtain 


(3.8) 


<-*!*> <3J>2  *  1  ’ 


d2x^  -v  i  /dt\2 

~2~  ~  *  r44  W 

ds 


(3.9) 


-  -islJ 


ai°s(-gU4) 

dx* 
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At  large  distances  from  the  sources,  we  can  make  the  further 
approximations  5*  *  ,  ds  %  dt  ,  by  virtue  of  (3.I)  and  (3.5)  •  Thus, 

the  components  of  acceleration  of  a  slowly  moving  test-particle  at  large 

distances  from  the  sources  of  the  field  (3-1)  are 


(3.10) 


d2r  clA 

dt2  ^  '  * 


d2z  ^  dA 
dt2  ^  '  *Z 


%  0 


The  conditions  under  which  (3.IO)  is  a  valid  approximation  are 
precisely  those  for  which  Newton's  theory  should  give  an  accurate  description. 
We  are  thus  led  to  conclude  that:  at  large  distances  from  the  sources, 
where  the  Newtonian  approximation  is  meaningful 9the  function  A(r,z)  is 

identical  with  the  Newtonian  potential  of  the  distribution.  We  note  further 
the  remarkable  fact  that  in  Weyl's  canonical  co-ordinates,  A(r,z)  satisfies 
throughout  free  space  the  equation  (3«5)s  which  is  formally  Laplace's  equation 

in  cylindrical  co-ordinates. 

These  facts  suggest  the  following  very  convenient  method  of 
generating  and  "labelling"  particular  Weyl  fields. 

(i)  Within  the  framework  of  Newton's  theory,  consider  an  axially 
symmetric  distribution  of  matter  J-jl  having  the  exterior  Newtonian  potential 
V(r, z )  . 

(ii)  Set  A(r,z)  =  V(r,z)  .  Then  (3*3)  is  satisfied.  Compute 
v(r,z)  from  (3.2). 

We  thus  obtain  a  Weyl  field  (3.1)  which  at  large  distances  is 
indistinguishable  from  the  predictions  of  Newton's  theory  concerning  the 


W  ,  ri  ;  :i;  a ;:.i“  t-'.  ;q  ■  b  :«U 

“a  l; 
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external  field  of  the  distribution 


in  Euclidean  space. 


At  smaller  distances,  we  have  to  reckon  with  distortion  of  the 
geometry  caused  by  the  gravitational  field.  Thus,  we  might  set  A 
formally  equal  to  the  Newtonian  potential  of  a  spheroid  and  label  the 
resulting  Weyl  field  as  that  of  a  "spheroid",  but  as  a  description  this  is 
only  good  at  large  distances.  In  a  curved  space,  we  experience  difficulty 
in  even  defining  a  spheroid.  The  point  emerges  most  strikingly  in  the 
case  of  spherical  symmetry.  If  we  start  from  A  =  -  m/(r  +z  )2  , 
corresponding  to  a  spherically  symmetric  Newtonian  distribution  <JD  ,  the 
resulting  line-element  (j.l)  is  not  spherically  symmetric  (see  §  3»3)> 
although  it  approximates  to  spherical  symmetry  at  large  distances.  To 
obtain  the  Schwarzschild  spherically  symmetric  solution  for  a  particle, 
one  has  to  set  A  formally  equal  to  the  potential  of  a  rod ,  as  Weyl  already 
pointed  out  in  his  fundamental  paper  of  191?  [3]«  Each  Weyl  field  generated 
by  steps  (i)  and  (ii)  ought  to  be  re-examined  a  posteriori  to  discover 
what  distortion  of  the  originally  contemplated  Newtonian  distribution  is 
produced  by  the  curvature  of  space.  Such  a  re-examination  involves 
difficulties  both  of  principle  and  of  practice,  and  has  been  carried  out 
only  in  very  few  cases,  such  as  spherical  symmetry.  For  the  most  part  we 
shall  be  content  in  this  chapter  x?ith  the  uncritical  procedure  of  generating 
and  "labelling"  particular  Weyl  fields  according  to  steps  (i)  and  (ii). 

An  initial  assault  on  the  problem  of  intrinsic  characterization  of  Weyl 
fields  is  due  to  Ehlers  and  Kundt  [44,45  Jo 


. 
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§  3*2  Fields  of  unbounded  distributions  of  mass . 

We  begin  our  catalogue  of  special  solutions  with  some  physically 
unrealistic  situations  in  which  the  distribution  of  matter  extends  to 
infinity  and  the  condition  of  asymptotic  flatness  (3*5)  is  violated. 

In  an  attempt  to  reproduce  the  field  of  an  infinite  plane  in 
general  relativity,  it  seems  natural  to  start  from  the  Newtonian  potential 
V  =  gz  ,  where  g  is  a  constant.  Applying  steps  (i)  and  (ii),  we  obtain 
after  an  elementary  calculation 

(3.11)  A  =  gz  ,  v  =  -  ^g2r2  . 

However,  the  Weyl  field  (3.1)  with  A,  v  given  by  (3. 11)  does  not  have  the 
property  of  homogeneity  which  ought  to  characterize  the  field  of  an  infinite 
plane;  i.e.,  it  does  not  admit  a  two-parameter  group  of  spatial  translations 
orthogonal  to  the  axis  of  symmetry.  A  Weyl  field  having  this  property  is 
given  by 

(3.12a)  A  =  ilog£||,  v  =  Jlog 

where 

2  2  2 
p  =  r  +  z 

and  p  is  a  constant.  It  is  easily  verified  that  this  is  a  solution  of 
(3.2),  (3.3).  The  conformal  transformation  of  the  space  co-ordinates 

2  2 

=  r*  z*  >  2$z  =  i(z*  -  r*) 


(3.13) 


2pr 
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converts  the  line-element  (3.1),  (3.12a)  into 


(3.12b) 


*  =  dr l  +  r2  dcp2  +  dz2  -  dt2 


In  this  form,  the  homogeneity  of  the  manifold  is  directly  evident. 


The  spatial  transformation 


(3.14) 


1  2  2 

irS  z*  =  1  +  §z 


-  2 

combined  with  the  scale-transformation  of  time  t  =  —  t  throws  (3.12b) 

§ 

into  the  form 


(3.12c) 


(l+gz)  dt2 


The  gravitational  field  of  an  "infinite  plane"  was  obtained  in  the  form 
(3.12c)  by  Chou  [25]  in  1931j  using  a  somewhat  different  approach,  and  still 
earlier  by  Whittaker  [^6]  who  derived  it  as  a  limit  of  the  Schwarzschild 
solution  \tfhen  the  gravitating  centre  is  removed  to  infinity  and  its  mass 
increased  proportionately  [cf.  Chapter  V  ]. 

Neither  Whittaker  nor  Chou  seemed  aware  that  (3.12)  represents  a 
flat  manifold.  To  see  this,  consider  the  Minkowskian  line-element 


,  ,*2  ,2  ,u2  ,  2 

9  =  dt,  +  dr]  +  d£  -  dT 


and  apply  the  transformation 


|  =  r^  cos  cp  ,  H  =  r*  sin  9 


(3.15) 


■ 


(  .  I . . ) 


. 
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we  then  obtain  (3.12b),  We  observe  from  (3. 15)  that  each  point  with 
fixed  space  co-ordinates  r^  ,  z^  ,  cp  executes  hyperbolic  motion: 

(*  =  constant  ,  r|  =  constant  ,  £  -  t  =  const.  Hence  (3.12b)  may  be 

interpreted  as  the  metric  of  flat  space-time  viewed  from  a  uniformly 

accelerated  reference  frame  or  (according  to  the  principle  of  equivalence) 
as  a  homogeneous  gravitational  field.  The  line -elements  (3.12a)  and 
(3.12c)  are  open  to  the  same  interpretation,  since  they  are  connected  with 
(3.12b)  by  purely  spatial  transformations  x^hich  do  not  affect  the  motion 
of  the  frame. 


We  turn  to  what  could  be  labelled  the  field  of  an  infinite  rod 
or  cylinder : 


(3.16) 


A  =  2m  log 


v  =  4m^  log  — 


A  result  equivalent  to  this  was  given  by  Wilson  [19]  In  1920,  and  before 
him  by  Levi-Civita  [4]  in  1918. 


In  the  case  of  a  hollow  cylindrical  shell  of  radius  a  ,  the 
exterior  metric  (3.16)  can  be  joined  to  the  flat  interior  solution 


2  2  2  2  2 
0  =  dr  +  dz  +  r  d<p  -  dt 


with  continuity  of  the  metric  tensor  g  ,  The  jump  discontinuity  in  the 
first  derivatives, 


'8* 

aVl 

1>J 

r=a 

[8rJ 

(3.17) 


r=a 


1 
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gives,  analogously  to  Newtonian  formula 

(5a8)  ^ -  *  kn<J 

a  measure  of  the  surface  density  of  the  shell.  Comparison  of  (3.I7)  and 
(3.18)  appears  to  indicate  that  m  is  to  be  interpreted  as  the  mass  of 
the  shell  (in  gravitational  units)  per  unit  length.  This  can  also  be 
checked  by  looking  at  the  geodesics  (paths  of  test-particles)  in  the 
manifold  (3.I),  (3.1-6).  For  slowly  moving  particles  [cf.  (3.9)],  we 
have 


(3.19) 


d2cp 

ds2 


0, 


£* 

ds2 


dfi 

ds2 


2m 

r 


where  dr  =  \/g  dr  is  the  element  of  arc-length  along  the  parametric 
lines  of  r  .  Equation  (3. 19)  is  in  agreement  with  the  Newtonian  result 
for  the  acceleration  due  to  an  infinite  cylinder  of  mass  m  per  unit 
length . 


§  3-3  Static  one-body  solution;  Curzon  and  Schwarzschild  particles . 

2 

The  formally  simplest  solution  of  the  potential  equation  V  A  =  0 
which  satisfies  the  asymptotic  boundary  condition  (3.5)  corresponds  to 
the  field  of  a  monopole  in  the  Newtonian  analogues 


4  = 


leading  to  v 


2  2 


m  r 

2p4 


(3.20) 
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2  2  2 

where,  as  usual,  p  =  r  +  z  .  The  manifold  (3.1),  (3-20)  is  not 
spherical  symmetric  [47],  i.e.,  it  is  not  related  to  the  Schwarzschild 
solution  by  a  co-ordinate  transformation.  We  shall  refer  to  it  as  the 
field  of  a  Curzon  particle.  The  curious  nature  of  the  singularity  at 
p  =  0  has  been  remarked  upon  by  Mysak  and  Szekeres  [48],  who  have  studied 
it  intrinsically  by  the  introduction  of  geodesic  co-ordinates. 

The  Weyl  field 

A 

(3.21) 

v 

where 

,,  00\  2  2  2  ,2  2  /  ,*2 

(3.22)  P  =  r  +  z  ,  p’  =  r  +  (z  -  b)  , 

has  as  its  Newtonian  analogue  the  field  of  a  uniform  rod  of  mass  m  extending 
along  the  axis  of  symmetry  from  z  =  0  to  z  =  b  [49]  .  One  readily  verifies 
that  the  elementary  flatness  condition,  v  =  0  when  r  =  0  ,  is  satisfied 
for  z  outside  the  range  (0,b). 

If  we  set  m  =  ib  ,  (3.21)  is  equivalent  to  the  well-known 
Schwarzschild  spherically  symmetric  solution.  The  co-ordinate  transformation 
(r,z)  ->  (R ,9)  ,  defined  by 

b  i 

(3.23)  r  =  R(l  -  — )2  sin  0  ,  z  -  ^b  =  (R-  A-b)  cos  9 
reduces  (3.I),  (3.2),  (3.21),  (with  m  =  ^b)  to 


P  +  P 


-  b 


m 

b  108  p  +  p'  +  b 


2  -p  log 
b 


(p+p1)2  -  b2 

4pp  * 
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2 

(3.24a)  0  =  +  R2(d02  +  sin2  G  dcp2 )  -  (l  -  £)dt2 

1--  R 
R 

(see  Appendix  II) 

A  discussion  of  this  line-element  is  included  in  every  standard  text,  and 
we  need  not  dally  over  it  here.  However,  there  is  one  remark  which  should 
be  made,  because  it  will  be  important  in  our  later  discussions.  It  has 
been  shown  by  many  authors  [  50,  5I,  52,  53  ]?  that  the  singularity  of  the 
form  (3.24)  at  R  =  b  is  not  a  geometrical  singularity  of  the  manifold, 
but  it  is  due  to  the  nature  of  the  co-ordinates  used.  Kruskal  [54]  has 
exhibited  a  co-ordinate  system  which  covers  the  entire  Schwarzschild 
manifold  without  singularity,  apart  from  the  geometric  singularity  at 
R  =  0  .  Consider  the  line-element 

(3.2ltb)  4  =  +  R2(d92  +  sin2@  dcp2) 

(l+p)eP 

2  2 

where  R  and  p  are  functions  of  £  -  t  defined  by 

(3.25)  R  =  b(l+p)  ,  4b2peP  =  -  t2  . 

? 

2  2  4b 

This  line-element  is  regular  in  the  domain  £  -  t  >  -  — —  .  If  we 
make  the  co-ordinate  transformation  (£,t)  -»(R,t)  where  t  is  defined 
by 

(3.26)  tanh  (™)  =  ~ 

on  the  subdomain  £  -  r  >  0,  we  obtain  the  Schwarzschild  form  (3.24a), 

Kruskal' s  manifold  (3.24b)  is  therefore  an  analytic  extension  of  the 
Schwarzschild  manifold.  It  can  be  shown  that  the  Kruskal  manifold  is 
maximal,  in  the  sense  that  every  geodesic  either  runs  into  the  proper 
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singularity  R  =  0  or  else  can  be  continued  to  infinite  values  of  its 
geodesic  arc  parameter.  It  is  further  immediately  clear  that  (3.24b) 
represents  a  solution  of  Einstein's  vacuum  equations  R  =0  throughout 


2  2  4b 1 

its  domain  of  validity  £  -  t  >  -  - 


ij 


since 


(i)  we  certainly  have  R. ,  =  0  for  the  subdomain  of  the 

ij 

2  2 

Schwarzschild  metric  t,  -  t  >  0  ; 


(ii)  R. .  is  an  analytic  function  of  the  Kruskal  co-ordinates  in 
1 J  p 

.  *2  2  4b 

the  region  £  -  t  >  -  — —  . 


A  detailed  and  very  interesting  discussion  of  Kruskal 's  manifold 
has  recently  been  given  by  Fuller  and  Wheeler  [55]. 

Oblate  spheroidal  shells. 

Chou  [25]  and  Misra  [27]  have  obtained  the  Weyl  field  whose 
Newtonian  analogue  is  a  spheroidal  homeoid,  i.e.,  a  thin  spheroidal  shell 
xtfhose  surface  distribution  of  mass  is  such  that  the  potential  over  its 
surface  is  constant. 

Let  us  consider  the  Newtonian  situation  for  a  moment.  We  recall 
that  oblate  spheroidal  co-ordinates  |  ,  r)  ,  cp  f  are  defined  by 

(3.27)  r  =  b/(l+|2)  (1-T)2)  ,  z  =  b£r)  ,  9  =  9 

in  terms  of  cylindrical  co-ordinates,  r,  z,  (p  «  Spheroidal  co-ordinates 
form  an  orthogonal  net,  as  one  immediately  sees  on  writing  (3.27)  in 
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complex  form; 


(3.28)  r  4-  iz  =  b  cosh[sinh  +  i  sin  ^  r)] 

This  shows  that  the  transformation  (r,z)  ->  (sinh*  1  ?3,sin-1  q)  is  conformal. 
The  surfaces  t]  =  constant  are  the  hyperboloids  of  one  sheet 


(3.29) 


2 

r 


b2(l-r]2) 


1 


The  surfaces  =  constant  are  the  oblate  spheroids 


(3,30) 


b2(l442) 


1 


In  particular,  the  surface  |  =  0  is  the  interior  of  a  disk  of  radius  b 
in  the  xy  plane*  r  =  0  ,  z  <  b  ,  and  the  surface  t)  =  0  is  the  region 
of  the  xy  plane  exterior  of  this  same  disks  r  =  0  ,  z  >  b  .  The 
surfaces  r)  =  +1  and  r)  =  -1  are  respectively  the  positive  and  negative 
segments  of  the  z-axis. 


to  [  56] 


In  terms  of  oblate  spheroidal  co-ordinates,  the  Laplacian  reduces 


(3.31) 


2 

S/  V  = 


,2  /g  2  2  \ 

b  (9  +rl  ) 


8  &2s  8V  8  2 v  8v 

Sf  (1+?  >  5f  +  ^(1-i  }  5^ 


for  an  axially  symmetric  function  V(^,rj).  Separable  solutions  of  Laplace's 
2 

equation  V  =  0  may  be  shown  to  be  of  the  form 


' 
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(3-32)  V(t,n)  =  Unrn(i!)  +  bnQn(ifJ][anPn(,,)  +  PnOn(t))]  . 

th 

Here,  is  the  n  Legendre  polynomial  and  Q  the  Legendre  function 

of  the  second  kind,  defined  by 


(3.33) 


Q>) 


dp. 


(p2-l)P 


2 

n 


Particular  values  are 


/  PQ(p)  =  1,  P1(p)  =  p  ,  ?2(v)  =  i(3n2-l),  - 

Qq(^)  =  ilog  =  coth'V  =  -cot"1^) 

I  Qx(n)  =  iulog  -  l 

^Q2(p)  -  i(3,2-l)  log  H±i  .  |n 


To  obtain  a  vacuum  potential  which  is  constant  over  one  of  the  spheroids 
(3.30),  we  have  to  choose  a  function  (3. 32)  which  depends  on  £  only, 
and  tends  to  zero  at  spatial  infinity  (^  -><»),  We  thus  arrive  at  the 
Newtonian  potential  of  an  oblate  spheroidal  homeoid; 

(3.35)  v  =  — Q0 ( i| )  PD(r,)  .  -  2  cot'hl)  . 


At  large  distances  we  have  from  (3.27) 

(3.36)  p2  =  r2  +  z2  %  b2^2  ,  cot  ^  ,  (|  -»°°) 


00,. 


- 


:  ■  '  ■'  I  5  :■  ■■  :  sJVSi!  ?  #  (,  , 


:■  »■  fti  .".I,-:  s  ii  ,  p  .  .  :  .  v, 
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and  hence 

V  ^  -  —  (p  — » oo) 

showing  that  the  constant  m  in  (3.35)  is  the  total  mass  of  the  spheroidal 
shell . 


To  obtain  the  corresponding  Weyl  field,  it  is  convenient  to  re¬ 
write  (3.2)  in  terms  of  spheroidal  co-ordinates.  This  is  done  in  Appendix 
III  ,  and  similar  calculations  have  been  carried  out  by  Misra  [27]. 

The  result  is 


(3.37a) 


(jV_ 

a? 


(l^2)(l-n2) 

?2  +  ,2  w 


g(l-ig)2  (&f„  2i(l4g)(l-n2)  »  ax 

«V  W  4  v 


(3.37b) 


av 

Br) 


1(1 42)2 


n(i-ii2) (i-t42)  /aA2 
/ft2  2v  lay 

(s  +ri  )  X  17 


2(1- r]2)i(i-t42)  aA  8a 

"  d2+n2)  ^  ^ 


With  the  aid  of  these  formulae,  we  easily  obtain  the  Weyl  field  whose 
Newtonian  analogue  is  the  field  of  an  oblate  spheroidal  homeoid; 


(3.38)  A  =  -  £  cot"1!  ,  v  =  log 

b  2b2  1  4|2 

Let  us  now  consider  the  question  of  joining  this  solution  to  a 
Weyl  manifold  representing  the  gravitational  field  in  the  hollow  interior 
of  the  spheroidal  shell  |  |  [See  ( 3 « 30 ) ] .  If  we  impose  the  reasonable 

requirement  that  A  is  to  be  continuous  across  the  shell,  then  the  interior 


solution  is 


‘ 


•"  :  ’  i  :  -  ' 
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(3.39) 


-v  m  -1* 

A  =  -  —  cot  t  ,  v  =  0  , 
b  o  7 


since  the  only  harmonic  function  A  which  is  constant  over  a  closed  surface 
is  constant  everywhere  inside.  We  observe  from  (5*38)  and  (3*39)  that: 


(i)  The  elementary  flatness  condition,  v  =  0  when  r  =  0  , 

2 

is  satisfied  everywhere  on  the  axis  rj  =  1  . 


(ii)  v  is  discontinuous  across  the  shell  -  i.e.  Weyl’s  canonical 
co-ordinates  are  "inadmissible"  in  the  sense  that  the  metric  tensor  experiences 
jump  discontinuities  across  a  surface  distribution  of  mass. 

(iii)  If  we  replace  ^  by  ^  ,  b  by  be  and  let  e  -» 0  ,  we 
obtain  from  (3.3O)  and  (3*27) 

(3.40)  P  =  bt  ,  r  =  pj  1-T]2  ,  z  =  prj 

so  that  we  obtain  in  the  limit  spherical  polar  co-ordinates  p,  cos  , 

cp  in  the  Euclidean  map.  In  this  limit,  (3*38)  and  (3.39)  become 


(3.41) 


^  ,  v  =  -  (1-T]2)  (Exterior) 

P  2p^ 

=  -  —  ,  v  =  0  .  (interior) 


These  results  should  be  compared  with  our  solution  (3.20)  for  the  field  of 
a  "Curzon  particle".  The  field  (3.41)  might  be  called  a  "Curzon  Shell". 

(iv)  Space-time  is  flat  inside  a  homeoid.  This  is  in  harmony 
with  Newtonian  theory,  according  to  which  the  interior  gravitational  field 


vanishes . 


: 

~  ;  -  j  pio  :r  ■  '•:/  -  nv  r  iv i. :.  i 8i 


■  .-J5  s.  3  : 


•  O  ■  ..  :  !'v  .  :  :>  '  I  -li  /•„.  =  ■  ■  .  , ,  d i 


'.  V  <  .  j  b  •  T.  ■  ;  id  '■  ‘  ..i  '  .  -  ••  1} 
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' 
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Prolate  spheroidal  shells. 


The 


replacements  b  ib  ,  £  -*■  f*  convert  (3.27)  and  (3.29) 


to 


(3.42) 


r  =  b/(l-$2)(n2-l)  ,  z  =  b^rj  , 


(3.43) 


b2(n2-l) 


,2  2 

b  r) 


=  1 


The  surfaces  r)  =  constant  are  prolate  spheroids.  Solutions  of  Laplace's 
equation  have  the  form 


(3.44) 


V  =  [a  P  (£)+b  Q  (£)][a  P  (ri)+(3  Q  (q)]  . 

n  n  ^  n  n  ^  n  nv  1  n  n'  1 


The  Weyl  field 


A  = 


?  <U<i)  -  ^  log 


m 


mi 


(3.45) 


V  = 


v_ 


b  ov  1  2b  &  q+1 


2  2  . 
m  n  -1 

2b2  °\z.e 


corresponds  to  a  Newtonian  field  whose  potential  is  constant  on  the  surfaces 
T]  =  constant.  It  may  therefore  be  labelled  the  field  of  a  prolate  spheroidal 
homeoid.  Noting  that 


(3.46) 


P1  +  P2  *  _  P2  ~  P1 
2b  9  *=  2b 


where 


,  ){  '  ,,-Ci 


■  Z:  itOB  ;  .O  _t  ;fX 

■ 

- 


.  .  . . . 


'  '  cf 


>; 
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/ _  i  _ \  2  2  /  ,  \ 2  2  2  /  ,\2 

(3.47)  Pj  =  r  +  (z-b)  ,  P2  =  r  +  (z+b)  , 

and  comparing  with  (3. 21),  we  recognize  (3.45)  as  representing  also  the 
field  of  a  uniform  rod  of  mass  m  and  extending  from  z  =  -b  to  z  =  b 
along  the  axis  of  symmetry. 

Two  special  cases  are  of  interest. 

(i)  Choosing  m  =^b  we  obtain  the  Schwarzschild  solution, 
representing  the  exterior  field  of  a  spherical  shell. 

(ii)  Replace  b  by  eb  and  r)  by  ,  and  let  e  —>  0  .  We  then 
obtain,  as  before,  the  external  field  of  a  Curzon  shell. 

Solid  spheroids. 

The  Weyl  manifolds  whose  Newtonian  analogues  are  the  exterior 
fields  of  solid  homogeneous  spheroids  have  recently  been  given  by  Misra 
[27],  and  were  independently  computed  by  the  author.  Details  of  the 
calculations  are  given  in  Appendix  IV  .  Here  we  quote  only  the  results. 

For  an  oblate  homogeneous  spheroid 

(3.48a)  A  =  A[4iQo(i?3)Po(q)-4iQ2(i|)P2(q)] 

(A  is  constant) 

=  A[4  cot"1!  +  (3q2~l) {(l+3!2)cot‘“1!-3|} ] 


leading  to 


g- 


: 


■ 


(  ■  "  n-r.  'W  a  ,t  a) 
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(3.48b)  v  =  A2^  (l+£2)  (l+9£2)  (cot  1|)2-(l4^  +l8£^)cot”*l*  +9£2+4] 

-  2A2t]2[ (l+£2) (l+5£2) (cot  1$)2-(6$+10$5)cot"1£  +5^2] 

+  A2[  (l-i^2)  (cot  )2+2(f- -£^)cot  1f'+f-2-4] 

For  a  prolate  homogeneous  spheroid 
(3A9a)  *  =  A[UPo(S)Qo(t1)  -  4P2(£)Q2(r,)  ] 

=  A[4  coth  -  (5|2-1)  ((3T]2-l)coth  ''t]  -  Jr]] 

leading  to 

(3.49b)  v  =  A2^[  (rj2-l )  ( 9rj2 - 1 )  (coth_1r))2  +  (l4T1-l8ri3)coth"1T]+9ri2-4] 

+  2A2|2[  (rj2-l)  (1-5t]2)  (coth  ^ t] ) 2  +(lOT)^-6r})coth  %“5r]2] 

+  A2[  (t]2-1  )2(coth  1t1)2  -  2(r)^+ri)coth  ^t)  +  r)2  +  4] 

Field  of  a  ring: 

One  of  the  first  Weyl  fields  given  explicitly  was  the  field  of 
a  ring,  worked  out  by  Bach  [8]  in  1922.  For  details,  we  refer  to  Bach's 
paper . 

Field  of  a  mass  dipole; 

The  following  result  is  new. 


For  X  take  the  harmonic  function 


•  —  ••  •*»•••  •  -  ■-** •-  . .  . ~ ... . 
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(3.50a)  A  =  22| 

P 

leading  to 

2  2  2  b 

(5.50b)  v  =  -  2m--^  +  ^  m  q-  (See  Appendix  V  ) 

P  P 

Because  of  the  peculiar  character  of  the  "force"  between  a  free  positive 
and  a  free  negative  mass  in  general  relativity  (see  Chapter  V  for  a  discussion) 
the  internal  mechanical  stresses  holding  this  dipole  in  equilibrium  must  be 
of  a  quite  remarkable  nature:  a  pressure  at  the  negative  pole  and  a  tension 
at  the  positive  pole. 

Erez  and  Rosen  [26]  have  given  the  field  of  a  mass  quadrupole 

(5.51a)  A  =  ~{[1  +  iq(5q2-l)(3|2-l)l  log  +  |qri(3^2-l)} 

(3.51b)  v  =  q2^{(9q^~10q2+l)(log^“)2  +  ( 36^-281-])  log^j—  +  36r)2-l 6} 

+  ^2(^2q2(“5qi++6ri2“l)(lo8^^)2  +  (~20rb>+^'rl^loS^J 

+  3q  +  %  q2(“20q2+|^)}  +  (|-  q2+q+|)  log 

+  i5q2(Ti^“2n2+1)(log^)2  +  {j^i2(9Ti2“i5q)-|qq)iog 

+  f^q2(ii2+i)+3q 


where  q  is  an  arbitrary  constant. 


•'  ^  lb  I 


- 


■ 


■  ."i:  :  bck,  ■  >(tf .  J  ■ 
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§  3*4  Field  of  n  bodies:  general  formulae. 


In  the  case  of  a  Weyl  manifold  whose  Newtonian  analogue  is  the 
static  field  of  n  distinct  bodies  B^,  B0,  . . . ,  B^  ,  we  may  set 


(3.52) 


A  = 


A 

i 


where  V  is  the  external  Newtonian  potential  of  B^  ,  and  thus  satisfies 


(3.53) 

V2X. 

1 

=  0  . 

To  obtain  v 

from  (3.2) 

,  we  compute 

v.  . 
ij 

from 

dr 

/dA. 

=  2r(^ 

dA. 

__JL 

dr 

dA. 

1 

dz 

(3.54) 

dv.  . 

0  z 

/dA. 

=  2r(^r 

dA. 

_! 

dz 

dA . 
1 

dz 

and  then 

n  n 

(3.55) 

V  = 

11  ^ 
1— »L 

11  r- 

1— 

v .  . 

rj 

The 

integrability  condition 

of  (3-54)  is 

dA. 

1 

dz 

v  xj  +  sr^ 

A 

=  0 

and  is  satisfied  by  virtue  of  (3.53). 


. 
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§  3*5  Static  n-body  solutions  formal  results. 

We  begin  by  deriving  a  formal  n-body  solution  of  the  equations 
(3*53)  and  (3*54),  deferring  until  later  the  question  of  physical  inter¬ 
pretation. 


Let 


(3.56) 


where 


(3.57) 


m. 

l 


log 


pi  +  pi 


p.  +  p!  +  b. 
l  i  l 


2  2  2  ,2  2  .2 

pi  =  r  +  Zi  ,  p3!  =  r  +  z. 


(3.58) 


zj[  =  z  -  (a1+|bi)  ,  z|  =  z  -  (ai-ib.) 


Formally,  A  =  EA  represents  the  Newtonian  potential  of  a  set 

ttl 

of  collinear  rods  on  the  axis  of  symmetry,  the  i  rod  extending  from 

z  =  t0  z  ai+'2bi  ,  and  having  mass  ,  It  will  be  convenient 

to  define 


(3-59) 


E(i, j) 


2 

p  .p  .  +  z.  z  .  +  r 

1  3  1  3 


E(i’,jO  =  E(j,i') 


P.P 


j 


«  2 
+  z .  z  .  +  r 

1  J 


E(i!,j')=  p!pj  + 


1  t  2 

z .  z  .  +  r 

1  J 


The  integration  of  (3.54)  is  easily  performed  with  the  aid  of 


identities  of  the  type 


(3.60) 


dr  + 


dz 


-1 


d  log[r"  E(i,j')]  = 


2 

r  -  z  z 
i  1 


Zi+Z1 


(for  details  see  Appendix  Vi),  we  obtain 


(3.61) 


ij 


m ,  ra , 
-L.J 

b.b  . 
i  J 


log 


E  ( 

[i'fj 

L)E(l,j') 

E| 

.i,  j) 

We  have  thus  obtained  a  Weyl  manifold  representing  n  bodies,  given  by 

(3.32),  (3»33)»  (3»36)  and  (3.6I).  This  result  is  given  here  for  the  first 
time . 


If  in  (3.56)  and  (3. 61),  we  let  each  b±  f  we  obtain  in  the 
limit  a  solution  which  is  formally  the  field  of  n  collinear  Curzon  particles? 


(3.62) 


m. 

1 


v  = 


v,  .= 

ij 


n  n 


m.m.  z.z,  +  r  -  p  .  p  . 
i L LJ .  !.„! 


p.p . 
1  J 


(z.-z.)4 

v  1  J 


(i  ^  J) 


(details  are  given  in  Appendix  VII).  A  solution  equivalent  to  this  was 
given  by  Hoffmann  [21].  For  n=2  this  solution  was  given  by  Curzon  [12] 


(  !->•.?) 


■  :  ...  ..  .■■-Vi  y/  f  ,/■; 

"" . . 


1 


t  -A.-  A  A 

----- 


"A  •  5  ' 
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in  1924,  by  Silberstein  [15]  in  1936,  and  by  Chazy  [10]  in  I923. 

To  Chazy  is  due  also  a  form  of  the  solution  representing  two 
Schwarzschild  particles,  and  equivalent  forms  were  given  by  Palatini  [11] 
and  (earliest  of  all)  by  Bach  [8],  Bach  gave  his  solution  in  essentially 
the  following  form: 


where 


(3.64) 


«g  -  a!  +  i(b2  -  V  • 


Chazy' s  2-body  solution  was  expressed  in  terms  of  spheroidal 


co-ordinates.  Define  x^,  y^  and  x^,  y^  by 

(3.65)  r  +  i(z-a^)  =  sinh(x1+iy1)  ,  r  +  i(z-a2)  =  sinh (x2+iy2 ) 


The  variables  x^,  y  are  related  to  our  prolate  spheroidal 
co-ordinates  of  (5.42)  by 


r)  =  cosh  ,  £  =  sin  y^  (assuming  =  0) 


Chazy  takes  [cf.  (3.4-5)] 


' j 

C  •  .  ■  . .  - .  t 


±  1  ■  •*  -  ,  :  •:  .  •-  ( 'M:  -v 
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( 3 • 66a ) 


X1  X2 
A  =  lc^  log  tanh  —  +  log  tanh  — 


and  finds 


(3.66b) 


"I  T 

i=l 


sinh^  x. 


log 


2  2 

sinh  x^+cos 


+  ktk2  log 


1  + 


2  sin  yx  -  sin  y2‘ 


a^-a2  cosh  x^+cosh  x2 


All  the  results  of  this  section,  both  for  Schwarzschild  and  for  Curzon-type 
bodies,  are  of  course,  subsumed  in  our  general  formulas  (3.56),  (3.61).  In 
what  follows  we  shall  therefore  concentrate  our  attention  on  these. 

§  3*6  Field  of  n  bodies;  the  elementary  flatness  condition. 

With  the  solitary  exception  of  Bach,  the  authors  of  the  formal 
two-body  solutions,  we  have  reviewed  in  the  previous  section  consistently 
misinterpreted  their  results  as  the  field  of  two  particles  at  rest  on  the 
axis  of  symmetry.  This  was  because  they  ignored  the  vital  condition  of 
elementary  flatness  on  the  axis  of  symmetry.  We  shall  now  examine  the 
consequences  of  this  condition  for  the  general  formal  n-body  solution  (3*56), 
(3.61). 

We  shall  proceed  on  the  assumption  that  the  field  (3.5 6),  (3*61) 
is  produced  by  an  axially  symmetric  system  of  n-bodies,  each  of  which  is 
separated  from  the  others  by  an  empty  region  completely  (no  struts). 

Let  us  also  assume,  for  the  moment,  that  the  n  rods  in  the 
Euclidean  map  be  arranged  in  non-overlapping  sequence  along  the  z-axis: 
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(3.67) 


b£  >  0  ,  +  Jbi  <  a 


i+1  '  -,bi+l 


(i  =  1,2,3, ...,n-l) 


th 


(This  may  not  be  possible  if  the  (i+l)  body  is  hollow  and  encloses  the 
itb  body). 


Elementary  flatness  requires  that  v  =  0  on  sections  of  the  z-axis 
which  are  free  of  material  singularities.  According  to  (3*5*0,  we  have 

=  0  and  hence  =  constant  on  segments  of  the  z-axis  not  occupied 

th  th 

by  the  i  or  j  rods  in  the  Euclidean  map.  Explicitly, 


v.  .  =  0  (r=0,  i  <  i ,  z  >  a  ,+i-b  .  or  z  <  a .  -±b  .  ) 
ij  *  ~  J  2  J  i  2  i 


(3.68) 


(o) 

V  =  V 

Ij  ij 


=  2 


m.m . 
i  .1 

bibj 


log 


(zi~z\ )(z[~z j) 


z.-zTTT z!-z’ 

1  J  1  j 


(r  =  0,  i  <  j  ,  ai+^bi  <  z  <  a  -■§bj) 


(details  are  given  in  appendix  VIII). 


Thus  the  elementary  flatness  imposes  (n-l)  discrete  conditions.  These 
can  be  written  down  by  expressing  the  fact  that  the  jump  in  v  when 
crossing  the  itb  rod  is  zero: 


i-1 


(3.69) 


» 

ij 


(o) 
v.  . ’ 

ij 


=  0 


j=l 


j=i+l 


In  the  Euclidean  map,  v[°.^  represents  precisely  the  Newtonian 
th  th 

force  between  the  i  and  j  rods  [Appendix  IX,  the  result  is  also  obvious 


"  :r"  t(  -  '  '  ■  ■-  -  -■  .«&/ h. .  ■  : 


r  1  ^  .‘i  .f-iAfp., 


•j;,. 


:v  :  ;  y  ,  , 

;v 
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from  our  general  theorem  of  Chapter  II,  page  52],  Hence  (3.69)  may  be 

interpreted  as  the  relativistic  analogue  of  the  condition  that  the  resultant 

tbi 

gravitational  force  on  the  i  body  is  zero. 


For  two  bodies,  (3*69)  reduces  to  the  single  condition 


log 


blb2 


(a1-a2)2.i(b1-b2)2 


=  0  , 


which  cannot  be  satisfied  for  any  non-zero  values  of  the  masses  m^m^ 


For  three  bodies,  (3*69)  becomes 


(o)  (o) 

Vs  '  =  v '  7 

12  23 


,(°) 

13 


or 


(5.70) 


mlrfl2 

blb2 


log 


blb2 


(al~a2^  '^bl"b2^' 


m2m5 

b2b3 


log 


b2b3 


(a2”a5)  “Kb2“b3^ 


mm 

1  3 

blb3 


log 


1  + 


blb3 


L  (ai-aj)2-i(b1+b3)2 


'•■or-: 


(OT.C). 
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If  we  insist  that  all  masses  are  to  be  positive,  then  these 
conditions  cannot  be  satisfied.  (This  is  in  line  with  our  general 
theorem  of  Chapter  II).  But  they  can  be  satisfied,  and  in  infinitely 
many  ways,  if  we  permit  negative  masses.  The  restriction  now  is  that 
m^  and  m  have  the  same  sign,  opposite  to  m^. 


We  shall  now  prove,  using  a  method  whose  basic  idea  is  due  to 
Hoffmann  [21]  that  (3.69)  has  no  solution,  for  non-vanishing  values  of 
the  masses,  when  n  =  4  . 


Define 


(3.71) 


a.  . 

ij 


j 


1 

m,m . 
1  J 


(j  <  i) 

(j  =  D 
(j  >  i) 


then  the  elementary  flatness  condition  (3=69)  can  be  written  as 


(3.72) 


n 


j  =  l 


A  necessary  condition  that  (3.72)  be  satisfied  for  non-vanishing 
values  of  the  masses  is 

det (a, . )  =  0  . 

ij 


(3-73) 


'  '•  '  "■  .  :J  ■■  ||;jt  !'  ;.V-  T 


■:J''  ■■■  *  h$J.  i  3:.:  -A'-  k  A  : X  0  13  jjv  fa&z  0 
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This  condition  is  trivially  satisfied  if  n  is  odd,  for  it  is  well  known 
that  the  determinant  of  an  antisymmetric  matrix  vanishes  if  it  is  of  odd 
order.  If  the  determinant  is  of  even  order,  it  can  be  written  as  a  perfect 
square : 


0,  n  odd 


(3. 7*0 


det (a.  . )  = 

ij 


/ 


L 


,  i.  .  „ 

1  1  n 

-  e  a.  .  a,  . 

L2mml  1112  13i4 


l2 


,a. 

1  , 1 

n-1  n-1 


When  n  =  4,  the  necessary  condition  for  the  elementary  flatness,  (3-73) > 
therefore  becomes 

(3.75)  a12  -  a13  =  0  . 


From  (3.68)  and  (3*7/ )  we  have 


i(al3-a12)  = 


log 


1  + 


bib^ 


(ara3 


)24(b1+b  f 


+  trir*  log 
blb2 


b  b 

1  2 


(a1-a2)2-^(b1-b2)£  J 


< 


(a1-a5)2-i(b] +b^)2  (ai‘a2)2"^(b1‘b2)2 


(ara2)2^(ara3)2^(bl+b3)2^(brb2)2 

[ (a1“a5)2“i(b1+b3)2][(a1-a2)2-i(b1-b2)2] 

b  -b2 

”(a3“a2)  (a2+a5“2a1)+|(b2H-b3)  (b1+ 

[(al"a3)2"Kbi+b3)2H(ai-a2)2^(bi-b2)2] 


i=2m) 


- 


■  \  :  >  '  B 


where  we  have  used  the  Inequalities 


log(l+x)  <  x  (x  >  -1) 

ai"aj  >  i(bi+bj)  (1  >  J) 


the  second  of  which  follows  from  (3.6V) .  Similarly,  we  can  show  that 

*24  -  aj4  <  0  • 


It  follows  that 


a12  a34  -  a13  a24  +  a14  a23  >  0 

and  hence  (3.75)  cannot  be  satisfied.  We  have  shown  that  no  static  solution 

is  possible  for  n  =  4. 

The  results  we  have  obtained  lend  support  to  a  conjecture  we 
shall  call  Hoffmann *s  conjecture; 

Consider  an  axially  symmetric  configuration  in  which  n  bodies 
with  masses  of  both  signs  are  arranged  sequentially  along  the  z-axis .  If 
n  is  even  the  configuration  cannot  be  in  equilibrium^  if  n  is  odd,  there 
are  an  infinite  number  of  equilibrium  configurations. 


We  now  drop  the  sequential  condition  (3.67).  This  admits 
configurations  of  the  type  considered  by  Marder  [20],  The  simplest  situation 
of  this  type  enables  us  to  reinterpret  the  two-body  solution  given  in  §  3»6. 


3:  <•,- 


>■ 


Figure  4, 

(Body  inside  Body  B^ ) 

In  the  Euclidean  map  let  the  first  body  B^  be  thin  hollow 
spheroidal  shell  with  mass  and  foci  at  +  -|b^.  We  assume  the 

equation  of  the  shell  to  be  +  p|  =  c^  where  >  b^  .  Let  the  rod 
B^  extending  from  (a^  -  -Ib^)  to  (a2  +  ^2^  al°n8  the  axis  representing 
a  second  body  enclosed  within  the  shell. 

For  the  field  exterior  to  the  shell,  we  have  the  two-body 


solution  of  §  3«6: 
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C,  .  Pl+Pl~bl 

-  bi  °8  P1+P|+b1 


m  °2+02'b2 

+  r*  log  -f  ■  ^ 

b2  p2+P2+b2 


(3.76)  < 


2  /.  ,v2  u2  2 


,2  .  2 


v  = 


m-  (P1+P1)  ’bl  .  m2  .  (p2+ppC-b2 


2  bf  108  bplpi 


+  2  —  log 

bi 


4p2°' 


„  mlm2  ,  E( l1 ,2)E( 1.21 ) 
+  2  b1b2  log  E(1(2)e(1\2') 


(pl+pi  >  Cj) 


The  field  in  the  interior  of  the  shell  is  given  by 


r 


Vbi 

c  i+b  i 


+ 


P2+P2'b 

P2+P2+b 


2 

2 


(3.77)  H 

V 

L 


log 


(p2+Pp)-b2 

^2P2 


( p i+p i  <  ci^ 


We  turn  to  the  elementary  flatness  condition,  v  =  0  when 

r  =  0.  It  is  easy  to  check  that  this  is  satisfied  on  segments  of  the 

z-axis  exterior  to  the  shell,  for  which  p.  =  z.  ,  p!  =  z!  ,  r  -  0 

1111 

(upper  segment)  or  p^  =  -z  ,  p^  =  -z|  ,  r  =  0  (lower  segment). 

For  segments  of  the  z-axis  interior  to  shell,  (3.77)  shows  that  y  =  0 
except  on  the  section  actually  occupies  by  the  rod?  P2  +  P2  =  b°  * 

In  (3»76) »  (3*77)  we  thus  have  a  two-body  solution  satisfying 
elementary  flatness  and  in  which  the  masses  of  the  two  bodies  are  quite 


arbitrary. 


■ 

:  :  :■  .... 


.  .  ■  ' 


.  »;  .  .  ,i 
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The  solution  given  above  differs  from  that  originally  given  by 
Marder.  Marder  considered  a  thick  hollow  shell,  and  used  a  Curzon-type 
solution,  which  is  more  specialized  than  the  one  given  above.  However, 
both  solutions  bring  out  the  essential  features  equally  well. 
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CHAPTER  IV 

STATIONARY  AXIALLY  SYMMETRIC  VACUUM  FIELDS 

§  4.1  Introduction. 

In  Chapter  II,  we  defined  what  is  generally  meant  by  a  stationary, 
axially  symmetric  space-time  manifold  (pagel'7).  We  showed  that  co-ordinates 

123  4 

x  ,  x  ,  x  =  (p,x  =  t  can  be  introduced  in  terms  of  which  the  metric 
form  reduces  to 

(4.1)  $  =  e^(dx^  +  dx^)  +  idCp^  +  2mdcpdt  -  f  dt^  , 

where  all  metric  coefficients  are  functions  of  x^,  only.  The  form 

of  this  metric  is  clearly  unchanged  by  any  conformal  transformation  of 
x^,  x^  .  The  Einstein  field  equations  associated  with  (4.1)  have  been 
studied  by  Andress  [30],  Lewis  [34],  Van  Stoekum  [35],  Papapetrou  [36]  and 
Tiwari  and  Misra  [37] •  To  Lewis  is  due  the  important  discovery  that  the 
vacuum  field  equations  can  be  significantly  simplified  by  the  introduction 
of  so-called  'canonical  co-ordinates"  (x^,  x^).  (These  co-ordinates  are 
similar  to  the  canonical  coordinates  of  Weyl  for  the  static  case  (page  24), 
and  in  fact  reduce  to  them  if  m  =  0.)  Even  with  this  simplification, 
solving  the  vacuum  field  equations  still  entails  finding  solutions  of 
two  coupled  non-linear  partial  differential  equations  in  two  unknown 
functions.  To  date,  the  only  result  obtained  along  these  lines  which  seems 
to  have  even  a  remote  physical  significance  is  the  one -dimensional  solution 


. 
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found  by  Lewis  and  discussed  at  some  length  by  Van  Stockum.  This  presumably 
represents  the  external  field  of  a  rotating  infinite  cylinder.  Van  Stockum 
has  also  given  a  corresponding  interior  solution,  representing  a  rigidly 
rotating  cylinder  of  dust  held  in  dynamical  equilibrium  by  gravitation  and 
centrifugal  force.  Other  solutions  have  been  given  by  Misra,  but  their 
interest  appears  to  be  mainly  mathematical. 

A  few  months  ago,  R.  P.  Kerr  [38]  reported  briefly  on  a  stationary 
axially  symmetric  vacuum  field  which  he  had  obtained  along  quite  different 
lines,  and  which  appears  to  represent  the  external  field  of  a  rotating 
sphere  or  ring. 

The  plan  of  this  chapter  is  as  follows. 

In  §  4,2,  we  derive  the  vacuum  equations  for  the  metric  (4.1), 
and  introduce  the  canonical  co-ordinates  of  Lewis.  In  §  4.3,  we  reduce 
the  problem  of  solving  the  vacuum  equations  to  the  solution  of  two  coupled 
partial  differential  equations.  Some  of  the  developments  of  this  section 
are  new.  In  §  4.4  we  give  a  new  and  very  much  simplified  derivation  of 
Lewis’s  solution  for  the  external  field  of  a  rotating  cylinder.  The  final 
section,  §  4.5,  is  devoted  to  a  very  brief  discussion  of  Kerr’s  solution. 

§  4.2  Field  equations.  Canonical  co-ordinates. 

Consider  the  metric  form 

(4.2)  <b  =  ga[3  dxa  dxP  +  g^  dXA  dxB 


■ 


-  ■'  ■  • ■ ::: 

■  r?;.1  1  ,  >;> 
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where 


(4.3) 


8aB(VX2>  -  gAB  '  SAB(X1’X2) 


(a, 3,7  etc.  =  1  or  2  ;  A,  B,  C  etc.  =  3  or  4)  which  is  a  slight  generalization 
of  (4.1).  We  denote  by  a  stroke  covariant  differentiation  with  respect  to 
the  two-dimensional  metric 

(4.4)  gap  dxa  dxP 


and  a  comma  indicates  partial  differentiation,  e.g. 

SAB,(i  8AB  ° 


Note  that  gA„  is  a  scalar,  and  gAm 

Ais  '  "  n" '  “  Ai> ,  |i. 

in  the  2-space  of  the  metric  (4.4). 


a  vector  with  respect  to  transformations 


The  components  of  the  Ricci  tensor  corresponding  to  the  line- 
element  (4.2)  have  been  given  in  an  elegant  and  concise  form  by  Van  Stockum 
[35]  t  whose  presentation  we  follow  here. 


Write 


(4.5)  C2  =  -  detail 

(4.6)  llg^l  -  HgjT1  ,  ||gaP||  =  ||gap|| 


Then  the  Christoffel  symbols  for  the  line-element  (4.2)  are  as  follows: 


j  .  ;  •>  I  ;  5  ;  i  , 


:  .  ~ '  Vi 


h 
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4  -  0 


(4.7) 


<X 


’  4  -  0  ’  7a  ■  0 


Ba 


,  AC 

=  #8  8 


BC,a,  AB 


jA  l  ap 


fS  g 


AB.  3 


For  the  components  of  the  Ricci  tensor,  we  then  find 


(4.8) 


RP  -  +  D'1(V  *  V'lX** 

^  ^D‘1(D8ACSBC)a)|psaS 


where  is  the  Ricci  tensor  of  2-space  (4.4) 


We  now  specialize  to  the  particular  metric  (4.1),  i.e.  we  take 


(4.9) 


lls 


'AB1 


11  a3 


11*  “  I 

"m  -f' 


D  =  if  +  m 


2lp 

.e  o 
'o  e 


2*11  • 


Since  space-time  is  normal-hyperbolic,  it  is  necessary  that  D  defined 
by  (4.9)  be  positive,  as  we  have  anticipated  by  the  notation. Equations 
(4.8)  now  reduce  to 


(4.10)  g  bJ  =  D(A^)£l  +  D  |Q  -  I^T1 


irD  ~(i  f0-f-i0f_+2m^mQ) 


*P  "v  r/  P  ,ajp  *  '%a‘,3  '  "4 P‘, a  ' 


(4.11)  •f-g  = 


i3a 


f  i  +tn  m 

,a  ,a 

D 
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(4.12) 

(4.13) 

(4.14) 


>T7*3 

4 


rfm  -  mf 
fg ,a 


D 

L 

- 

rmi 

£m  -i 

fg 

D 

rif  + 

,a 

mm  ~i 

,g 

D 

L 

where  A^  = 


f  +ii 

r,aa  >.2  2 

1  ox^  dx2 


We  proceed  to  explore  the  consequences  of  the  vacuum  field 


equations 


(4.15) 


Rm  =  0 

n 


From  (4.11)  and  (4.14)  we  have 


(4.16) 


0  =  -/^i(Rj  +  r£) 


=  4-3. 


g 


(if+m2) 


g 


AD 


Hence  r  =  D(x^,x2)  is  a  harmonic  function  In  vacuo  (more  generally! 

1  2 

medium  whose  energy  tensor  satisfies  =  0)  .  Let  z(x^,x2) 

the  conjugate  harmonic  function,  so  that  the  analytic  transformation 
(x1#x2)  ->  (r,z): 

r  +  iz  =  f(x^  +  ix2) 


in  a 
be 


;  ..  :•’)  . 

* 
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is  conformal.  This  transformation  leaves  the  form  of  the  metric  (4.1), 
and  hence  the  form  of  equations  (4.10)  to  (4. 14), unchanged . 

We  can  thus  rewrite  our  previous  equations  in  terms  of  the 
"canonical  co-ordinates"  (r,z)  .  We  have 


C^-17) 


3>  =  e 


2^r'Z)  (dr2  +  dz2)  +  i(r,z)  dep2 


+  2m(r,z)  depdt  -  f(r,z)  dt  , 


with 

(4. 18)  £$,  +  m2  =  r^ 

The  Ricci  tensor  for  this  metric  is  given  by  (4.10)  to  (4.14),  with  D  =  r. 
Equation  (4.16)  is  now  an  identity. 

The  introduction  of  the  Lewis  canonical  co-ordinates  r,z  and  the 
simplifications  attendant  upon  their  introduction,  are  always  possible  in 
vacuo  or  in  a  medium  whose  energy  tensor  satisfies 

1  2 

Tx  +  T2  =  0  , 

since 

(t\  +  T~)  =  -  K[R*  +  R2  -  i&jR  -  ^52R] 

=  K(R^  +  r£)  ' 

=  —  AD 

g" 


'  i  .  rjo 1  ba<; 

■ 

-  '  r';  ■  .  '  .  of 

-  ■■■.  tr  ,  ;  co  >  : 
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by  (4.16).  We  observe  from  (4.17)  and  (4.18)  that,  when  m  =  0,  r,  z 
reduce  to  the  canonical  co-ordinates  of  Weyl  for  the  static  case  [c.f. 
Chapter  II,  page  24]. 

We  proceed  to  unpack  the  remaining  equation  of  (4.15).  From 

(4.11),  and  (  .14), 


which  yields 


^g(R  \ 


K)  =  # 


a 


U 


it 


0 


(4.19)  fA*£  =  iA*f 

where  (for  any  function  ^(r,z)) 


(4.20) 


1, 

P,i 


h2* 

ar2 


1  ajj/ 

r  dr 


+ 


i± 


The  equations 


(4.21) 

(4.22) 


5  4 

RJ  =  0  ,  R,  =  0  yield 

mA*i  =  iA*m  , 

fA*rn  =  nA*£  . 


Any  one  of  the  three  equations  (4.19),  (4,21),  (4,22)  is  a  consequence  of 
the  other  two. 

2  2 

Operating  with  A*  on  the  relation  if  +  ra  =  r  ,  and  noting 
A*r^  =  0  ,  we  find 
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(4.23) 


iA*f  +  fA*i  +  2Vf  •  VI  +  2inA*m  +  2 (V  m)"  =  0 


where,  for  any  two  functions  <p(r,z),  ^(r,z),  we  define 


(4.24) 


Vcp  -Vf  »  +  Sg> 


5r 


6z 


From  (4.19)  and  (4.23) 


f  A*  i  =  i(f  A*  i  +  i  A*  f) 


m 


A*  m  -  Vf  •  Vi  -  (Vm)‘ 


j-  A*  i  -  Vf  *  Vi  -  (H2  by  (4.21) 


2  2 

Hence,  using  fi  +  m  =  r  , 


(4.25) 


A*  £  +  °  \yi  +  (vm)4  =  0 


From  (4.19)  and  (4.21),  it  then  follows  that 


(4.26) 


r2  A*  f  +  f  [  V  f  *  Vi  +  (Vm)2]  =  0 


(4.27) 


r2  A*  m  +  m[Vf  »  Vi  +  (Vm)^]  =  0 


We  turn  to  the  equations  =  0  (a, 3  =  1  or  2).  We  have, 


since  D  =  r 


D  =  1  ,  D 

,  i  } 'd 


D 


,a|3 


ap 


0 
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The  equations  R  =  0  give 


’,1|2  -  47  (i,l  f,2  +  l,2  f,l  +  2m,l  ra,2> 


=  0 


which  can  be  written 

i 


(4.28) 


1_  df  di 

4r  dr  ’  dz  +  dz 


df  0dm  dnu 

57  +  2S7  •  S7' 


1  2 

The  equations  =  0  lead  to 


D 


1|1  '  U,2|2 


77 


2il  f  +  2m 

>  ^  » *-  ) 1 


or 


(4.29) 


dr  ~ 


i  r^i  af  af  af  /dm\2  ,anu2 

4r  °  ar  az  az  +  ar  'az 


The  integrability  condition  of  (4.28)  and  (4.29)  is 


(4.30) 


M  a*  f  +  jr  +  2^  A*  m  =  0 

dz  dz  dz 


and  is  identically  satisfied  because  of  (4.18),  (4.21)  and  (4,22). 


1  2 

It  is  routine  to  check  that  the  equations  =  0  are 

consequence  of  (4.28)  and  (4.29),  so  that  the  system  (4.23)  to  (4.29) 

2  2  m 

ff  +  m  =  r  ,  is  fully  equivalent  to  the  vacuum  equations  R  =  0 


2™?2 )  *  ° 


a 

,  with 


'  ■  U  ■'  v:*  <J')i 
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§  4.3  Reduction  of  the  vacuum  field  equations. 

We  have  already  pointed  out  that,  of  the  three  equations  (4.25), 
(4.26)  and  (4.27)  for  the  three  functions  i,  f  and  m  ,  one  is  redundant. 
This  can  be  regarded  as  a  consequence  of  the  identity 

(4.31)  if  +  m2  =  r2  . 

We  shall  now  make  this  redundancy  more  apparent  by  reducing  the  three 
functions  i,  f,  m  to  effectively  two. 

First  let  us  write 

(4.32)  i  =  ru  ,  f  =  rv  ,  m  =  rw 
so  that  (4.31)  becomes 

(4.33)  uv  +  w  =  1  .  We  have,  for  any  function  tp( r,z) 

A*(r ,ip)  =  r  A*  ip  +  ip  A*  r  +  2S7  ip  •  7r 

=  r  A-^  -  ^  +  2  | 


=  r  V2  ip 


where 


(4.34) 


V2  ip  (r,z) 


2 

y'  ip  is  formally  the  Lapacian  for  cylindrical  co-ordinates  r,  z,  cp  m 
Euclidean  space.  Equations  (4.25),  (4.26)  and  (4.27)  now  become 


V 
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(4.35) 


(4.36) 


(4.37) 


V2. 


w  +  Pw  =  0 


where 


(4.38) 


rP  =  -1  +  vf  '  V&  +  (<?™)2 


r2[Vu  *Vv  +  (V«)2] 


We  have  made  use  of  (4.33)  in  simplifying  the  expression  for  P  . 

We  now  replace  i,  f,  m  by  two  algebraically  independent  functions 
p(r,z),  0(r,z)  in  such  a  way  that  (4.33)  is  identically  satisfied. 

2  2  2 

Case  I  ;  Assume  m  <  r  ,  i.e.  w  <1  , 

Let 

(4,39)  u  =  e”^sin  0  ,  v  =  e^  sin  0  ,  w  =  cos  0  . 

This  satisfies  (4.33)  for  all  P  and  0  . 

We  then  obtain  from  (4.38) 


(4.40) 


P  =  r[(ve)2  -  Sin2  0  (Vp)2] 


Equations  (4.35),  (4.36)  and  (4.37)  reduce  to 


?  2  2 
sin  6>[v  p  +  2(cot  0)vp  •  y  ^  1  ~  c°s  0[V  0  +  sin  9  cos  9  (v'p)  ]  =  0  , 


sin  9[S^\x  +  2(cot  0>?P  *^0]  +  cos  0  [V2©  +  sin  0  cos  0  f/p)2]  =  0 


“  77  - 


(4.41)  V^0  +  sin  0  cos  0  (y  p)‘  =  0  . 

Thus ,  (4.35),  (4.36)  and  (4.37)  are  equivalent  to  the  two  equations  (4.4l) 
and 

(4.42)  V2 [i+2  cot  e(V|a  •  V@)  =  0 

Equations  (4.29)  and  (4.30)  can  be  written,  in  view  of  (4.32) 
and  (4.40), 


(4.43) 


dr 


1_ 

4r 


.  2  .  ,  2  2  v. 

sm  0  (|a  -  P  o)] 

,1  ,  CL 


(4.44) 


r 

2 


0  •  ^ 

9  L 


We  thus  have  the  following  prescription  for  obtaining  stationary 

axially  symmetric  vacuum  fields  with  metric  (4,17),  in  the  case  where 

2  2 

m  <  r  .  Solve  the  coupled  equations  (4.4l)  and  (4.42)  for  p.(r,z)  and 

0(rtz).  Integrate  (4.43)  and  (4.44)  to  get  r,z),  The  coefficients 

£,  mt  f  in  the  metric  (4.17)  are  then  given  by 

(4.45)  £  =  r  e’^sin  0  ,  f  =  r  e1^  sin  0  ,  m  =  r  cos  0 

Unfortunately,  the  system  (4,4l),  (4.42)  is  in  general  intractable. 

Two  special  cases  are  of  interest. 

(i)  If  we  set  0  =  ^  ,  we  come  back  to  the  static  case, 


m  =  0 


. 

,  -  :  '7 ■  -f 
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Equations  (4.41)  and  (4.42)  reduce  to 


(4.46) 


V  |i  =  0  , 


while  (4.43)  and  (4.44)  become 


(4.47) 


dr 


1  r  /  2  2  \ 
4r  +  4  ^2^ 


(4.48) 


dz 


2  ^,1  4,2  ' 


These  are  precisely  the  Weyl-Levi-C ivita  equations  (2,47)  and 

(2.48)  of  Chapter  II,  if  we  set 


p  2A  ,  , 

re  =  e  ,  ^  =  y  -  A 


(ii)  Set  p  =  log  K  (constant).  Then  (4.4l)  and  (4.42)  reduce  to 

(4.49)  V20  =  0 

and  (4.43),  (4.44)  become 


(4.50) 


d£ 

dr 


1_ 

4r 


] 


(4.51) 


Q 


,2 


The  corresponding  metric  is 

(4.52)  $  =  e^(dr2  +  dz2)  +  r  sin  @(Kdcp2  -  dt2)  +  2r  cos@  d9  dt 

Since  (4.49)  is  formally  Laplace's  equation  in  cylindrical  co-ordinates,  a 
large  class  of  vacuum  manifolds  of  the  type  (4.52)  can  immediately  be  found. 


' 
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This  class  of  solution  was  first  isolated  by  Tiwari  and  Misra  [37 J. 

2  2  2 

Case  II:  We  now  suppose  m  >  r  ,  i.e.  w  >  1,  corresponding  to  a  field 
in  which  rotational  effects  are  very  strong. 

Equation  (4.33)  is  automatically  satisfied  if  we  set 

(4.53)  u  =  e  ^sinh  0,  v  =  -e^  sinh  0,  w  =  cosh  0 

2 

We  observe  that  this  makes  the  coefficient  of  dt  positive,  or  else  the 

2 

co-efficient  of  dcp  negative.  However,  the  normal  hyperbolic  character  of 

the  metric  is  left  inviolate,  since  it  is  entirely  determined  by  the  positive- 

r  2 
ness  of  if  +  m  . 

Equation  (4.37)  yields 

(4.54)  P  =  -r[(y0)2  -  sinh20(via)2]  . 

Equations  (4.35)>  (4.36)  and  (4.37)  become 

p 

sinh  0[^(j+  2  coth  0y  p  •  y0] 

2  2 

-  cosh  0fy  0  +  sinh  0  cosh  0  (vp)  ]  =  0  , 

p 

sinh  0[v  |i  +  2  coth  0yp  •  V03 

2  2 

+  cosh  0[y0  +  sinh  0  cosh  0  (yp)  ]  =  0  , 

p  p 

(4.55)  y  0  +  sinh  e  cosh  0  (vpO  =  °  » 

and  are  therefore  equivalent  to  the  two  equations  (4.55)  an^ 

y2p  +  2  coth  0  Vp  •  V0  =  0  . 


(4.56) 
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Equations  (4.28)  and  (4.29)  can  be  written,  by  virtue  of  (4.32)  and  (4.53)> 


(4.57) 


1  r  ra2  2  .  .2  af  2  2  v, 

57  +  4  [0,1  '  e,2  -  Slnh  e(tx,l  -  “.a*1 


(4.58) 


|(Sirh2e  M  IX  -  e  e>s) 


In  summary,  stationary  axially  symmetric  vacuum  fields  (4.17)  may 
be  obtained  by  solving  the  system  (4.55)  to  (4.58),  and  then  setting 


(4.59) 


"M-  [1 

f  =  re  sinh  G  ,  &  =  -re  sinh  Q  ,  m  =  r  cosh  6 


A  special,  tractable  class  of  solutions  may  again  be  obtained  by 
setting  p  =  log  k  .  ‘The  metric  (4.17)  then  reduces  to 

(4.60)  $  =  e^/(dr^+dz^)  -r  sinh@(^  d<p2  +  Kdt^ )  +  2r  cosh  6  dcp  dt  , 

K 

with  @( r,z)  and  ^/(r,z)  determinable  from 


(4.61) 


71  =  o 


(4.62) 


dr 


h  +  f(^2i  -  e22) 


4r  T  4W,1 


(4.63) 


dz 


G  G 


1  ,2 


The  above  reduction,  and  the  isolation  of  the  special  class  of 
solutions  (4.60),  have  not  previously  been  given  in  the  literature. 


V'  .  ar.?3  r^j. 
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§  4.4  Lewis's  solution.  A  new  derivation. 


In  the  special  classes  of  solutions  (4.60),  (4.52)  which  we 

isolated  in  the  previous  section,  let  us  assume  0  to  be  a  function  of 

2 

r  only.  To  satisfy  V  0  =  0  ,  we  have  to  take 


(4.64) 


o  =  cl  log  r  +  log  cg 


2  2, 

For  the  case  m  >  r  ,  (4.62)  and  (4.63)  then  yield 


2 

ip  =  -^(c^  -  1)  log  r  +  const, 


2  2 

while  for  the  case  m  <  r  ,  we  have  from  (4.50)  and  (4. 51) 


2 

ip  =  ~  A(C;l  +  1)  log  r  +  constant 


We  thus  arrive  at  the  two  stationary,  cylindrical  symmetric  vacuum  metrics 


2  .  2,  w 

(4.65)  *  =  Ar  1  <dr  +dz  }-i(c2r 


1+C1  1  1“C  1  O  O 

-  ”  r  )  (|  dcp2  +  K  dt2) 

c2 


14C  t  1-C 
+  (c  r  +  r 

^  C2 


)  d  9  dt 


-Me  +1) 

(4.66)  $  =  Ar  1  (dr2+dz^)  +  r  sin(c1  log  r  +  cQ)  (K  dcp2  -  ~  dt2  ) 


+  2r  cos(c^  log  r  +  cp)  d  cp  dt 


The  metric  (4.65)  is  the  solution  originally  given  by  Lewis  [34] 
here  obtained  by  a  new  route.  The  metric  (4.66)  was  given  by  Van  Stockum 


■  «r. :  1 b  !;? ■  ixl 

'  "-!  ■■  !  ■  mo  ■  v-ij'iq  ari s  eti  bssjsio&j: 


!  '  •  ■  :■  .  ■  'ifi'.: ■• '  '  »■  l  ~xt-  saw-  ■  -.vv 
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[  3Cl  ]  who  has  interpreted  both  solutions  as  the  external  fields  of  rotating 
infinite  cylinders,  and  has  shown  how  they  can  be  ■joined  to  suitable  interior 
solutions.  We  shall  not  pursue  these  questions  further  here. 

§  4.5  Kerr's  solution. 

The  remarkable  and  important  solution  recently  obtained  by 
R.  P.  Kerr  [38]  was  derived  by  a  method  which  bears  no  resemblance  to  that 
developed  for  stationary  axially  symmetric  fields  in  §§  4,2  and  4.3,  It 
may,  in  fact  be  fairly  described  as  an  accidental  and  very  fortunate  by¬ 
product  of  an  investigation  not  directly  connected  with  the  problem  of 
axially  symmetric  fields.  Details  of  Kerr's  calculations  have  not  yet 
been  published,  and  it  will  be  impossible  to  give  them  here.  But,  before 
presenting  his  explicit  solution,  I  shall  attempt  to  indicate  in  a  qualitative 
way  the  kind  of  problem  he  was  concerned  with. 


An  algebraically  special  vacuum  field  is  a  space-time  free  of 


matter  and  electromagnetic  radiation  (R, .  =  0)  in  which  there  exists  at 

^  J 


,a 


least  one  null  vector  field  K  satisfying 


(4.67) 


R  n  r?-  K  ,KP  K7  =  0 

apyLS  eJ 


Goldberg  and  Sachs  [57l  have  proved  that;  a  vacuum  field  is 
algebraically  special  if  and  only  if  it  contains  a  shear-free  null  geodesic 
congruence . 

A  null  geodesic  congruence  is  a  family  of  null  geodesics  exactly 
one  of  which  passes  through  each  point  of  the  manifold.  The  congruence  is 


"  '  '  jS  '  ,  '  ;•  . 

O  .*■  .  :  •'  ?  •:  ;  •  '  '  ;  «' 

^  -  .  :■  "  '  j>  ? 

i  J  ;  ■  ■  ■  •  )■,  i  ;  '■  ;■  ■;  .  / 

" 
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said  to  be  shear-free  if  the  field  of  null  tangent  vectors  to  the  geodesics 
satisfies 

<Ka|P  +V  ^  • 

The  geometrical  significance  of  this  condition  has  been  explained 
by  Sachs  [38].  According  to  the  above  theorem  of  Goldberg  and  Sachs  the 
null  vector  field  satisfying  (4.68)  in  an  empty  space-time  also  satisfies 

(4.67). 


In  studying  algebraically  special  fields,  it  is  natural  to  focus 

attention  on  the  null  vector  field  K  and  its  variation  along  the  null 

a 

geodesic  rays.  This  is  accomplished  by  choosing  K  as  one  of  the  vectors 
of  an  orthogonal  tetrad.  Let  m^  be  a  second  null  vector,  satisfying 

a  »  „a  . 

mm  =  0,mK  =1. 
a  a 

In  the  space-like  2-flat  orthogonal  to  the  2-flat  containing  m^, 
choose  a  pair  of  orthogonal  vectors  Y^  ,  Z  ,  which  thus  satisfy 

a  ~  ’  a  ”  ’  a 

Y  ma=Z  ma  =  Y  K<X=Z  Ka  =  0  . 

a  a  a  a 

It  is  then  e asy  to  see  that  we  can  express  the  metric  tensor  in 

the  form 


g  _  =  m  K  +  m  K  +  Y  Y  +  Z  Z. 

aP  a  p  Pa  ap  aP 


(4.69) 


1 
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for  we  merely  have  to  check  that 

g  is  defined  by  (4.69).  (Since 

ap 

the  four-dimensional  manifold,  it 
and  three  similar  relations.) 


g  n  =  A  for  all  vectors  A 
aP  a 

the  vectors  m^, 


is  enough  to  verify  that  g  trv~ 

exp 


^  when 

span 

=  m 
a 


It  is  customary  to  introduce  the  complex  null  vectors 


:a  =  —  (Ya  +  iZa) 


ta  ,  -i-  (Ya  -  iZa) 

/2 


so  that 


1 


and  (4.69)  can  be  written 


g  „  =  m  K.  +  m  K  +  t  t  +  t  t 
aP  aP  Pa  ap  ap 


The  method  used  by  Kerr  and  others  in  treating  algebraically  special 
vacuum  fields  is  to  re-write  the  field  equations  R„  .  =  0  and  the  condition 

ij 

(4.67)  in  terms  of  the  tetrad  vectors  K,  m,  t,  t  .  I  shall  not  attempt  to 
go  into  the  complicated  details  of  this  procedure. 

In  the  course  of  his  investigation  of  algebraically  special  fields 
Kerr  found  the  following  particular  vacuum  field: 


(4.70)  $  =  (R^+a^cos^@ )  (d@^+sin^@  dcp^)+2(du+a  sin^Odcp)  (dR+a  sin^@  dep) 


2mR  \  1  j  ,  2„  j^s2 

- — - )  (du  +  a  sin  9  dep) 


~  ( 1  - 


2  2  2 
R  +a  cos  9 


’  '  ■  ;  r 

. 

. 

‘  1  ;  .  ’  :  "!■  X 


I 
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It  is  clear  that  this  line-element  is  axially  symmetric  (metric  co-efficients 


independent  of  cp)  and  "stationary"  in  the  generalized  sense  that  the  metric 


co-efficients  are  independent  of  u  .  Whether  it  is  stationary  in  the 
narrower  sense  of  Chapter  II  (page  17 :  invariance  under  simultaneous 

inversion  of  "time"  and  reversal  of  the  sense  of  rotation)  is  not  clear 
from  the  form  (4.70).  In  fact,  the  entire  relationship  of  Kerr's  solution 
to  the  Lewis -type  approach  remains  mysterious  and  would  be  an  interesting 
topic  for  detailed  investigation. 


If  in  (4.70),  we  set  the  parameter  a  equal  to  zero,  we  obtain 


(4.71)  0  =  R2(de2  +  sin2e  dcp2)  +  2  du  dR  -  (l-  ™)  du2  . 

K 


This  line-element  represents  an  analytic  extension  of  the  Schwarzschild 
manifold  [c.f.  Chapter  V]  which  resembles  Finkelstein' s  well-known  extension 
[52].  To  convert  (4.7*)  into  its  more  familiar  form,  introduce  a  new  co¬ 
ordinate  T  ,  defined  for  R  >  2m  by 


(4.72) 


u 


R  +  T  +  2m  log  (~~  -  l) 


so  that 


0*.73) 


du 


R 


The  line-element  (4.71)  becomes 


n2/ ,_2  „  2  *  j  2 •.  2ms  , -,2 

R  (d@  +  sm  0  dcp  )  -  (1-  tt-)  dT 

K 


ni 
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Thus,  Kerr's  solution  includes  Schwarzschild 1 s  as  a  particular  case  (a=0 ) 

Kerr  has  indicated  another  way  in  which  the  general  solution 
(4.70)  can  be  expressed  in  terms  of  asymptotically  Galilean  co-ordinates 
x,y,z,t  .  Make  the  transformation 


r*  - 


(u.74) 


x  =  (R  cos  cp  +  a  sin  cp)  sin  9 


S  y  =  (R  sin  cp  -  a  cos  cp)  sin  9 


V 


z  =  R  cos  9  ,  u  =  R  +  t 


We  then  obtain  from  (4„70), 


(^.75) 


A  1  2  ,2  ,2  ,2  2imr  2 

©  =  dx  +  dy  +  dz  -  dt  +  — ^ — g— g  ^ 

r  +a  z 


where  q  stands  for  a  Pfaffian  form  defined  by 

2  2  2  2  2 

(4.76)  (r  +  a  )rq  =  r  (xdx  +  ydy)  +  ar(xdy  -  ydx)  +  (r  +  a  )(zdz  +  rdt) 

where  the  function  r  is  defined  by 

/ 1  \  4  ,  2  2v  2  2  2  2  2  2  2 

(4.77)  r  -  (p  -a  )r  -  a  z  =0,  p  =x  +  y  +  z 

For  r  »  a  we  have 


r  =  p  +  O(^)  (§  -*  °°)  . 

p  a 

In  this  co-ordinate  system  the  solution  is  regular  everywhere  except  at 
p  =  a  ,  z  =  0 


'H  -*,f ' 


rx:  . 


•  '  1  •'  ;■  C 


/  .  , 

A  ' 
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Kerr  reports  that  he  has  compared  the  field  (4.75)  at  large 
distances  with  known  approximate  results  for  the  exterior  field  of  a 
rotating  sphere.  There  is  agreement  to  the  third  order  of  the  Einstein- 
Infeld-Hof fmann  approximation  if  m  is  taken  to  be  the  mass  and  ma 
the  angular  momentum  about  the  z-axis. 

Very  recently  Newman  and  Janis  reported  to  the  January  1964 
meeting  of  the  American  Physical  Society  in  New  York  on  a  new  interpretation 
of  the  Kerr  solution.  The  abstract  [ 39 3  of  their  (as  yet  unpublished) 
paper  reads  as  follows? 


"It  is  shown,  that,  by  means  of  a  complex  co-ordinate 
transformation  performed  on  the  monopole  or  Schwarzschild 
metric,  one  obtains  a  new  metric  (first  discovered  by 
Kerr).  It  has  been  suggested  that  this  metric  be  inter¬ 
preted  as  that  arising  from  a  spinning  particle.  We 
wish  to  suggest  a  more  complicated  interpretation; 
namely,  that  the  metric  arises  from  a  ring  of  mass 
that  is  rotating  about  its  axis  of  symmetry.  The 
argument  for  this  interpretation  comes  from  3  separate 
places ; 

(1)  the  metric  appears  to  have  the  appropriate  multipole 
structure  when  analyzed  in  the  manner  discussed  earlier; 

(2)  in  a  covariantly  defined  flat  space  associated  with 
the  metric,  the  Rieraann  tensor  has  a  circular  singularity; 

(3)  there  exists  a  closely  analogous  solution  of  Maxwell’s 
equations  that  is  the  field  of  a  rotating  ring  of  charge." 
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CHAPTER  V 

BONDI'S  PROBLEM 


§  5-1  Introduction. 

Under  "normal"  conditions,  Newtonian  theory  is  a  remarkably  close 
approximation  to  general  relativity,  apart  from  one  or  two  minor  points 
(advance  of  perihelion  of  Mercury  etc.).  Indeed,  our  discussion  up  to  now 
has  relied  quite  heavily  on  the  analogy  between  the  Newtonian  and  relativistic 
theories.  On  the  debit  side,  the  predictions  of  Einstein's  theory  under 
"normal"  conditions  can  lead  only  to  results  which  have  been  essentially 
familiar  for  250  years.  The  most  interesting  applications  of  general 
relativity  are  to  those  "exotic"  situations  which  classical  theory  is 
inadequate  or  powerless  to  deal  with.  It  is  largely  for  this  reason  that 
fields  remote  from  present  day  experimental  varif ication,  such  as  gravita¬ 
tional  waves  and  cosmology,  are  so  heavily  cultivated  by  research  workers. 

Our  object  in  this  chapter  is  to  deal  with  one  such  "exotic"  application. 

In  an  essay  which  was  awarded  the  first  Gravity  Research  Foundation 
Prize  in  1957;  Hermann  Bondi  [7]  considered  the  following  problem.  Two 
free  particles  of  masses  +m,  -m  are  simultaneously  started  from  rest. 

What  happens? 

The  conclusions  of  a  naive  Newtonian  analysis  of  this  problem 
are  rather  startling.  According  to  Newton's  law  of  gravitation  if  r^(t)  , 
r^(t)  are  the  position  vectors  of  the  negative  and  positive  masses,  and 
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— r  — r 

if  F^,  F0  are  the  gravitational  forces  experienced  by  them,  then  we  have 


(5.1) 


_  _  (-■”)(-»")  /■?  _  ->  v 

->  ,3  K  2  1' 


|r2"rl 


Further,  according  to  the  second  law  of  motion, 


(5.2) 


V  v  - >  y  v 

^  =  (-m)  rx  ,  Fg  =  (+m)  r. 


We  thus  find 


(5.3) 


=  r. 


m 


lr2'ri 


The  negative  and  positive  masses  are  equally  accelerated  in  the  SAME  direction 

-  the  negative  mass  chases  the  positive  mass. 


A  little  critical  reflection  shows,  however,  that  the  above 

argument  has  no  secure  foundation.  We  have  assumed  that  the  constants 

(+m),  ( -m)  entering  (5*1)  are  the  same  as  those  entering  (5.2).  The  only 

support  for  this  assumption  is  observational  (the  Eotvos  experiment)  and 

observation  is  (at  least  so  far)  confined  to  positive  mass.  We  have  failed 

to  distinguish  between  gravitational  and  inertial  mass.  Inertial  mass  is 

the  constant  m„  characteristic  of  a  particle  which  measures  its 

inert 

kinematic  reaction  to  a  given  applied  force: 


(5.4) 


F  = 


HI,  r 

inert 


Gravitational  mass  is  the  constant  (analogous  to  charge  in  electrostatics) 


.  ~ 
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entering  Newton’s  law  of  gravitation.  It  is  in  fact  possible  to  distinguish 
two  kinds  of  gravitational  mass.  Active  gravitational  mass  of  a  particle  is 
the  constant  which  measures  the  gravitational  field  generated  by  it,  while 
its  passive  gravitational  mass  determines  the  force  exerted  on  it  in  a 
given  external  gravitational  field.  Thus! 


Gravitational  force  on  particle  2  due  to  particle  1 


(5.5) 


(mx) 


(m2) 


active. gray  pass, gray  /-»  -*  \ 

3  r2  "  rl' 


lr2  "  ri 


The  situation  is  now  much  less  clear-cut.  It  is,  for  instance,  possible  - 

3 

to  mention  just  one  of  2  possibilities  -  that  a  particle  with  m 

J  r  r  active. grav 

negative  has  positive  m  and  positive  m.  _  .  To  speculate 

pass, grav.  r  inert. 

further  along  these  lines  seems  profitless. 


Let  us  now  turn  to  consideration  of  Bondi's  problem  from  a  pseudo- 
relativistic  point  of  view.  Built  into  the  foundations  of  Einstein's  theory 
is  the  principle  of  equivalence.  For  our  purposes,  this  may  be  stated  in 
the  following  restricted  form;  All  free  particles  at  a  given  point  of  a 
gravitational  field  fall  with  the  same  acceleration.  According  to  (5-^)  and 
(5.5)  this  is  only  possible  if  the  passive  gravitational  mass  of  every 
particle  is  equal  to  its  inertial  mass.  We  find 


(5.6) 


(im.  ) 

1 act. grav. 

I-*  — *  1 3 

lr2-rli 


(r2  '  rl) 


The  arbitrariness  and  uncertainty  inherent  in  the  Newtonian  argument  has  been 
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eliminated  by  the  principle  of  equivalence. 

In  considering  the  motion  of  a  set  of  particles  under  their 
mutual  gravitational  influence,  the  only  type  of  mass  which  needs  to  be 

considered  is  active  gravitational  mass.  Particles  with  positive  active 
gravitational  mass  attract  all  other  particles;  particles  with  negative 

active  gravitational  mass  repel  all  other  particles. 

If  we  reconsider  the  problem  of  two  particles  with  equal  and 
opposite  active  gravitational  masses  from  the  present  point  of  view,  we 
are  led  to  precisely  the  same  conclusion  as  we  were  on  the  basis  of  our 
naive  Newtonian  argument.  If  the  particles  are  started  from  rest,  it 
appears  that  they  will  move  as  a  rigid,  uniformly  accelerated  unit  -  a 
"Bondi  dipole"  -  with  the  acceleration  directed  from  the  negative  to  the 
positive  mass. 

The  conclusion  remains  startling,  and  it  would  be  desirable  to 
have  confirmation  by  the  construction  of  an  exact  solution.  Do  Einstein’s 
field  equations  admit  a  vacuum  solution  which  represents  a  Bondi  dipole 

and  which  is  free  of  singularities  except  along  the  world-lines  C  ,  C+ 
of  the  two  opposite  masses  (Fig,  b)?  To  settle  this  question  would  appear 
to  be  a  straightforward  matter,  since  the  field  is  static  and  axially 
symmetric  when  viewed  from  a  co=moving  uniformly  accelerated  frame.  How¬ 
ever,  it  is  at  present  still  an  open  question  whether  a  global  solution 
of  this  type  exists.  There  is  indeed  no  difficulty  in  obtaining  a  vacuum 
line-element  whose  domain  of  validity  is  the  sector  A  0  B  of  space-time 
bounded  by  the  light-like  asymptotes  to  C  ,  C  j  what  remains  undecided 
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is  whether  it  is  possible  to  extend  this  line-element  smoothly  across 
the  null-barrier  into  the  rest  of  space-time. 

Bondi  was  however  able  to  prove  the  global  existence  of  another 
type  of  solution  involving  four  particles.  This  solution  possesses 
reflexional  symmetry  and  represents  two  identical  Bondi  dipoles  moving  in 
radially  opposite  directions.  Bondi's  result  was  purely  an  existence 
theorem.  He  did  not  give  the  solution  explicitly.  This  is  accomplished 
in  §  5.3  of  this  chapter. 

One  of  the  basic  objectives  of  this  chapter  is  to  disclose  a  remarkable 
connection  between  Bondi's  extendability  problem  and  a  fundamental  unsolved 
problem  [48]  of  general  relativity:  the  extendability  of  the  Schwarzschild 
manifold  of  a  particle  in  the  presence  of  an  external  gravitational  field. 
Penetrating  Bondi's  null  barrier  is  tantamount  to  penetrating  the  "Schwarzschild 

singularity"  of  a  particle  with  infinite  mass. 

The  procedure  is  readily  understood  in  terms  of  a  semi-Newtonian 

analogue , 


Consider  a  system  of  three  collinear  particles  with  masses 
m©,  m^,  m^  .  If  one  of  the  outer  masses,  m^  >  0  ,  and  its  position  are 
arbitrarily  assigned,  the  masses  and  locations  of  the  other  two  particles 
can  be  adjusted  (in  infinitely  many  ways)  so  that  the  system  is  in  (unstable) 
equilibrium  and  static.  We  then  have,  necessarily,  m^  <  0  ,  m^  >  0  .  Now 
let  the  mass  m^  and  its  distance  from  m^,  m^  simultaneously  become 
infinite  in  such  a  manner  that  its  limiting  field  is  homogeneous.  Since  a 
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homogeneous  gravitational  field  can  be  transformed  away  by  changing  to  a 
uniformly  accelerated  frame  (c.f.  Chapter  III,  p.  39)  we  expect  to  be  left 
with  a  uniformly  accelerated  system  of  two  particles  of  opposite  mass. 

Bondi's  null  barrier  corresponds  precisely  to  the  history  of  the  Schwarzschild 
sphere  of  the  mass  m 

o 

The  plan  of  this  chapter  is  as  follows.  In  §  5*2  ,  we  deal  with 
the  problem  of  extending  the  Schwarzschild  manifold  of  a  particle  imbedded  in 
a  static,  axially  symmetric  field,  with  particular  reference  to  the  case  of 
n  collinear  particles. 

Slightly  generalizing  a  result  of  Mysak  and  Szekeres  [48],  we  show 
that  a  trivial  modification  of  Kruskal's  transformation  [54]  yields  an 
analytic  extension  of  the  manifold.  The  reflexional  symmetry  of  Kruskal's 
transformation  entails  a  duplication  of  all  singularities  in  the  extended 
manifold . 


Finally,  in  §  5«3  >  we  carry  out  the  limiting  process  just  described 
for  the  extended  line-element  representing  three  collinear  particles.  In 
the  limiting  manifold,  Kruskal's  co-ordinates  are  asymptotically  Galilean. 

The  reflexional  symmetry  means  that  the  limiting  manifold  represents  a 

solution  of  Bondi's  4-body  (rather  than  his  2-body)  problem. 

If  it  is  true  that  all  simply  connected  analytic  extensions  of 

a  given  analytic  manifold  are  isometric,  then  our  result  points  to  the 

conclusion  that  there  is  no  solution  of  Bondi's  2-body  problem  which  is 

analytic  everywhere  outside  the  world-lines  C  ,  C  .  A  uniqueness  theorem 

-  + 

of  this  type  has  never  been  proved,  although  the  isometry  of  the  various 
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analytic  extensions  of  the  Schwarzschild  metric  obtained  independently 
by  Lemaitre  [50],  Robertson  [51],  Finkelstein  [52],  Fronsdal  [53]  and 
Kruskal  [54]  provides  empirical  eveidence  for  uniqueness,  at  least  in 
the  spherically  symmetric  case.  In  any  events,  since  the  history  of  the 
Schwarzschild  sphere  is  a  null  3-sPace  -  i.e.,  a  characteristic  hyper¬ 
surface  of  the  field  equations  -  there  is  no  compelling  reason  to  insist 
on  analytic  extensions.  Relaxation  of  the  analyticity  requirement  would 
admit  a  great  variety  of  extensions,  and  perhaps  lead  to  a  satisfactory 
global  solution  of  Bondi's  2-body  problem. 

The  presentation  in  this  chapter  is  based  in  part  on  a  paper  by 
Israel  and  Khan  [59]  which  is  due  to  appear  shortly  in  Nuovo  Cimento,  In 
a  very  recent  (Feb.  1964)  paper,  Bonnor  and  Swaminarayana  [24],  attacking  the 
problem  from  a  different  point  of  view,  have  displayed  a  global  solution  of 
Bondi's  problem  representing  four  Curzon-type  particles.  Their  solution  is 
included  in  the  more  general  one  given  in  §  5*3* 

§  5.2  Analytic  extension  of  the  Weyl  manifold  representing  a  Schwarzschild 

particle  in  an  external  field. 

To  pass  from  Weyl's  static  line-element  9  for  three  collinear 
masses  [Chapter  III  page  59]  to  a  global  solution  of  Bondi's  problem,  it  will 
be  necessary  to  extend  9  within  the  Schwarzschild  sphere  of  one  of  the 
positive  masses. 

We  shall  begin  by  considering  the  problem  of  extension  for  the 
general  case  of  a  Schwarzschild  particle  at  rest  in  an  external  static, 


'V.  .w;.f 
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axially  symmetric  field.  If  the  external  field  is  restricted  to  be  regular 
everywhere  on  the  axis  of  symmetry,  then  as  shown  by  Mysak  and  Szekeres [48], 
a  direct  application  of  Kruskal's  transformation  yields  an  analytic  extension 
of  the  manifold.  Our  interest  is  in  the  situation  where  this  restriction 
is  dropped. 

We  recall  from  Chapter  II  that  every  static,  axially  symmetric 
Einstein  vacuum  metric  can  be  reduced  to  the  canonical  form 


,c  ,  2(v-A)  /  2  ,2.  2  -2*  ,,-2  2 A  2 

(5*7)  $  =  e  v  '  (dr  +  dz  )  +  r  e  d9  -  e  dt  , 


where  the  functions  A(r,z),  v(r,z)  satisfy 


Sv 

Sr 


r/S^\2  /SAv2t 
r  1  > 


(5.8) 


Sv  p  SA  SA 

S  =  2r  37  '  'SZ  ■ 


The  integrability  condition  of  equations  (5*8), 


(5*9) 


_  S2A 

V  A  =  — -  +  —  -c-  + 


Sr 


1  SA 
r  Sr 


s2a 

Sz2 


=  0  , 


is  formally  the  Newtonian  potential  equation. 


In  (5*7)  let  us  write 


(5*10) 


A  =  A  +  11/  ,  v  =  v  +  &  , 

o  ’  00  7 


where 


‘  !Vi  O  ;  I  :  ..  ■  . ;  j)  ;■  ■  ,  t  v  ,  ;i,  } 

•  ■  •  ■  ■  .  1  ■  ■  .  ■ 


. 


. 


-  96  - 


(5.11) 


A 

o 


ilog 


+  p; 

+  po 


b 

o 


+  b 

o 


represents  the  field  of  a  Schwarzschild  particle  of  mass  m  >  0  ,  v 

o  '  oo 

(computed  from  A^  by  (5.8))  is  given  by  an  expression  of  the  form  (3.61) 
(where  m^  =  ^b^)  ,  and  ip(r,z),  5(r,z)  represents  the  external  field  and 
its  interaction  with  the  particle.  Since  v^q(0,z)  =  0  for  |z-aQ|>  ^bQ  , 
elementary  flatness  requires 

(5.12)  6(0, z)  =  0 


for  portions  of  the  z-axis  not  occupied  by  matter.  On  the  segment 

|z-a  I  <  4b  ,  Tp  is  regular  and 
1  o'  2  o  > 

8A  8Aq 

limit  r  •  =  0  ,  limit  r  =  1  , 

r  -*  o  r  —>  o 


From  (5.8)  it  then  follows  that  6  is  well-behaved  on  this  segment  and  that 

s(s-g n  ,  s(g-gj 'i  ,  0 

or  oz 


(5.13) 


i.e.  6=2^=  constant  =  -2k  ,  (r  =  0,  lz“a0l  <  Tp) 


We  pass  to  spherical  polar  co-ordinates  R,  0,  9,  T  by  the 


trans  format  ion 
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e  *Cr  =  R(l  -  sin  0  , 

(5.14)  e  ^(z-aQ)  =  (R~|t>)  cos  0  , 

lc  rn  -k 

et  =  T  ,  b  =  b  e 
’  o 

The  metric  (5*7)  takes  the  form 


(r  d)  2(5-^  +k)  rH  bv-1  ._2  d2  ,n2,  -2i//+  2k  2  .  2_  .2 

(5.15)  *  =  e  v  r  [  ( 1  )  dR  +  R  d0  ]  +  e  r  R  sin  0  dcp 

K 

-  e2«"k>  (1  -  |)  dT2  , 

K 

and  according  to  (5-13)  ?  we  have 

(5.16)  5=2^+  2k  =0  when  R  =  b 

Now  define  p,  h,  t  by  Kruskal's  formula  [54]  (c.f.  Chapter  III, 

page  43) 

(5.17)  1  +  P  =  f  ,  h2  =  4peP  , 


t  =  bh  cosh  (■  — ) 
2b 


r  =  bh  sinh  ( — -) 
2b 


The  transformation  R,  T  -»  £  ,  r 
(5.18)  4  =  (l+p)"1  iP  e2,*-k) 

+  — p“ - =”pT  [' 

4b  p(l+p)e  P 


yields  (for  details  see  Appendix  X) 
(d£2-dx2)  +  e2(°“^+k)  R2(d02+e2°sin20  dcp2 

2(6-2*ak)_1]  m  _  TdT)2  _ 


) 
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Because  of  (5.16)  there  is  no  singularity  in  this  line-element  at  the 
Schwarzschild  radius  R  =  b,  p  =  0  .  We  assume  ip(R,6),  5(R,@)  extended 
as  analytic  functions  of  their  arguments  within  the  Schwarzschild  sphere. 

Tine  domain  of  validity  of  ( 5 . 18 )  is  certainly  contained  in 

p 

p  p  hh 

(5.19)  r  -  T  >  -  *  o  <  e  <  k  ,  o  <  cp  <  2« 


owing  to  the  singularity  at  p  =  -1  ,  R  =  0  .  This  singularity  will  be 
the  only  one  encountered  in  the  interior  of  the  Schwarzschild  sphere 
(R  <  b)  if  the  analytically  extended  if/( R,@)  and  S(R ,0)  are  real  and 
regular  in  this  region,  and  if  the  elementary  flatness  condition 

(5.20)  &(R,0)  =  &(R,jt)  =  0  (0  <  R  <  b ) 

is  satisfied. 

Since  (5.18)  is  invariant  under  the  transformations  £  — » and 
T  -T  ,  we  have  complete  reflexional  symmetry  in  every  2-space  Q  =  constant., 
cp  =  constant. 

The  line-element  (5.15)  is  now  to  ber  regarded  as  describing 
2  2 

the  submanifold  £  -  t  >  0  of  (5. 19)  }  and  is  connected  with  (5.I8) 

through  the  transformation  (£,t)  (R,T)  defined  by 

(5.21)  R  =  b(l+p),  Ub2peP  =  f  -  t2  , 


tanh 


i 


(5.22) 


if  -  T2  >  0)  . 
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We  proceed  to  consider  the  specific  form  assumed  by  (5.I8)  in 
the  case  where  the  external  field  is  produced  by  two  particles,  with 
masses  m^  <  0  ,  m^  >  0  at  rest  on  the  z-axis.  Owing  to  reflexional 
symmetry,  the  2~space  9  =  0  exhibits  four  singular  curves  of  constant 
R  ,  in  addition  to  the  singularity  at  R  =  0  (Fig  6). 


Setting  n  =  3  in  (3*52)  and  (3. 55),  we  obtain  from  (5.10) 


(5o23) 


t  =  \  +  \  , 


5  = 


i=l  j=l 


)  v„  .  +  2(v  ,  +  v  ) 

Z_j  ol  02 


where 


(5.24) 


V  =  •“ 

1  b  „ 
1 


log 


p .  +  p : 
1  1 


p .  +  p !  +  b . 
1  x  1 


(R,0) 


It  can  easily  be  verified  that  the  analytically  extended  functions 
,  pg  '  (R,0)  are  real  and  positive  throughout  the  region  (5.19)> 


and  that 


(5.25) 


"k  D  ,  2  G 

e  p  =  R  -  b  cos  ~ 
~  o  2 


R 


.  .  2  9 

b  sin  -  , 


where  k  is  computed  from  (5. 13)  and  (5.24)  to  be 


(5.26) 


log 


(zo  - 


Zo) 


-  z. 


(For  details  see  Appendix  Xl)„ 
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It  follows  that  ip(R,Q)  and  £>(R,0)  are  real  and  regular  in 
the  region  0  <  R  <  b  ,  and  (5. 20)  is  satisfied.  Thus  the  only  true 
singularities  produced  in  the  manifold  by  the  Schwarzschild  particle  bQ 
are  a  pair  of  curves  with  equation  R  =  0  .  Similar  remarks  apply  to 
the  other  two  particles. 

§  5-3  Transition  to  Bondi's  problem. 

We  preface  this  section  by  indicating  how  a  "homogeneous 
gravitational  field"  (or  equivalently,  a  uniformly  accelerated  reference 
frame)  can  be  characterized  as  a  limit  of  the  Schwarzschild  solution. 

In  Chapter  III  (page  39)  we  saw  how,  starting  from  the  Minkowski 
line-element  of  flat  space-time, 

(5.27)  d^2  +  dr)2  +  d£2  -  dr2  , 

the  co-ordinate  transformation 


(5.28) 


^  =  r^  cos  Cp,  r)  =  r^  sin  <p 

t  t 

£  =  z^  sinh  (“)  ,  t  =  z^  cosh  (— )  (3  =  constant) 


leads  to  the  form 

(5.29)  =  dr*  +  r2  dtp2  +  dz2 


2 


Since  each  point  with  fixed  space  co-ordinates  r^  ,  z^,  <p  executes 
hyperbolic  motion,  0  ^  may  be  interpreted  as  the  metric  of  flat  space- 
time  viewed  from  a  uniformly  accelerated  reference  frame.  The  motion 
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of  this  frame  is  not  affected  by  a  conformal  transformation  r#  ,  z.  r  , 
z  of  the  space  co-ordinates,  defined  by 


(5.50) 


20r  = 


r*  2* 


2(3z  =  J(z*  -  r;) 


This  has  the  effect  of  casting  cp  into  the  Weyl  canonical  form 

ua  J 

(5.7)  ,  v;ith  A  =  A  ,  v  =  v  given  by 
'  '  '  ’  ua  ’  ua  J 


(5.31) 


ua 


(p  +O 

1 1  00 

=  flog  - 

2P 


_  1 


ua 


iiog 


p  +z 

o  o 

2p 


where 


2  2  2 

p  =  r  +  z  ,z=z 
o  o  '  o 


If,  for  comparison,  we  write  down  the  Schwarzschild  spherically 
symmetric  solution  $  for  a  single  mass  extending  from  z  =  0  to 

z  -  -b  in  the  Euclidean  map  ofWsyl's  co-ordinates,  we  find  A  =  A^  , 
v  =  v  where  (c.f.  (3*21)  with  m  =  Ab) 


(5.32) 


p  +  p'  -  b 

-v  li  o  o  o 

A  =  flog  - ■ — -  ,  v 

p  +  p'  +  b 
o  o  o 
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=  ilog 


P  D  1 
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P  o' 
o  o 


Now  make  the  constant  scale  transformation  $  „  ,  -»  ,  ,  defined  by 

schw .  Schw .  * 


8ij  ^sij  =  £®ij  -  «uu  -*8fo  - 


-1 

€  g 
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then  $  I  ,  is  also  a  solution  of  the  Einstein  field  equations.  Further, 
Schw . 
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b  =  e”1^ 
o 

and  let  e  ->  0  keeping  (3  fixed.  Since 


limit  (A 

o 

e  ->  o 


ilog  e )  =  A  ,  limit  v  =  v  , 
d  ua  '  oo  ua  7 

€  ->  O 


we  obtain 

(5.33) 


limit 
€  -»  o 


4> ' 

Schw. 


4> 

ua 


Yet  another  procedure  would  be  to  throw  ^gc^w  into  Kruskal's 
form  [Chapter  III,  page  1+3]  before  taking  the  limit  e  -►  0  .  This 
would  lead  directly  to  the  Minkowski  form  (5,27).  We  shall  use  a 
completely  analogous  method  to  transform  the  3-particle  solution  obtained 
in  §  5.2. 


In  (5.15)  and  (5.18)  let  us  set 
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r  =  e  r 


z  =  e  z 


a.  = 
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€  V 


m.  =  €  m. 

X  X 


where  e  is  a  small  parameter.  Then  V  and  v  (i  =  1,2)  are  given 
by  the  same  formal  expressions  in  the  barred  as  in  the  original  unbarred 
quantities.  We  shall  further  write 


k|32p 


2  0 
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so  that,  according  to  (3.15), 
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(2P?)2  =  r2  p*  +  0(e)  , 

2pz  =  J(p*  -  r2)  +  0(e) 


In  place  of  0  ,  we  shall  introduce  as  co-ordinate  in  the  line- 

element  (5.18). 

We  keep  (3  and  all  barred  parameters  fixed  and  let  e  -»  0  , 
Then  k  -» 0  by  virtue  of  (5.2 6).  From  (5.I8)  we  obtain  a  line-element 
which  continues  to  satisfy  the  Einstein  field  equations: 

(5.3*0  ®  =  e^(d £2  -  dT2)  +  e2^°"^^(dr2  +  e"2&  r2  dq>2) 

+  e^[e2^-l] 

r  - 

where 

(5.35)  &  =  5  (r2  ,  p*)  ,  f  =  f(r 2  ,  p*)  , 

and 

(5.36)  P*  =  i2  -  r2  , 

o 

according  to  (5.21).  Explicitly,  we  have  (dropping  all  bars) 

(5.37)  f  =  \  +  \  , 

2  2 

(5.38)  8  '  Z  Z  Vlj  +  2(v01  +  v02)  ’ 

1=1  j-1 


TdT  )2 

2 


'  5  1  ■'  •  -  »  -v:  ,  -T  :V‘  „  r:  .  :  ?  /,  f 

■ 


••  ^  £j  •„* ,  ■;  /-;b  j  ....  v,-".  C  ;■ 

... 


"<  ■  "  '  '  "  '  :s.  . 


- 


' 
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with 


m. 


sr 


log 


pi  +  P1  -  bi 

P1  +  Pi  +  b± 


V 


ij 


Ti 

bibj 


log 


E(i, j)  E(i',j') 


and 


m. 


(5.39) 

V  .  = 
01 

=  -  log 

2bi 

(5.40) 

o0) 

pi 

=  [(2S)': 

(5.41) 

z(') 

zi 

=  (40)'1 

(5.42) 

23  Z 

0 

-  -  r' 

E(0, i 1 )  Pi  +  2i 

E(0,i' )  p!  +  z!  ’ 

r*  P*  +  z|  (1=1,2) 

(P*  "  r*)  -  (+)ihi  (1=1,2) 
,  2Ep0  =  |(p*  +  r2)  . 


2  2 

Since  6  and  ip  tend  to  zero  when  £  +  r*  for  fixed  t  , 

the  co-ordinates  r^  ,  cp,  £  are  asymptotically  cylindrical.  We  recall 
from  (5.13)  that 


b  -  2y  =  0  (p*)  Cp*  -*■  0)  , 

2  2 

so  that  there  is  no  singularity  on  the  null  cone  £  -  t  =  0  .  In  fact, 

the  only  singularities  involved  in  (5.3*+)  are  the  histories  of  the  two 
"rods"  on  the  z-axis: 


r*  =  0 ,  1  (its)-1  a2  -  2) 


|  <  HbJ  (i  =  1,2,). 


(5.43) 


-  a, 
1 


'  '  w 


.  ,  j- r/"1':''  *  .'.nvt 

;b  r.  (?;i  ;}  .  - 

■  -  ■  •  '  r  ■  :  ■  '  .  ! 

r 
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One  checks  immediately  that  the  elementary  flatness  condition  (6=0 
when  =  0)  is  satisfied  everywhere  outside  these  singular  regions, 
provided 


(5.1*1.) 


m. 


log 


ai  +  2^ 


2-i  » 

=  r  log 


'o  - 


2m  ira  2 
"^2 


log 


1  - 


bl  b2 


( a i -ag)  -te>1-b£)' 


This  condition  follows  if  we  write  the  elementary-flatness 
condition  (5.70)>  for  three  particles,  numbered  0,  1,  2  and  let 
mQ  =  -|b  =  -a  ->-oo  (for  details  see  Appendix  XII). 

The  arguments  of  §5.2  show  that  each  of  the  singularities  (5.4j) 
is  removable  by  a  co-ordinate  transformation.  The  only  true  singularities 
are  the  four  time- like  curves 


r*  =0,  £2  -  t2  =  4(3a.  (i  =  1,2)  . 


It  is  now  clear  that  (5. 54)  represents  an  explicit  analytic 
solution  of  the  Bondi  4-body  problem  described  in  §5.2.  We  observe  from 
(5.44)  that  the  masses  of  the  two  components  of  Bondi  dipole  are  not 
precisely  equal  and  opposite,  but  proportional  to  the  inverse  square  root 
of  the  respective  accelerations  (as  measured  by  the  curvature  of  the 
world- lines) . 


. 


■ 

-  .  '  J  ''  ,  . 

■  -  '  -  :  :>r 
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APPENDIX  I 

CALCULATION  OF  RICCI  TENSOR  CORRESPONDING 
TO  THE  METRIC  FORM  (2.33). 

Following  are  the  detailed  calculations  of  the  components  of  the 
Ricci  tensor  corresponding  to  the  line-element  (2.33)  (in  Chapter  II  page  22) 


We  have 


(1) 


a  2  A2r  t  j  1\2  /.  2\2 ,  t»2/j  3\2  _2  /  ,  4\2 

ds  =  A  [(dx  )  +  (dx  )  ]  +  B  (dx  )  -  C  (dx  ) 


1  2 


where  A,  B,  and  C  are  functions  of  (x  ,x  )  only. 


(2) 


811 

S22  A  * 

S33 

a 

and 

gbb 

11 

g 

22  -2  . 

=  g  =  A 

s33  = 

b"2 

and 

bb 

g 

R.  .  =  R. 
ij  ijk 


Ricci  tensor  is  defined  by 

(3) 

Explicitly  we  have 

(1+)  R.  . 

ij 


g^1"1  R,  .  .  =  R  . .  . 

kijh  Jl 


ra  _  ra  +  ra  rb  -  ra  rb 

ai,j  ij,a  +  bi  aj  ij  ab 


a  „b 


'  iUog(-g)](1J  -  rlj)3  -  ir1j[iog(-g)]^a  +  rbl  raj 


■  ■  i  . 

'■  •  ' '  ■;  v 

Iff,;-,  iuwi  Lv..v:  :  :v  e.4.)  j  v; 


\  :V  ■"  ,  '  '• 'S  0  {  '  \  ,  :  .  ..  '■  , 


.  <  ^vrfv  ■■■ 


^  ” 

,  . 

* 

\r4 
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The  non-vanishing  components  of  Ricci  tensor  corresponding  to 
(l)  are  as  follows. 

Rn  -  Vx  V>  +  Vi  s2a)  +  (i  4®  + 1  4C> 

+  X  S2A(5  32®  +  5  32C>  -  I  31®  +  £ 

®22  =  31(A  dlA)  +  32^A  32A-1  +  3 22®  +  C  322C^ 

+  X  31A(|  31®  +  h  31C)  -  I  32A(I  32B  +  h  32C)  • 

(5)  R12  =  1  sf2B  +  i  4c  -  i  s2A(i  axB  +  i  axc) 

-  X  31A<I  32®  +  £  32c) 

R33  -  \  (311B  +  4®  +  ^31B31C  +  32B32C)] 


\k  =  $  [311C+  4c  +  f<3l®  31C  +  32®  32c>] 

A 


Since 


g1J  R.  .  =  R* 

ij  i 


Rj  =  g55  r33  =  b'2  r33 


(6) 


■  'ft  --  ■  *  .  “•  '  ■ 

.  :  •  t  K-,  ■  :  f ) 


":Yr  F  . *!>  y  '  YY  -  U,6  Y/  -  r:: 


. 


<3  .. 

"■  /  '  '  "l 


r¥;:*  ::  y  7 


•  '  -  .  .  ]z: 


:  t; 


-  io8  - 


(7) 


M  44 


-2 


\  ~  8  R44  "  “C  R44 


Hence 


(8) 


R5  +  R4 


^  33 

'  B2 


R_  R 


44 

2 


-7T-  [&?,  (BC)  +  a2  (BC)  ] 

A  BC 


In  a  matter  free  domain,  the  field  equations  are 


Hence  (8)  becomes 


R.  .  =  0 

ij 


R3  +  <  -  ^  '4  <BC>  +  4  <BC”  ■  0 


1  2. 


i.e.  (BC)  is  a  harmonic  function  of  (x  ,x  )  . 


•s _ -s _ 

dl  ~  1  J  2  =  .  2  * 

ox  ox 


a2  -  a2  -  4 

11  s  2  ’  22  '  i  2 

ox^  dx^ 


I  '  M* 

' 

l 


■  •  '  ■  -  .  ■  ..  c  ■  ■ .  .  .  r- 


■ 


■  ‘ 
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APPENDIX  II 

CO-ORDINATE  TRANSFORMATION  OF 
THE  METRIC  (3.1). 

Following  are  the  detailed  calculations  for  the  co-ordinate 
transformation  of  the  line-element  ( 3 • 1 )  (c.f.  Chapter  III,  p.1+3)* 

From  (3.I),  we  have 

/,x  _  2(v-A),,  2  ,  2\  2-2^  2  2h  2 

(1)  cp  =  e  '  '(dr  +dz  )+r  e  dcp  -  e  dt 

where 

(2)  *=£  log  p  tPV.-.b 

p  +  p'  +  b 

/  \  2m2  ,  (p  +  p 1  )2-b2 

(3)  v  =  -3-  log  ^ ^ - 

b^  4pp ' 

/|  \  2  2  2  12  2  /  r\2 

(4)  P  =  r  +  z  ,p  =  r  +  (z-b) 

Set 

(5)  r  =  R(1  -  |)i  sin  9 

(6)  z  -  -l-b  =  (R  -  ib)  cos  0 


(7) 


m 


"  ■  '  .  •  His  /  :  .  ..  1 


■ 


3b  •  . 
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Hence 


(8)  dr  =  sln.9.(‘jR~b)  dR  +  R(x  _  |)i  cos  g  dg 

2R(1-  §)?  R 


(9)  dz  =  cos  0  dR  -  (R  -  ib)  sin  0  d 0 


Squaring  and  adding  (8)  and  (9)  we  have 

?  b2 

0  0  R  -  bR  +  4-  sin  0  r  _2 

(10)  (dr2  +  dz2)  =  - ^ - 


R 


dR“  D2,fl2 

+  R  d© 


L  h-!> 


From  (2)  and  (7)  we  get 


(11)  A  =  ilog 


p  +  P1  -  b 

p  +  p'  +  b 


(12) 


v  =  -Uog 


(p  +  p* )2  -  b' 

W 


From  (5)  and  (6) 


(13) 


2  2 

p2  =  (R2  -  bR  +  )  +  ^r-cos2©  +  b  cos  0(R  -  -i-b) 


(R  -  -i-b)2  +  cos2©  +  b  cos  ©(R  -  -i-b) 


( 14 )  p '  2=  (R  -  ^-b)^  +  —  cos2©  -  b  cos  ©(R  -  |-b) 


2  ’2  lu \2  b  2_s2  ,2  2_,_  lu v2 

pp  =  [(R  -  -±b)  +-j—  cos  0  J  -  b  cos  0  (R  -  ^-b) 


'V  ’>  :  ''v';  ,  :  ■  ;  '  V’A  ■: 


■  - 

. • 


' 


■ 


; ; 

■•'i  '■  1  '  1  1  ■ 


.  fr. .  „  :  .  ..  ; 


Ill 


(15)  p^p  =  { (R  -  £b)2  -  cos'-  e}2 

2 

(16)  pp'  =  {(R  -  ^b)"  -  cos20} 

From  (13),  (14)  and  (16),  we  have 

(p+p')2  =  b(R  -  |)2 

(17) 

(p+p')  =  2(R  -  |)  . 


Hence  (ll)  and  (12)  can  be  written  as 


(18) 


 p+p 1 


p+p'  +  b 


(1  -|) 


(19) 


2v  _  (p  4-  o'f  -  b‘ 

Itpp' 


RS(i 


>2  ,  _  b  .  2  -v 

(R  -  bR  +  -r—  sin  9) 


From  (l8)  and  (19) 


(20) 


2(v-A) 


R 


2 

(R2  -  bR  +  y-  sin20) 


Hence  from  (10),  (l8)  and  (20 ) 


?a|cr 


'  -  '  -t-  ;• 

. r  " 

-  .  . . . 


) 
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(21) 

and 

(22) 

Hence 


2(v-*),.  2,2,  ( 

e  (dr  +  dz  )  =  ( 


dR 


(i--) 

'  R' 


P  P 
+  R'  69  ) 


2  -2A 
r  e 


R 


sin 


e 


(i)  can  be  written  as 


e2(v-A) 


+  dz2) 


2  -2 A 
+  r  e 


dR 


d-1) 


+ 


P  P  P 

R  (d©  +  sin  0 


dcp2) 


(1  - 


POlcr4 


■  : 


, 

<  •  '  v  ■.  . 


■ 
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APPENDIX  III 

VERIFICATION  OF  EQUATION  (3.37). 


Following  are  the  detailed  calculations  for  the  verification  of 
equations  (3*37a)  and  (3<37b)  (c.f.  Chapter  III,  page  46), 


Using  the  Levi-Civita  equations 


(1) 


dv 

dr 


=  r 


61 :) 


d4x\ 

dz  ) 


2 


(2) 


dv  _  „  dA  dA 

dz  r  dr  dz 


we  write 


(3) 


dv 


dv 


=  r 


r  /dA 

dr 


=  r 


dA 

-  V  dr 


dA\  1  dr  0  dA  dA  dz 

d^jj^  +  2rd7  dl  df 


dA\  dr  dA  dA  dz 

dzy  J  dr)  +  X  dr  dz  dq 


or 


(4) 


s’  ■■[(*§ 


=  r 


v  ^  +  Tv  a 


dzy)_ 


dr 


+  2r 


and 


---  ..t  .  j  ,  .v? 

:  '■  •  ■  :  . : 1  ■  ,  ? )  ■  «v:  U;r.*  .  ■,  >  i,;  ^ 


■  r 


.... 


■ 
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From 

(6)  r  =  b(£2  +  l)i(l  -  T,2)i  ,  2  =  b  i  T!  , 


we  have 


(7) 


dr  = 


b£(l  -  t i2)2 

(?2+iF 


dl 


.  bn(g  +  1) 
o  1 


p  1  drl 

(  1  -  T]  )2 


and 


(8) 


dz  =  bt]  d^  +  dr]  . 


Hence  from  (7)  and  (8) 


(9) 


8r  ,  /.  ( 1  -  ) 2  8r 

s?  =  b|  ^ 


b  ($2  +  1)2 
-  br)  ■ ~nlT 

(1  "  rf)2 


(10)  |f  “  brl  *  |f  =  ^  * 


Setting  dz  =  0  in  (8),  from  (7)  and  (8)  we  have 


' 


- 
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(11) 


dr  =  b?  (1  -  '1  d?  -  bn  v?  + 1  dn 

(e2  +  d5  a  -  i2)s 


0  =  br)d^+b^dr) 


Solving  (ll)  for  dl-  and  dr]  we  have 


dz=o 


£(g2  +  i)2_(i,  -  n2)' 

b(H2  +  n2) 


(12) 


dz=o 


•  i(g2  +  l)4l  -  n2)' 

b(?2  +  n2) 


Next  by  setting  dr  =  0  in  (8)  we  have 


(13) 


n(62  +  1) 

b(t2  +  r,2) 


1(1  --_,n2L 

b(|2  +  t,2) 


From  (4)  we  get 


dv 

si' 


d£  dz 


+ 


dr]  dr?  dz 
dr  dz 


-2rA^A 


d  T)  d$\  dz 
+  dr  3zy  d£ 


.  *  ■ ' 


;  ;  J  :  -.l  ri  v; 


' 

•  •  ■ 


>'■  >  ■  {  ■  ' 


O--X0' 


"  ■  *“•-  \i'C- 


>.  ;  .  ■  ■  «  ' 

'■■•'6  .  '  '  ,:V,  , 


.'t. 


'  1 o 

'  i  . 

•  •  X  '  ■■■  .r,  ■ 
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From  (9),  (10),  (12)  and  (13),  we  get 


(14) 


.  /^ifl  3r  _  j(f,2  +  1)(1  -  T,2)(g2  -  g^e2  -  1} 


bzJ  J  di 


,c2  2^2 

(£  +  n  ) 


(15)  2r 


i  M  Sj  .  ggng(e£  -1  1)2(1  -  I2) 
Sr  Sz  S£  ,,2  2.2 

a  +  n ) 


Hence 


(16) 


/s*\2 


"aaTi  6r  2  6^  m  6z 
a$  +  d  6r  6z  ae 


_  eo-2  ±  1 ) (1  -  42) 

,*2  2v 

(£  +  n  ) 


Also 


(17) 


2  2  >2x 

T1  -  %  ) 


(18)  2r  ^  ^  & 


Sr 


or)  oz 

3z 


2nH(i 1±,.0(A  -If 

(t2  +  2) 


Hence  from  (17)  and  (]8) 


(19) 


r[f  _  AU  6r  ,  p  67  6r]  6z 

i^ar J  \<dz )  f  as  ar  az 


.  %&—£t 

if  +  I2) 


Again  from  (9),  (10),  (12)  and  (13) 


(20) 


Si  Sa  S£.  Si^Sr  2£2ti(£2  +  1)0  -  H2)2 

(f2  +  n2)2 


z  az  dr  dry  a^ 


(  ;  ('5'-  ■" 
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and 

(21) 


+  aj  |i\dz 
Sr  Szy  cij; 


+  i) (i  -  i2)(e2  -  2^n2  -  n2) 


Hence  from  (20 )  and  (21) 


(22) 


S£  SrA  Sr 
Sr  Sry  S^ 


Sjj  S^\Sz 
Sr  cte y 


-2nd2  +  i)(i  -  n2) 

(Z2  2^ 

(£  +  4  ) 


Finally  from  (4),  (16),  (19)  and  (22) 


(23) 


Sv 

5 


(eV)  w 


$(l-n2)2  f&f 

(^2+n2)  W 


2u(i-+1l(.1.:.Ti  ..).  ^ 

(,2  2\ 

(I  +4  ) 


SA 

Sr) 


In  exactly  the  same  way  it  can  be  shown  that 


Sv  _  -nd2*!)2  /3A\2  n(l-42) (£2+l)  /7SA"f  2a(£2+l)(l-n2)  SA  SA 

^  "  ft2  2 v  IS tj  “  ,ft2  2 <.  ISrJ  /*2  2v  * 

1  (4  +4  )  X  7  U  +  4  )  X  7  (I  +4  ) 


' 


'  x’.. 

■  •  V-c  ■  - 

T&-  St  TO 


%  X  '  .  . 


•  .... 


APPENDIX  IV 


FIELD;,  OBLATE  AND  .  PROLATE  HOMOGENEOUS  SOLID  SPHEROIDS. 

With  reference  to  Chapter  III  page  '30,  following  are  the  detailed 
calculations  for  the  field  of  oblate  and  prolate  homogeneous  solid  spheroids, 

I.  Oblate  spheroidal  co-ordinates  are  defined  by 
r  =  b(l+|2)i  (1-t,2)S 
z  =  b  £  T)  . 

dr  =  -  r)(0>l)  dl) 

(r+l)2(l-tF 

dz  =  b{t]  d  ^  ^  d  t|) 

If  A  represents  the  Newtonian  potential  of  an  oblate  homogeneous  spheroid 
then 

(5)  A  =  A[4  cot"1!  +  (jjrf-l)  {(1+3!2)  cot  X!  -  3|] 

where  A  is  a  constant  connected  with  the  mass  of  the  spheroid 


(1) 

(2) 

Hence 

(3) 

(4) 


-  '  ■ . 


•t"  ■  !  ■  -  i-  .  ‘f  ■  •  '  t  .  .  .!  ■  <■»*>  xoi,  ,  • 


■  ' 


■ 
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(7) 


dv  _  dx  ax 

^"z  r  dr  ctz 


From  (6)  and  (7) 


(8) 


,  dv  ,  dv  . 

dv  =  t—  dr  +  ■?—  dz 

or  dz 


r/dXv2  ,dAx2,  dA  dA 

rC  (^7)  -  (37)  ]  dr  +  2r—  —  dz 


'dr  dz 


But 


r{  (^7)^  "  dr  can  be  written  as 


(9) 


,  [,§,= . 


dr 


=  A2r^3(l-T12)(l-5T12)  +  ed-T]2)2}  (cot"1!)2 
L 


2{|2(l»5T]2)(l-r12)}  (cot"^)  + 


1* n  A  43(i-3Ti2)(i-n2)-^T12(i-n2 


(£2+D 


dl| 


A2r(|2(|2+l)r,(l-5r,2)  +  n(l-n2)(|2+l)}  (cot'1?)2 


-  2(?n(52+1)(l-5n2))  cot  X|  +  (|2r|(l-5T]2)-4r|d 


d  T] 


and 


'  ' 
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(10) 


8A 


dz 


=  "  A2r4^T]2(|2+l)(l-T)2)(cOt"1|)2-l+Tl2(l-T12)(2^2+l)cot"1^ 

-  A2r4|2n(^2+l)(l-T)2)(cot“1^)2-4^(l-T|'L)(2^2+l)cot‘1^ 


(8)  becomes 


(ll)  dv  =  A2  |^{!3-10!3r)2  +  9^ ^  -  6|rj2  +  5^ T]^  +  !)(cot  ^)2 


+  [20^2t]2  -  l8|2^  -  2^2  +  4 T)2  +  4r)*")(cot  ^) 


Hence 


•  32  .  rt43  4 


4  in 


-  lot  n  +  9S,  n  -  8g,n  +  8^n  ) 

(?2  +  1) 


d!  + 


+  ^(9l\3  -  -  6^r)  +  10|2t^  +  rp  -  rj]  (cot"1!  )2 


+  (lO!5^  -  +  6|r]  -  l4!r)3}(cot  1|) 


+  (9$2t)3  -  5|2t]  -  4n3} 


dn 


(12)  =  a2 


!  3- io!  3  t]2 +9!  3  r/*-6!  r]2+5|'n^+! 


(cot  L!)2  + 


r,rst^-  2  ,0*2  ^  .i2  ,  2  «.  4 
20|  r)  -18!  4  ~2s  +4  9  -4r] 


(cot  li)  + 


[!3  -  lo!3!-)2  +  9^'V  -  8$ q2  +  8!ri^] 

(e2  +  i) 


+  4!ti2(1“t)2)  d! 
+  ^2T](l-92)  d9 
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and 


(13) 


dv  .2 

Tn  -  A 


9£  ^ T)3  -  -  6(-2r)  +  10 ^ t\  +  x\  - 


(cot"1^)2 


10|3r]  -  18£,3t]3  +  6£r)  -  Ibir? 


(cot"1^)  +  ^9i2rp-^i2T)+brp 


\ 


For  the  determination  of  v  from  (12),  let  us  consider  the  evaluation  of 
the  following  integrals 


(i) 


^(cot-1^)^ 


&4 

-i —  (  *■  ~  l/5’  \2  i 

k  (cot  V  +p 


,b  , -lt 

d? 

r+i 


£  (co t-\f  + 1  /  &-+  4  ■■  s2---1- *■  -1  ■  -1 

4  2  J  tr  + 1) 


cot  |  di| 


(cot-1!)2  +  i  J  fcot-H  d|  -  i  J  cot"1^  d|  +  i  J~  ^°g*  ^  d£ 


=  V(c°r4)2  +  |f  V  cot'4  +5/  d0  - 1  ^cot'4  +  fM 


e  +1 


-  rCcot-1?)2 


=  |-(coc'1d2+|  e3cot_1e-^  y^u-  |  d|  - 1  ?  cot'1^  - 1  iog(e2+D 


ifcot-^)2 


=  •i(|!*-l)(cot‘15)2  +|  |(?2-3)cot_1|  -  ilog(e2+l) 
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(ii)  In  exactly  the  same  way  we  get 


2 

I  ?(cot'1?)2d?  =  (■  g^Hcot"1!)2  +  ?  cot'1?  +i  log(?2+l) 


I2  cot"1?  d?  =  ^  cot"1?  +  |  ?2  -  |  log(?5+l) 


and 


r  -ii- 

J  h2+d 


i!  .  I 


log 


(e2  +  1) 


Hence 


(H) 


v  =  A 


r(*2+D 


^(I-IOt^t]4)  +  ( 1  ~  t]2  ) 2 


(cot-1^)2 


(l-6-n2-7r)4)  -  |2(l-10Tj2+9ri4) 


cot  ^ 


?2r  2  u 

+  tH  1  -  10  n  +  9n 


+  f(n) 


where  f ( rj)  is  an  arbitrary  function  of  t] 


Differentiating  (14)  partially  with  respect  to  t] 


■"  ■ 


T 


1 


-  123  - 


(15) 


av  2 

N-  =  A 

an 


ofi4  3  lnC2  3  ^2  3 

9q  r)  -5?  r|+10£  T]  -6t,  rj-f-rj  -r\ 


(cot_1£)^ 


L i8£  t^+io£  V  i4£  rp+6|  t) 


(cot-1!-  )  + 


9eV-5£% 


+  f 


From  (15)  and  (13)  we  have 


(16)  f'(r|)  =  bk2r?  . 


Integrating  we  get 


(17)  f  ( T) )  =  A2n4  +  a2c 


where  C  is  a  constant  of  integration, 


From  (14)  and  (17) 


(18)  v  .  a2||(|2+i) 


?2(1-1o7+9ti4)+(1-ii2)2 


(cot  ^|)2  + 


(l+6r)2-7il  ^)“^(l“10r)2+9il4) 


(cot  1|)  + 


1  -  10  T)2  +  9r\k 


+  T)  +  C 


The  constant  of  integration  C  can  be  calculated  from  the  condition  that 
the  line  element  becomes  flat  at  infinity  i.e.  v  -> 0  when  r  =  0, 
z  -*  00  corresponding  to  £  -» 00  ,  r|  =  1  .  This  condition  is  satisf  ied 
only  if  G  =  -1  „ 


(t|) 


' 


*■ 


; 


Hence, 


(19)  e‘(i-i0n2+9i|4)+(i-n2)2 

L 


(cot  ^)C  + 


2 

+  r(1+6r12-7T14)-^2(l-10T)2+9Tl4)  (cot"1^)  +  \  r  1-1°r12+9r)4 


+  9 


or 

V  =  9a2r),^(l+^2)(l+9^2)(cot‘1!;)2-(l4|+l8^5)(cot“1^)+942+U 

L 

-  l8a2r12l  (]+$2)(]+5|2)(cot'1^)2-(6^+l0^3)(cot"1|)+5^2 

L  J 

+  9a2r(l+^2f(cot“14)2+2(|-|3)cot"1^+^2-4 

L 

(where  a2  =  4A2)  . 

This  result  agrees  with  Misra's  [27]  result  of  oblate  homogeneous 

spheroid . 

II,  Prolate  spheroidal  co-ordinates  are  defined  by  (c.f.  Chapter 
III,  (3.1*2)) 

(20)  r  =  b(l-|2)i(7-l)i 

(21)  Z  =  b  4  T) 


Hence 
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dr  -  bU(l-QdTi  -  g(rt-l)  dj) 

d-rF(Ti2-iF 

dz  =  +  £dr)] 

If  A  represents  the  Newtonian  potential  of  a  prolate  homogeneous  spheroid 
then 

(24)  A  =  I<[4  coth  ^r|-(3^~l)  {(3'H^“l)coth  rj— 3r) }  ] 

where  K  is  a  constant  connected  with  the  mass  of  the  prolate  spheroid. 

As  before 


(22) 

(25) 


(25) 


,  (~/dAv2  ,dAv2 

dv  =  r'  (^r)  -  (— ) 


bz' 


dr  + 


-I 


^A 


dz 


where 


(26) 


.r (&)2 .  iSkf 

\ (br>  'dzJ 


dr 


=  k2!  13(i-e2)(i-5i2)-i,(i-2?2-)4lt) 

r  ti3(1-5£2)(1-62)+Ui£2(1-£2)1 

(t,2  -  1) 


(coth-1n)2-r2n2(i-52)(i-5?2) 


d^  - 


COth  ^  Tj 


k2  |[n2(Ti2-i)e(i-5?2)-?(Ti2-i)(i-?2) 


(cot  ^T])^ 


f2^(i-5e2)(n2-i) 

L 


(coth  %)  +  Tt]2^(1-5|2)+4^ 


df* 


' .  ■  .  ... 

;  '  -  -■  ■ 


. 


. 


;;;  ■'.>  1 

■ 

■ 


. 
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and 


(27)  2r  ^  ^  dz 


=  -K2|Ur^2(n2-l)(l-e2) 


(coth"\)2+l4|2(l-e2)(l-2r12) 

L 


coth  ^  T) 


+  U?2r,(l-|2) 


d  r) 


-K2  |U?Q2(1-|2)(t,2-1)  (coth'1n)2/lt^(l-|2)(l-2^2) 


coth  ^r) 


de  . 


From  (26),  (27)  and  (25)  we  have 


(28) 


dv  =  K2-|  i  r]^+97^^-10r)\2H-6r)^2-5T]^-r)  (coth  %)2  + 


T  ,.2  2  .. o  2.4  0  2  ,  .4  ,  4 
+  20^  r)  -lor]  ^  -2 r)  +4|  -4r] 

L 


(coth  ^  T) )  + 


+  [  x?-l0v^i2+9^ih^vt2-Qrjk]\  d 

(n2  -  1)  J 


2I  n2£  10  2*3  *3  , 


+  K  -j  9rj  |  -5r|  £+69  (*-10t]  "I 

L 


(coth  ^ t] ) 2  + 


+  aOT]^-6r,^-l8T]3t5+l4T]^3 


(coth-1!])  +|^9t]2|^-5t]2|-44^ 


- 


- 


’f  ■  '  ' 


■■  ' 


' 


.  .  • 


' 


: 
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Hence 


(29) 


dv 


-  K2  j  »i5(l+944-loe2)+ii(6?2-5?1*-l) 


(coth  %)2  + 


+  ^  n2(20i-2-i8£l+-2)+4£2(e2-i) 


(coth  ^r|  + 


and 


+  [i3(u9gu-ioe2)4geg^(i-e2)]l 

(i2  - 1)  1 


(50)  =  K2  jr9ri4i;3-5r)4?+6n2|-10Tl2e3+§3-5 


(coth  ^T))2 


+  10r)^-6^  rj-l8r]^^+l4r|^ 

L 


(coth  1r\)+[9r\2^-5T\2£-k€? 


In  order  to  determine  the  value  of  v  from  (29)  ,  we  consider  the  evaluation 
of  the  integral 


(31)  f  r|5(coth  1r])2  d r|  =  (coth-1ri)2  +  \  F 


k  ..  -1 

n  coth  ri  , 
“2 - 1  dn 

(r)  -1) 


4  2S  /  2 


=  -J-(coth-%)2+i  f  (1  -1  Ml  -1)U  coth-i 
4  J  (n2-i) 


dr) 


■ 
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4  -1  2 

(coth  r|)  +  ^ 


r }'  coth  %  dr]+|  /  coth  *r|  d r)+g  /  C-°^ - **  dr] 

J  rJ  (n  -i) 


*+  -1  2 
^  (coth  T])+| 


^-coth"1^ 


(ri2-l) 


dn 


-HjT]Coth  ^  rj 


-  idn 

,2-i 


-i(coth_1Ti)2 


n4  ,  ..  -1  v2  1  3  -1  1 

ij—  (coth  r) )  +^r|  coth  t)+^ 


+  — )  dr]  +  7}  coth  1r|  +  £log  (  t]2-1  ) 

1  -1 


"^■(c®th  ”  1  T]  )2 


=  ■MT/+-1)(coth"lTi)2 


+  \  n(il2+3)  coth"1  rj  +  \  log( r]2-l )  + 


Similarly 


(32)  J  rj(coth  '''rj)2  dr]  =  |r(  rj2 - 1 )  (coth  %)2+r|  coth  4r)  +  -|log(T]2-l ) 

(33)  f  r)2coth  ^r)  dr)  =  coth  ^"rj  +  ^  rj2  +  g  log  (rj2-l) 


and 


(3*0 


1) 


dr) 


2 

2 


+  i  log  (  r)2 - 1 ) 


From  (29)  to  (3*0  we  have 


Q  f, 

* 


•  ■ 


/ 


. 


;  '  .)  ;Ot  ;  '  -i  >  •  ;•  ;  ■ 


-  ;  V'.:  ;  ■  •  ’  ,v 
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(35) 


V  = 


2 

KT  -j  (l-10^+9|1+)j  ■J(r^-l)(coth"1T))2+  £  r)(rl^+3)coth~1  r|+ 


(6f*2-54^-l)^-( r)2-l)  (e°th  %)2+r]  coth 


-  2(l-10^2+9^)j^  -^coth_1n  +  ^  n2 


-  4^2(l-^2)^r)  coth  ^ Tj 

L 


+  -g-  ri-loe2+9^ 

^  i 


+  g(£) 


where  g(|)  is  an  arbitrary  function  of  £, 


Differentiating  (35)  partially  with  respect  to  |  we  get 


(36)  ||  =  K2|f9ri1|'?3-5r,lt?-t6T12|- lo7l3+43-S 


(coth  1T])2  + 


r 


lOrf'’!' -6£,  r)-l8r]^^+l4r)^ 


(coth  ^rj)  +  9ti^4^“5ti^| 


|+  S'  (£) 


Comparing  (36)  with  (30)  we  have 


(37)  S' 01)  =  -4kV 


Integrating  (37)  we  have 


g($)  =  -K2^  +  ,  where  D  is  a  constant  of  integration. 


From  (35)  and  (37)  we  have 


■ 


- 


($)& 


(38) 


n2(i-ioe2+9e4)-(i^2)2 


(coth 


-  £  f  r)2(l-10£2+9£4)  +  (l-6^2+7^4) 


L-1 

COth  T]  + 


r  2 


■J-  (l-10e2+9?U)  -  n  +  D 


As  before  the  constant  of  integration  D  can  be  calculated  from  the 
condition  that  the  line-element  becomes  flat  at  infinity,  i.e. 
v  -» 0  when  r  =  0,  z  -»  co  f  corresponding  to  r|-*°°,i*  =  l.  This 
condition  is  satisfied  if  and  only  if  D  =  1  in  (38). 


Hence 


(39)  v  =  K2|i(r,2-1) 


i,2(i-ioe2+9?4)-(i-e2)2 


(coth  ^rj)2 


_  H 

2 


rf^l-lofW^)  +  (1+6|2-7?U) 


(coth  T) )  + 


r  -jj-(i-io^2+9^4)  -  n4  +  1 


or 


=  9P2^4^  ( T)2 - 1 )  (9il2-l)  (coth  ‘*‘ii)2  +  (14t]-18t]^)  coth  ^  rj  +  9t]2 

L 


+  2p2^2  (r]2-l)  (l-5r]2)  (coth  1 T))2  +  (10 rp ~6r\)coth  %  -  5r]2 


+  P 


(r)2-l)2(coth  ^r])2  -  2(r^+r))  +  r|2  +  b 


(where  |32  =  4K2) 


'  ' 

•  ».  W  M 


' 
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APPENDIX  V 

FIELD  OF  A  MASS  DIPOLE 


Below  are  the  detailed  calculations  for  the  expression  v  ,  for 
the  field  of  a  mass  dipole  (c.f.  Chapter  III,  page  52). 

If  A  is  the  Newtonian  potential  of  a  dipole,  then 


(1) 


A 


mz 

3 

P 


.  2  2  2 

where  p  =  r  +  z 


Levi-Civita  equations  are 


(2) 


dv  _  f  /dA\2  /^A\2 

dr  1 ^  dry  'dz' 


dv  _  0  dA  dA 

dz  1  dr  dz 


From  (l) 


dA 

dr 


m 


1_ 

3 


mz 


3 

P 


(3) 


m 

~3 


3mz  z 

3  P 


m 


3mz 


dA 

dz 
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Hence 


dv 


•2r  (^ES)  (E-  -  & 

P  P  '  P  ' 


/  2  2 
6m  r  z 


(4 


3z 


or 


dv 

dz 


/  2  2  -o  2  2  3 

6m  r  z  18m  r  z 

8  +  10 


6m2r2z  I8m2r2z(p2-r2)  ,  2  2  2N 

- 3—  +  — — L  (as  p  =  r  +z  ) 

P 


00 


dv  _  12m2r2z 

dz 


l8m2r4z 

10 


Integrating  we  have 


(5) 


O  2  2  _  2  4 

2m  r  9  m  r  rt  \ 

v  .  .  — ¥  — g-  +  f(r) 

p  p 


Differentiating  (5)  partially  w.  r.  to  r. 


dv 

dr 


4m2  r 


12m2r5 


2  3 
9m  r; 


l8m2r5 

10 


+  f'(r) 


(6) 


4m2  r 


21m2r5 


l8m2r'5 

10 


f'(r) 
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Also  from  (2)  and  (3) 


(7) 


9ni2r\2 

Sr  10 

P 


2 

m  r 


n  2  4 

9m  rz 

10 


6m2 rz2 
+  8~ 


Comparing  (6)  and  (7) 


2  3  2 

f(r)  =2=4^ 

P 


2 

m  r 

T~ 

p 


9m2rz^  6m2rz2 
10  +  8 


4m2r  2lm2r^  l8m2r^ 

T"  8  +  10 

p  p  p 


0  2  3  2 

f(r)  =^F- 


*  2 
3m  r 

”3" 

P 


9m2rz[|'  6m2rz2 

“ 10~  +  8"" 


21m2r5 


l8m2r5 

10 


3m  r  ,22  4  ^4  n22  22  ,4, 

— Yq  [3r  Z  +  p  -  3z  +  2p  z  -  7r  p  +  6r  ] 


3m  r  r_  2  2  ,  2  2,2  _  4  2  2N  2  _  2,  2  2,  ,  4, 

-"Yq  l3r  z  +  (r  +z  )  -  3z  +  2(r  +z  )z  -  7r  (r  +z  )+6r  J 


3m  r  ,.22  4.22  ^  4  0  2  2  0  4  7  4  7  2  2  ,  4, 

-77-  [3r  z  +r  +2r  z  +z  -3z  +2r  z  +2z  -7r  -7r  z  +6r  ] 


=  0 


Hence  from  (l)  and  (5) 


A  = 


mz 

3 


v  = 


2m2  r 2  9  m2rl+ 

~6  +  4  ~ir 
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APPENDIX  VI 


DERIVATION  OF  FIELD  OF  n-BODIES, 


Following  are  the  detailed  calculations  for  the  function  v. 


n  collinear  axially  symmetric  bodies  (c.f.  Chapter  III,  page  55)- 


In  the  case  of  a  Weyl  manifold  whose  Newtonian  analogue  is  the 

static  field  of  n  distinct  rods  EL  ,  B^,  B^,  B  ,  we  set 

l’  2’  3  n 


(1) 


n 


"  X  At 


i=l 


where 


(2) 


m.  p ,  +  p !  -  b . 

i  i  ill 

A.  =  r~  log  — ■ — - — r r~ 

i  ba  °  p.+p  +b. 

i  ill 


is  the  external  Newtonian  potential  of  B^  of  mass  kl  extending  from 
z  =  a.  -  Ab„  to  z  =  a »  +  4b .  and  thus  satisfies 


(3) 


'7% 


=  0 


(4) 


2  2 

r  +  z.  , 


2  ,2 

r  +  z  ’ 

i 


(5) 


z  „  =  z 

1 


(a  + -§b  )  ,  z !  =  z  -  (a  -  J-b  )  . 


To  obtain  v  from  (3*2)  (Chapter  III),  we  compute  from 


for 


\ 
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(6) 


- 

dr 


r&i  -  !2i  ^ 
\br  c3z  c)z J 


ax.  ax 


(7) 


*>■ 


)A.  ax.  ax.  ax 


ST 


1  ^  +  sr  sr) 


we  have 


(8) 


2  *  ai  =  ibi  ?t  1i  >  r  =  -l)a(l-^)= 


where 


pi  +  p\ 

n,-  =  — k - 1  and  £,• 


11  >  bi 


p:  -  Pt 

b. 

i 


Similarly, 

(9) 


z-aj  ■  ?bj  ^  ’ij  '  r  =  iV^ 


where 


P  .  +  p  t 
n  =  -1—1 
nj  b, 


P  “  0  . 

IJ _ -U 

b. 

J 


From  (8)  and  (2)  we  have 


(10) 


m .  r)..  - 1 

X.  =  —  log  --~r 

l  b.  n.+l 

x  'i 


also 


(ID 


ar)±  2tj1(^  -D2(i-q)2  a^  2^-1) 

ar  “  t  /  2  e2v  *  ^z  ,  /  2  42v 

VVV  MV^i* 


C^i  sEJvD'tKi^  2^(1-?®) 

^  ‘  '  bt(faf)  ’  ^  '  bt(^) 
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From  (lO)  we  have 


(12) 


S\ 

2m^ 

1 

4m 

Sr 

bi 

(7-d 

Sr 

bi 

(v^X-ii-D* 

SA. 

1 

2m. 

1 

1 

4m 

1 

4m.  £ 

xi 

Sz 

"  bi 

( v1) 

5i”  = 

b. 

1 

Vv# 

Hence 


(13) 


ax.  ax. 

_ l 

Sr  Sz 


) 


1+mi 

bi  (v^HviF 


^  i  4m .  rj .  ( 1  -i  . ) 2 

<V«?  (i^)(^-1)i 


S, 


<v«?> 


Also  from  (8)  and  (9) 


(14) 


P  JL 

(vi)2 


2r 

b.(i-|?)i 

1 


(v1)* 


2r 


b.(l-e2)2 

J  J 


Hence  from  (7),  (13)  and  (14)  we  have 
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(is) 


ax  ax  dx  ax^ 

2r (  ^ ^ 

or  dz  oz  or  J 


16111^^^  .(l-g^)  l&ryn^g^l-g2) 


(bibj)bj(i?-<?)(^j)  r  NVMv^V^ 


l6m. m . 

L 1 


b.b 

1  'i  ^i/v  'j  ^2 


41  ^ +  ^ 


l6m.  m 

1  1 


b.b  .(r)2-i2)(i]2-i2) 
1  2  'i  1  'j  j 


b. 

1 


b  . 
J 


b. 

1 


b .  b 
1-  J 


Pi^i .  z±l£  plh  !iZil 

b .  |  b .  b.  b.  b.  b.  J 

1L  j  1  - 


b  .  b  . 
J  J 


b. 

1 


16m.  m 

L L 


b.b.(n2-e2)(^2-e2) 

1  2  '!•  1  'j  J 


r  (p  . +p  !  )  (p  .“P  .  ) 
v  1  1  .1  2 

b2b. 

J *  1 


2 ( z-a . ) (p! -p . ) (p ' -p  . ) 

\  L  /  \r-x  r-x>  \r- J  r- j 

2  2 
b  bf 

1  J 


(p  ,+p (p !  -p .  ) 
J  J  1  1 

b2b  . 

1  J 


2(z-ap(p'-p.)(p^-p.) 

-  - 

bf  bf 

1  J 


b  b  m.Tn. 
p .  p  .  p !  p 


b.  b2 
1  J 


2(z-a.)(p!-p.)(p*-pj) 

2  2 
bf  bf 
1  J 


(p  +p' (pl-p^ 

■  ■  ■■mL  1  1 

b2  b. 

1  J 


2(z-a  )(p  -p  )(p!-p  ) 

J  J _ J _ l _ L_ 

2  2 
bf  bf 

1  J 


r 


■  ,.v  .  .  ...... 


i  ;■'<  3  ■  * 


"■  ■ .  -• .  ■  •  •  . .  •  - 

<! 


•  ’■ 


-  1^8  - 


dv 


y.  =  (yyV,) 


pi  pj  pi  °j 


bi^pipj"pipj+pipj_pipj  J  "2(z"ai)(PiP1"Pi,Pj-PiPj+PiPj 


+b  j  [pjp  ! -p.p  }-2(z-a^  )  {p|p^-p  .p^-pjp^+p.Pj  ) 


(m,m./b  .b  .  ) 
i  J *  1  J 

p .  p  .  p !  p  \ 

i  J  i  J 


-p.p,  {b.  +  2(z-a.)  +  b.  +  2(z-a  )] 


1  j  1 


1  J 


+p.p'.  fb.  +  2(z-a.)  -  b.  +  2(z-a.)) 
i  J  i  i  J  J 


-p!p.  [b.  -  2(z-a.)  -  b.  -  2(z-a.)} 
i  J  i  i'  j  yJ 


+p!p'.  {b.  -  2(z-a.)  +  b.  -  2(z-a.)) 
i  J  i  i  J  J 


(m.m  ,/b  .  b  .  )  r 

p^pj-pf-^  |/ptpj(2zi  +  2zP  +  pipjt£zi  +  2zj] 


■p!p.{-2z.  -  2z[ )  +  p.'p*.  {-2z.  -  2z  . ) 
1  1  1  J  1  J  1  j 


d  v.  . 

<l6> 


2m. m.  r  (z* +z I ) 

x  j  i  y 

b  b  ”  p »  p ' 

i  j  L  i  j 


( z !  +z  . ) 

i  y 

p !  p . 

i  J 


( z .  +z '. ) 
i  J 

p .  p  I 
i  J 


(zi+z5) 

Pi  Pj 


By  straightforward  but  tedious  calculations  using  (6),  (8)  and 
(9),  we  find  that 


dv.  .  2m. m. 

(17)  Z  -JJ.  =  -1-I 

Ur;  8r  b.b. 


i  J 


z ' z  1  -r 

XJ  _ 

rpiPj 


z ! z  .  -r 
i  J 

rp !  p  . 
i  J 


z .  z -r 

rp .  P ! 

i  J 


z , z  .  -r 

_LJ _ 

rp  .p  . 
i  J 
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In  order  to  get  the  value  of  v_  from  (16)  and  (I7),  let  us 
establish  the  following  identity 


(18) 


r^-z ! z  . 

d  log[r"1E(i',  j)]  =  J  dr  + 

r 


z!+z  . 
- *1  dz 

pipj 


Set 


-1, 


(19)  f  =  log[r~  E (i ’ , j  )  ]  =  log 


-p !  p  .  +z !  z  .  +r 


then 


M 

dr 


rr 


(p!p.+z!z„+r  ) 
'  i  j  i  j 


8p  dpj  2 

K  5T  +  pjgr~  +  gr }-(p'p.«'2.+r  ) 


(p!p.+z!z„+r  ) 
i  j  i  j 


r[pi  +  pj  FT  +  grHpiPj+zjZj+r1  ) 


1 


(p!p.+z!z.+r  ) 
i  j  1  j 


22222  22  2 

r  pi  +r  p.+2r  pip. -pi  p.- zlz.plp.-r  pip. 

rp  I  p  . 
ij 


(p!p.+z!z.+r  ) 
'  i  j  i  j  ' 


2/2  ,2\  2/2  2\  2.  /  2  ,  2  \  /  2  2  \  ,  , 
r  r  (r  +Zl  )+r  +Zy+T  pjpj~(r  +21  )(r  +z  ^zizjplpj' 


r  p !  p  . 
i  J 


(plp.+zlz.+r  ) 
'  i  j  i  j  ' 


k  2  ,  ,22 

r  +  r  p  p  -z  z  -  z  z  p  p 
-1-  J. . 1  -J 

r  p !  p  . 
i  J 


■  ,  \ q 


' 

'  >  ■  •••  '  I-  •  ..  . . .... 
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df 

£>r 


(p!p  .+z.'z  .+r  ) 
i  J  r  j 


(rg-2.z,)(Plp,  +z|Zj  +r2) 


r  p !  p  . 
J-  J 


i.e. 


r  -  z .  z 


Li 


r  Pi  Pj 


(20)  log[r_1E(i 1 , j ) ] 


r  -  z !  z  . 
LJ. 

r  p  !  p  . 
i  J 


Again  from  (19)  we  have 


df 

31 


( p !  p  .  +z !  z  .  +r  ) 
i  J  i  j 


3p.  dPi 
r  rfo'  ——-l  +  1 

rnpi  dz  +  J  dz 


dz .  dz ! 

I  <7  I  I  <7  m.m  III  1 

+  2i  dz  +  2,i  dz  J 


(p!p„+z!z„+r  ) 
'  i  j  i  j  ' 


Pi 


z  . 


lBj 


+  P  .  — T  +  Z  •'  +  Z  . 

J  1  J 


(p!p.+z!z.+r  )4 
i  J  i  J 


2  2 
r  p !  z .  +  z !  p  „  +  z !  p !  p  .  +zp!p.  ~\ 

-1 — J- _ LJ _ JL±-J _ 3. ■  ■L-  J 


p!  P . 
i  J 


( p !  p  .  +z !  z  .  +r  ) 
3-  J  1  J 


/2,2s  ,/22\  .  ,  , 

(r  +2;  )z  +  zj(r  +Zj)  +2jpip.  +Z.pjp. 


p !  p . 

i  J 
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(P:pi+z-zi+r  ) 


2 

r  ^zi'fzp,fztz.1^zl't,i!.l),fpip1(zi+zl) 


P1  PJ 


or 


i.e . 


df 


(pip .+z,'z  ,+r  ) 
1  j  l  j 


(2i  +  2j)(pipj  +zi2i  +r2) 


p!  P  • 
i  J 


(21)  log[r”1E(i 1 , j )  ]  = 


z  I  +  z  . 
-L__ _ 1 

pi  p . 

1  J 


Hence  the  identity  (l8)  is  proved. 


Exactly  the  same  way  we  can  show  that 


d  log[r“1E(i, j  ’ )]  = 


r  -  z .  z  \ 
1  J. 


r  p  .p  . 
1  J 


t  dr  + 


*i +  zj 

pipj 


dz 


(22) 


d  log[r'1E(i, j)] 


r  -  z .  z  .  z .  +  z  . 

-i-J.  dr  +  -i— — 1  dz 


r  p.p  . 
1  J 


p.p  . 
J-  3 


d  log[r-1E  (i  ' ,  j  '  )  ]  =  '  dr 


r  pip! 
1  J 


z 1  +  z  1 

+  -V— idz 


PiPj 


From  (16)  to  (22)  we  have 


2m . m .  r 

(23)  2,  dv  =  ^ 

1  J 


d  log{r  1E(i',j))+d  log (r"1E(i,  j ' ) }-d  log{r'1E(i, j ) } 

-d  log {r_1E (i ' , j ' ) ] 


: 

•  ^  -v1.:, .  a. 


. 


•'  *},W’  M  {  '•  .  a)  ^  k?  i  ^  £(f 
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Hence 


(24)  ^  ViJ 


2m  m . 

bibj 
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APPENDIX  VII 

LIMITING  PROCESS  OF  SCHWARZSCHILD  PARTICLES. 


In  this  appendix  we  show  that  when  each  ->  0  ,  the 

(3.61)  is  formally  the  field  of  n  collinear  Curzon  particles 

III,  p.55). 


(1) 


From  equation  (3*56),  Chapter  III,  we  have 


log 


P.  +  p! 

1  1 


p.  +  p: 

1  X 


where 


(2) 


2  2  ,2  2  ,2 
r  +  z.  ,  p!  =  r  +  z]  ,  m.  is 
11  1  }  1 


(3) 


Zo  =  z  -  (a.  +ib„),  zl  =  z  -  (a.  -  ±b 
1  ^  1  *  1  V  ^ 


In  (i)  set 


(4) 


f  (p. 


f  p  i ^  i 


p  „  +  p  .  -  b  „ 

=  m.  log  -■  ■; —  7 

1  pi  +  pi  +  bi 


and 


(5)  gO^)  =  bi 

When  b  =>0  ,  f(p  ,p!,b  )  -*  0  ,  g(b  )  ->  0  .  So  V 

JL  III  1.  J 

indeterminate  when  b„  =  0  . 

1 


solution 
c.f.  Chapter 


constant 

)  * 


is 
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Using  L'Hospital's  rule, 


(6) 


df_ 

db, 


=  m 


1 


dp,  dp!  dp,  dp 

+b  (pi'1<)i+bi)t3b7  +  SbT  -l)-(Pi-H;.i-bi)(Sr+  +1} 


(p.  +  p^  +  b  ±Y 


From  (2)  and  (3) 


(7) 


dp 
db . 


20.  ’ 


dp^ 

db_^ 


z  l 


From  (7)  and  (6),  we  have 


df  1  f  ~z±  zl 

db.  "  ™i  (p.+pl-b.  )(p.+p!+b.  )  "I  ^pi+pi+bi  ^2p .  +2p!  ^ 

1  1  1  i/v  1  1  x'  ^  11 


-(p.+pl-b.)^  +  2^  +  1) 


m.  r  zx  z.  'j 

(p.+p! -b.  )(p.+p!+b.  )  bi  ^p]"  "  p .  ^  “  2 ^Pi  +  Pi  j 
v  1  1  l  i  1  x'  ^  i  1  J 


df 

limit  rr — 

,  ob . 

b  o  1 

1 


2m.  (pi) 

. 7T - -  (since  p !  ->  p .  when  b.  ->o) 

rp  1  1  1 

2pf 

1 


Also 


m. 

1 

Pi 


(9) 


limit 
b .  —*  ( 


dg 

*bi 


=  1 
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(10) 

which  is 

(ID 

where 

(12) 

and 

(13) 

In  (11) 
and  v .  . 


From  (8)  and  (9) 


m. 


the  Newtonian  potential  for  n  "Curzon  particles”. 


Also  from  equation  (3.61),  Chapter  III,  we  have 


v.  .  = 

ij 


m.m. 

b.b. 
1  J 


m.m , 


log  HU* 


oJ- 


E(i, j)E(i’, j 


[1°8  E(i',j)+log  E(  i,  j  '  )-log  E(i,j)-log  E(i,,j')] 
i  j 


E(i',j)  = 
E(i,j')  = 
E(i,j) 
E(i',j')  = 


p !  p  .  +  z .'  z  .  +  r 
1  J  1  J 


p ,  p  !  +  z .  z  !  +  r 
1  J  1  J 


p.p ,  +  z.z ,  +  r 
1  J  1  J 


p  1  p  1  +  z  *  z  1  +  r 
1  j  1  J 


In  (11),  set 


9(p,z,b)  = 


m.m.  log 


lmj  iOS  E(i,j)E(l"r7 


b)  =  b.b.  . 
1  J 


if  we  keep  b^  fixed  and  let  b^  ->  0  ,  then  cp  ->  o  and  if/  ->  0 

indeterminate.  As  before  using  L’Hospital's  rule  of  differentiation 


' 

i  riu  ■  . 


S.)  • 


146  - 


we  have 


(14) 


dq> 


^  =  mimj 


E(l’’J)  +WTF)\  E(l'j,) 


-  KITjT  \  E(1’J)  -  i^TrT  \E(1’'J,). 


From  (7)  and  (12)  we  have 


(15) 


i 


\  E(i,j) 
i 


p  .  z +p !  z  . 

J  1  ^  e(i  *»)  = 

2p !  9  °b.  ' 

1  1 


p .  z  . +p . z . 
..JLJ.  -J.  ,1 


2p, 


P3zlWiZ3 

2pi 


5.  E(i ' ,3 ' )  = 


p  \  z !  +p !  z  1 
J.  1 

2PJ 


Similarly, 


\  E(i'j) 
j 


p !  z  . +z ! p  . 

1  i.  1  ]  ;  8,  e(i,3') 


2p3 


b  . 
J 


p .  z  \ +z . p  ! 
l.J . 1  J. 

2p3 


(16) 


6b  E(l,3)  = 


plz3^3zl  . 
2p3  ’ 


\  E(i',j') 
j 


p!z!+p:z: 
i  1  J  1 

2P3 


From  (14)  and  (15) 


(IT) 


1 .  . .  69 

limit  ^rr — 
,  db . 

b  o  i 

1 


m.m . 

-LJ. 

pi 


eTITTT  (p3zi+z3pi)-  e(T7F7  (pjzi^i-j) 


(since  b±  -*  o  ,  E(i',j)  =  E(i,j)  ;  E(i,j')  =  E(i',J')) 


Also 


(18) 


limit 

b^  —*■  o  1 


=  b 
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Hence, 

(19) 


C*P 

^bT 

limit  v, .  =  limit  ^ , 1 

iJ  V°  1 


m.m 

bjpi 


(pizi+eizi) 

e[i7)T 


(piZi'fpiZl) 

"E(I7jV~ 


In  (19)  if  b j  ->  o  ,  then  v_  again  becomes  indeterminate. 
Hence  ,  differentiating  the  numerator  and  denominator  with  respect 


to  b .  we  have 
J 


v.  .  = 


m.m . 

LJ. 


2  2 

— r-TT  [-(p.p .+z. z .) (p.p ,+z.z ,+r  )+(z.p ,+z .p. )  ] 
ll  /  2\2  ll  1  li  J  1  J  11  1  1 

J  p  .  p  .  ( p  .  p  .  +z .  z  .  +r  )  J  J  J  J  JJ 


1  J  1  J  1  J 


m  m . 


pipj 

(p.p ,+z.z ,+r2) 

1  J  1  J 

m.m . 

1  1 

2  2 
r  (-p.p ,+z.z ,+r  ) 

1  l.  1  1 

P.P. 

1  J 

2  2 

(p.p.+z.z.+r  )(r  +z.z.-p.p.) 
1  J  1  J  1  J  1  J 

m.m. 

1  J 

2  /  2  \ 
r  (r  +z.z .-p.p . ) 

1  j  1  y 

p.p. 

1  J 

/  2  s2  2  2 

(r  +z±z.)  -PiPj 

m.m. 

1  J 

2  /  2  v 

r  (r  +z.z  .-p.p  .) 

1  .1  1  r 

P.P  • 

1  J 

/  2  \2  /  2  2\ /  2  2\ 

(r  +z.z.)  -(r  +z. )(r  +z . ) 

1  J  i  r 

1 

(r2+z. z .-p .p . ) 

1  J  1  J 

p .  P .  / 

1  j  (z . -z  . ) 

i  J 


. 
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APPENDIX  VIII 

DETERMINATION  OF  v  BETWEEN  THE 

TWO  RODS  ON  THE  AXIS  OF  SYMMETRY, 


.  th 


th 

In  this  appendix  we  determine  the  value  of  v  between  i  and 
rods  on  the  axis  of  symmetry  (c.f.  Chapter  III,  page  58). 


From  equation  (3.6 l),  we  have 


(1) 


where 

(2) 


and 


(3) 


m.m . 

v  1  3 

log 

E(i'» j)E(i. j ' ) 

ij  b.b. 

1  J 

E(i,j)E(i,,j* ) 

m.m 

2  2 

[r  +z !  z  .  +p !  p  .  ]  [  r  +z.z!+p.p'] 
1  J  1  J  1  J  1  j 

 1  J 

log 

b0b# 
1  j 

2  2 
[r  +z.z  ,+p.p  .  ][r  +z!z'+p!p!] 
1  J  1  j  i  j  i  J 

2  2 

2  . 

,2  2  ,2 

+ 

t-i 

n 

Q. 

Zi  ' 

*1  =  r  +  1 

zi  =  z  "  (ai  +  >  zi  =  z  -  (ai  “  ibi) 


z.  =  z  -  (a.  +  ^b.)  ;  z '.  =  z  -  (a .  -  ib , ) 
J  J  2  J  ’  3  3  2  j 


When  z  <  a^  -  -|-b^  on  the  axis  of  symmetry, 
z\  are  all  negative  and 


z , 
1 


z !  ;  z ,  and 
1  '  J 


(4) 


p!p.  =  z.'z.  ,*  p.pl  =  z.z!  :  p.p.  =  z.z,  and  pjp!  =  zjz! 

1  J  1  j  *  *iHj  i  j  >  HiHj  1  j  KiHj  i  j 


I'- 
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Hence  from  (l)  and  (4) 

vij  =  0  (r  =  0,  i  <  j  ,  z  <  (ai-ibi))  . 

Similarly,  it  can  be  shown  that  when  z  >  a,  +  Ab,  then  z , ,  z!  ,  z.  and 
'  J  2  j  i’  i  '  J 

Zj  are  all  positive  and 

vij  = 0  (r  =  °»  1  s  j>  z  >  aj +  • 


th 

Let  us  now  consider  the  case  when  z  lies  between  the  i  and 


.th 

j  rods. 


th  th 

For  all  values  of  z  between  the  i  and  the  j  rods  on  the 


z-axis,  z ^  and  z!^  are  always  positive,  while  z^  and  zj  are  both 
negative,  v  does  not  vanii 
the  z-axis „  From  (l)  and  (2) 


negative,  does  not  vanish  in  between  the  i*^  and  the  j*'*1  rods  on 


(5) 


(r^+zjz  .+  pip  . )  =  [r2+z„*z.  +  \f  (r^+z  !*)  (r^+z'f)  ] 
i  j  r  j7  i  j  i  2 


r  2  ,  /T2"2  2/  ,2  27  4-, 

I  r  +z  l  z  .  +  v  z  z  .  +  r  ( z !  +z  .)  +  r  s 
i  3  i  j  i  j7 


/  ,2  2v 

2  (zi  +V 

r^+z  J  z  .  -  z !  z  .  -H+r  - —  + 

ii  ill  ,22 

J  J  z .  z .  z ;  z . 

1  J  1  J 


.2  2/ 


O  ,2  2 

2  ,  ,  f,  r2  zi  +Z1 

r  +zizj-zizj  |1+  -  -j2 

L  L  z !  4Z  . 

1  J 

2 

(neglecting  higher  powers  of  r  ) 

2  (z!  -  z  .  )2 

r  v  i  i7 

2  z 1  z  * 

i  j 


■+■  •  • 


. 
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Similarly  it  can  be  shown  that 


(r  +Vj+plPj)  " 


2  (z  -z \Y 
.  £_  1  3 


2  « 


(6) 


(r  +  +  PjPj)  = 


2  (y»,) 


‘i  *j 


(r  +  +  pjpj)  = 


2  (z!-z!)‘ 

r_  x  j' 


*i  j 


Hence  from  (l),  (5)  and  (6) 


(7) 


m.m. 

1  j 

i]  b  .b  . 
i  J 


log 


(zj-z.)2(z.~z!)2 

A...  J.L  a_jl__ 

(z. -z  .  )2(z! =  z 1 )2 
x  J7  1  j 


2m„m , 

b.b. 
i  J 


log 


(z!“Z.)(z.-z!) 

A  ll  -A  .J... 

(z. -z  ) (z! -z '. ) 
v  i  j'v  i  J 


(r=0,  i  <  j,  a.+ib.  <  z  <  a.-A-b.) 
\  >  Jf  x  2  x  J  2  J 
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APPENDIX  IX 


EQUIVALENCE  OF  v  AND  THE  NEWTONIAN 


FORCE  BETWEEN  THE  RODS. 


In  this  appendix  we  show  the  equivalence  of  the  Newtonian  force 
between  the  itb  and  jtb  collinear  rods  separated  by  the  region  of  the 
z-axis  in  the  Euclidean  map  and  the  relativistic  expression  v  in 
between  these  two  rods  on  the  axis  of  symmetry  (c.f.  Chapter  III,  page  58) • 

th 

Let  us  consider  any  small  length  dx^  on  the  i  rod  of 

length  b^  and  of  uniform  mass  nu  .  Let  the  distance  of  this  length  from 

the  centre  of  the  rod  be  x^  .  Similarly  dx^  is  any  small  length  on  the 
th 

i  rod  of  length  b.  and  of  mass  uniform  m.  .  x.  denotes  the  distance 

J  3  3 

th 

of  dx,  from  the  centre  of  the  i  rod. 

3 

If  F„ .  is  the  force  acting  between  the  two  rods,  then 
ij 


(1) 


ij 


m  m. 

b  ,b  . 
1  J 


^b: 

4b, 


-H=i-b 


dx. dx . 
1_ 3 


(r .  .  +  x .  -  x.  )' 
v  ij  J  1 


i  <  j 


where  r..  is  the  distance  between  the  centres  of  the  two  rods, 
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which  is  the  same  result  as  we  have  obtained  for  the  expression 


ij 


tVi  th 

between  the  two  i  and  j  rod  in  Appendix  VIII,  equation  (7)  . 
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APPENDIX  X 

VERIFICATION  OF  EQUATION  (5.18). 

Using  the  Kruskal's  formula  for  the  co-ordinate  transformation, 
following  are  the  detailed  calculations  for  the  metric  transformation 
(5.18)  (c.f.  Chapter  V,  page  97)  . 

From  equation  (5. 15)  Chapter  V,  we  have 
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Defining  p,  h,  £  and  t  by  the  Kruskal's  formula  (5.4)  we 
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APPENDIX  XI 

DETERMINATION  OF  K  OF  EQUATION  (5-26). 


Below  are  the  detailed  calculations  for  the  determination  of 
equation  (5.26)  (c.f.  Chapter  V,  page  99)* 

From  equation  (5. 13),  we  have 
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Also  from  equation  (3),  Appendix  VIII  we  have 
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Hence  from  equation  (l),  Appendix  VIII 
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For  any  value  of  z  between  (a  -  4b  )  and  (a  +  ib  ) 
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APPENDIX  XII 

ELLEME NTARY  FLATNESS  CONDITION  FOR  BONDI  U-DIPOLES 

Below  are  given  the  detailed  calculations  for  the  determination 
of  elementary  flatness  condition  for  Bondi  4-dipoles. 

Relabelling  the  three  particles  as  0,  1,  2  instead  of  1,  2,  3 
and  setting  mQ  =  =  -a^  -*  oo  f  the  equation  (3»70)  of  Chapter  III, 

page  59  takes  the  form  as: 
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In  the  limit  as 
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